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Abstract 


The  heterodyne  /  direct  detection  DIAL  comparison  (HD/DD  DC)  experiment  series  was 
conducted  at  Kirtland  AFB,  NM,  to  simultaneously  characterize  and  compare  the  radiometric 
and  chemical  detection  sensitivities  of  heterodyne  and  direct  detection  DIAL  systems.  The 
system  developed  by  the  Air  Force  Research  Laboratory  Directed  Energy  Directorate 
demonstrated  the  first  known  programmable  and  shot-to-shot  wavelength-agile  heterodyne 
DIAL  measurements.  The  experiments  studied  radiometric  issues,  speckle  mitigation  through 
spread  spectrum  (modelocked)  operation,  and  chemical  detection  sensitivities.  The 
measurements  were  performed  over  horizontal  paths  at  standoff  ranges  from  4  to  15  km,  using 
both  natural  and  man-made  targets.  Heterodyne  and  direct  detection  radiometric  and 
chemometric  results  are  presented  and  contrasted,  and  are  compared  with  predictions  from 
simulations  and  models. 
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1.  Introduction 

A.  Overview 

The  Heterodyne  /  Direct  Detection  DIAL  Comparison  (HD/DD  DC)  experiments  were  designed 
to  provide  a  direct,  simultaneous  comparison  of  the  radiometric  and  chemical  detection 
sensitivities  of  the  two  receiver  techniques.  From  October  1998  to  September  2003,  the 
HD/DD  DC  measurements  were  a  main  component  of  the  Laser  Remote  Optical  Sensing 
(LROS)  program  at  the  Air  Force  Research  Laboratory  Directed  Energy  directorate  (AFRL/DE). 
In  general,  it  is  expected  that  heterodyne  detection  will  have  a  significant  radiometric 
sensitivity  advantage  (by  a  factor  of  approximately  102  to  103)  over  direct  detection.1  A  key 
issue  for  heterodyne  detection,  however,  is  the  sensitivity  of  the  process  to  speckle,  which 
affects  the  chemical  detection  sensitivity.  For  a  direct  detection  system,  it  is  possible  to 
spatially  average  multiple  speckles  across  the  receiver  aperture,  thereby  reducing  the  speckle- 
induced  fluctuation  in  the  received  signal.  For  a  common  single-mode  heterodyne  receiver, 
the  requirement  of  having  a  coherent  phase  front  over  the  receiver  aperture  results  in  a 
speckle-limited  single  shot  SNR  of  unity.  For  the  HD/DD  DC  experiments,  a  speckle  mitigation 
technique  was  studied  which  uses  a  spread  spectrum  (modelocked)  transmit  waveform.  In 
this  technique,  multiple  laser  longitudinal  modes  are  transmitted,  each  of  which  is  speckled 
independently  when  returned  from  a  diffuse  target  with  sufficient  longitudinal  extent.  An 
illustration  of  the  SNR  characteristics  versus  range  for  each  of  these  techniques  is  given  in 
Figure  1.  The  values  in  the  Figure  were  not  calculated  for  specific  systems,  but  are 
representative  of  heterodyne  and  direct  detection  systems  with  similar  transmit  energies  and 
receiver  apertures.  A  common  characteristic  of  all  of  the  techniques  is  that  at  short  ranges  the 
SNR  will  be  speckle  limited,  and  therefore  mainly  independent  of  range.  At  some  distance,  the 
SNR  will  become  signal  limited,  and  will  decrease  proportionally  with  the  range  squared  (R2), 
or  faster.  The  heterodyne  radiometric  advantage  is  that  the  breakpoint  between  speckle- 
limited  and  signal-limited  operation  will  occur  at  significantly  greater  ranges.  The  SNR 
advantage  at  short  ranges  in  the  speckle-limited  regime  seen  for  a  direct  detection  system  in 
comparison  to  a  single-mode  heterodyne  system  can  be  reduced  or  nearly  eliminated  by 
utilizing  a  spread  spectrum  transmit  waveform.  The  cost  of  this  improvement  in  comparison 
with  the  more  common  single-mode  heterodyne  systems  is  slightly  more  complexity  in  the 
front  end  of  the  heterodyne  receiver,  where  a  higher  bandwidth  capacity  is  required.  This 
bandwidth  can  later  be  reduced  by  either  analog  or  digital  processing  to  a  level  similar  to  that 
of  a  single-mode  heterodyne  receiver.  The  complexity  of  the  spread  spectrum  (modelocked) 
transmitter  actually  tends  to  be  less  than  that  of  a  wavelength-agile  single-mode  transmitter. 
At  short  ranges  a  direct  detection  system  will  always  have  an  advantage  over  a  heterodyne 
system,  since  direct  detection  can  utilize  both  spatial  averaging  and  spread  spectrum 
operation  to  provide  speckle  mitigation,  while  heterodyne  detection  can  utilize  spread 
spectrum  operation  but  not  spatial  averaging.  The  heterodyne  system  advantage  will  always 
occur  for  longer  range  operation,  where  heterodyning  can  continue  to  make  measurements 
after  the  direct  detection  system  return  level  has  dropped  below  the  system  noise  floor. 
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Figure  1.  Illustration  of  spread  spectrum  heterodyne,  single-mode  heterodyne,  and 
direct  detection  SNRs  as  a  function  of  range. 


B.  Development  Timeline 

The  Heterodyne  /  Direct  Detection  DIAL  Comparison  (HD/DD  DC)  experiments  utilized 
components  from  two  separate  systems,  the  Laser  Airborne  Remote  Sensing  (LARS)  and 
Coherent  Remote  Optical  Sensing  System  (CROSS).  The  LARS  program  demonstrated  the 
capability  of  direct  detection  DIAL  measurements  from  an  airborne  platform,  achieving 
multiple  chemical  detection  at  a  one-way  standoff  slant  range  of  greater  than  30  km  (see 
Figure  2)2.  As  the  LARS  program  was  concluding,  the  CROSS  efforts  began  to  develop  a 
heterodyne  DIAL  capability,  with  the  eventual  goal  of  extending  the  operational  range  to  a 
one-way  standoff  slant  range  of  greater  than  80  km.  The  development  timeline  for  the  LARS 
and  CROSS  systems  is  shown  in  Figure  3.  The  primary  development  required  for  the  CROSS 
program  was  a  wavelength  agile,  stable  C02  local  oscillator  capable  of  high  pulse  repetition 
frequency  (PRF)  operation  (with  an  eventual  PRF  goal  of  greater  than  1  kHz).  The 
Wavelength  Agile  Local  Oscillator  (WALO)  was  developed  to  meet  the  goals  of  the  CROSS 
system3,  and  is  discussed  in  more  detail  later  in  this  report. 
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30  km  Detection  of  Multiple  Chemical 
Plume  from  Airborne  Platform 


LARS  Established  the  Current 
Record  for  Standoff  Range  with 
an  Airborne  Chemical  Detection 
Sensor 
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Figure  2.Airborne  chemical  detection  at  30  km  slant  range  using  LARS  direct  detection  system. 
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2.  Experimental  Information  and  Configuration 

A.  HD/DD  DC  Equipment  Description 

The  layout  of  the  equipment  used  in  the  HD/DD  DC  experiments  is  shown  in  Figure  4.  The 
various  subcomponents  will  be  described  in  more  detail  in  the  following  sections.  A  common 
transmitter  was  used  as  the  illumination  source  for  both  the  heterodyne  and  direct  detection 
receivers.  The  transmitted  light  was  passed  through  a  gas  absorption  cell,  and  the  optical 
return  signal  for  both  receivers  also  passed  through  the  gas  cell. 
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Figure  4.  Layout  of  transmitter,  direct  detection  receiver,  and  heterodyne  receiver  for 
Heterodyne  /  Direct  Detection  DIAL  Comparison  (HD/DD  DC)  experiments. 


B.  Transmitter  System 

The  Breadboard  Oscillator  (BBO)  C02  laser  developed  by  Textron  Systems  Corporation  was 
utilized  as  the  transmitter  for  the  Heterodyne  /  Direction  Detection  DIAL  Comparison  (HD/DD 
DC)  experiments.  The  BBO  laser  was  also  used  in  the  Laser  Airborne  Remote  Sensing  (LARS) 
program,  which  developed  an  airborne  direct  detection  DIAL  system  that  demonstrated 
simultaneous  detection  of  multiple  chemicals  at  ranges  of  greater  than  30  km  (see  Figure  2). 
The  laser  and  associated  cavity  optics  are  shown  mounted  on  the  left-hand  side  of  the  LARS 
airborne  optical  bench  in  Figure  5.  The  BBO  laser  is  an  RF-excited  transverse  electric 
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atmospheric  (TEA)  C02  laser  with  output  pulse  energies  of  greater  than  4  J  on  strong  C02 
transition  lines,  and  a  pulse  repetition  frequency  (PRF)  of  up  to  30  Hz.  For  the  HD/DD  DC 
experiments  the  laser  PRF  was  limited  to  approximately  1  Hz  by  the  capabilities  of  the  data 
acquisition  and  control  system.  The  laser  uses  a  mechanical  grating  to  change  the  output 
wavelength  in  a  programmed  sequence.  The  HD/DD  DC  experiments  commonly  utilized 
sequences  consisting  of  12  to  13  of  the  approximately  60  laser  lines  accessible  by  the  BBO 
using  the  12C1602  isotope  in  the  laser  gain  medium.  To  investigate  speckle  mitigation  the  laser 
operation  was  often  alternated  between  non-modelocked  and  modelocked  configurations  for 
successive  data  sets.  In  both  configurations  the  laser  had  multiple  longitudinal  mode  output, 
but  modelocked  operation  resulted  in  more  longitudinal  modes,  with  higher  and  more  stable 
amplitudes.  The  key  parameters  for  the  BBO  laser  are  given  in  Table  1. 


Figure  5.  Photograph  of  BBO  laser  and  system  on  Argus  flight  bench. 


C.  Direct  Detection  Receiver 

A  diagram  of  the  LARS  bench,  comprising  both  the  laser  transmitter  and  the  direct  detection 
receiver,  is  shown  in  Figure  6.  The  transmit  optical  path  is  shown  in  blue,  and  is  primarily  on 
the  left  side  of  the  optical  bench.  The  direct  detection  receive  signal  is  shown  in  purple  on  the 
right  side  of  the  optical  bench.  The  direct  detection  receiver  utilizes  a  16”  diameter  reflective 
telescope,  with  a  secondary  /  pickoff  mirror  obscuration.  The  receive  light  is  focused  onto  a 
0.500  mm  diameter  liquid  nitrogen  cooled  HgCdTe  detector  with  a  bandwidth  of  80  MHz.  The 
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signal  from  the  detector  is  sampled  at  60  MS  a/s,  and  stored  by  the  LARS  Acquisition  and 
Processing  System  (LAPS)  for  use  in  both  quick-look  and  offline  analysis. 


Table  1.  BBO  Laser  Key  Parameters 


Parameter 

Value 

C02  Isotope 

12c16o2 

13c16o2 

Pulse  Energy 

>  4  J  (strong  line,  using  12C1602) 

>  50  mJ  (weak  line) 

Pulse  Repetition  Frequency 

30  Hz  (maximum) 

10  Hz  (nominal) 

1  Hz  (HD-DD  DC  experiments) 

Wavelength  Tuning  Capability 

Single  shot 

Number  of  Accessible  Wavelengths 

~  60  (12C1602  isotope) 

~  50  (13C1602  isotope) 

Pulse  Type 

Gain  switch  spike  with  relaxation  tail 

Pulse  Width 

~  10  \is  (non-modelocked) 

~  6  qs  (modelocked) 

Beam  Divergence 

~  1.0  mrad  (from  laser) 

0.285  mrad  (transmitted  downrange) 

Figure  6.  Plan  view  of  laser  transmitter  and  direct  detection  receiver. 
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D.  Wavelength  Agile  Local  Oscillator  (WALO) 

The  Wavelength  Agile  Local  Oscillator  (WALO)  was  developed  for  stable  wavelength  agile 
operation  at  high  PRFs  (>  20  Hz).  The  WALO  uses  a  cavity  waveguide  for  the  C02  gain 
medium,  with  acousto-optic  modulators  (AOMs)  inside  the  laser  cavity  for  wavelength  tuning. 
Prior  wavelength  agile  C02  local  oscillators  had  used  mechanical  tuning  of  diffraction 
gratings,  which  tended  to  be  unstable  at  high  PRF  operation  (>  20  Hz),  due  to  the  mechanical 
settling  time  of  the  moving  grating  fixture.  The  WALO  utilized  paired  AOMs  to  achieve  high 
PRF,  stable  wavelength  agile  operation.  Since  the  AOMs  utilize  acoustic  waves  to  change  the 
wavelength  of  light  transmitted,  the  tuning  time  limit  was  determined  by  the  propagation  time 
of  the  acoustic  wave  in  the  transverse  direction  across  the  AOM,  and  the  stabilization  time  of 
the  C02  laser  into  CW  operation.  Tuning  rates  of  500  Hz  were  demonstrated  with  the  WALO, 
which  was  well  beyond  the  tuning  capability  of  the  transmitter  laser  used  in  the  HD/DD  DC 
experiments.  Tuning  rates  of  many  kHz  have  been  demonstrated  by  other  research  groups 
utilizing  AOMs  with  C02  waveguide  lasers,  so  the  tuning  rates  used  for  the  WALO  are  not  at 
the  limits  for  this  type  of  technology.  Since  there  are  no  moving  parts  in  the  design,  the 
WALO  is  relatively  robust,  and  well  suited  to  operation  in  both  field  and  airborne 
environments. 


E.  Heterodyne  Receiver 

The  CROSS  heterodyne  receiver  is  shown  in  Figure  7.  The  receive  signal  and  WALO  beams 
are  combined  using  a  beam  splitter,  and  are  focused  onto  a  0.100  mm  x  0.100  mm  square 
HgCdTe  detector  with  a  bandwidth  of  800  MHz. 
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Figure  7.  Layout  of  heterodyne  receiver  with  Wavelength  Agile  Local  Oscillator 
(WALO). 

F.  Gas  Cell 

The  gas  cell  used  in  the  HD/DD  DC  experiments  (see  Figure  8)  is  a  4’  x  4’  x  4’  aluminum  cube, 
with  IR-transmissive  polyethylene  windows  on  two  sides.  The  gas  cell  was  designed  for  both 
gas  and  liquid  chemical  insertion.  Gas  valves  are  incorporated  into  the  side  panels,  and  a 
funnel  mechanism  with  heated  Petry  dish  is  mounted  on  the  top  panel  to  evaporate  chemicals 
inserted  in  liquid  form.  Fans  are  mounted  inside  the  cell  to  provide  uniform  mixing  of  the 
chemicals  inserted.  The  cell  has  been  used  with  gaseous  sulfur  hexafluoride  (SF6),  gaseous 
ammonia  (NH3),  and  liquid  ammonium  hydroxide  (NH4OH).  For  the  HD/DD  DC  experiments, 
only  the  chemicals  in  gaseous  form  (SF6  and  NH3)  were  used.  Insertion  of  gaseous  chemicals 
was  performed  by  filling  sample  cylinders  with  pure  gas  to  above  the  ambient  atmospheric 
pressure.  The  sample  cylinders  were  connected  to  valves  on  the  side  of  the  gas  cell  (note 
sample  cylinder  attached  to  valve  on  right  side  of  cell  in  Figure  8  photo),  which  were  opened 
at  the  desired  point  in  the  experiment  to  insert  the  excess  pressure  from  the  sample  cylinders 
into  the  cell.  After  a  few  seconds  the  valves  were  again  closed,  to  prevent  all  of  the  gas  in  the 
sample  cylinders  from  leaking  into  the  cell.  Although  this  insertion  procedure  is  not  extremely 
precise,  it  provides  a  relatively  simple  method  to  achieve  reasonably  predictable  and 
repeatable  results.  From  viewing  the  DIAL  absorption  plots,  it  appears  that  the  internal  cell 
fans  cause  complete  mixing  of  the  chemicals  in  the  cell  within  10  sec.  The  quick- release  door 
panels  were  removed  to  rapidly  evacuate  the  inserted  chemicals  between  data  sets. 
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Figure  8.  Photograph  of  4’  x  4’  x  4’  gas  cell. 


G.  Target  Sites 

Target  sites  at  one-way  ranges  of  4.320  km  (4KS),  7.460  km  (7KS),  and  14.9  /  15.3  km  (15KS) 
were  developed  for  use  in  the  HD/DD  DC  experiments.  Other  target  sites  were  also  available 
at  approximately  2  km  (2KS),  5  km  (5KS),  and  6  km  (6KS),  but  were  not  used  in  the  HD/DD  DC 
experiments,  other  than  for  alignment  purposes.  Figure  9  shows  the  beam  paths  from  the 
transmitter  location  in  Building  770  (B770)  on  Kirtland  AFB,  to  the  various  target  sites.  The 
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Figure  9.  Topographic  map  with  beam  paths  (blue  arrowed  lines)  from  B770  to  2KS, 
4KS,  5KS,  6KS,  7KS,  and  15KS  sites.  The  elevation  profile  from  B770  to  the 
15KS  target  is  shown  below  the  map. 
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Figure  10.  Beam  height  and  elevation  profile  from  B 770  to  the  15KS  target. 
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propagation  paths  to  all  of  the  target  sites  were  within  5°  of  horizontal,  and  within  100  m  of  the 
ground  at  all  locations  along  the  paths.  The  elevation  profile  from  B770  to  the  15KS  target  is 
also  shown  at  the  bottom  of  Figure  9,  and  is  representative  of  the  general  elevation  profile 
characteristics  for  all  of  the  target  sites.  A  more  detailed  plot  showing  the  beam  height  and 
elevation  profile  from  B770  to  the  15KS  target  is  given  in  Figure  10.  A  weather  station  was 
operated  next  to  B770  to  provide  atmospheric  pressure,  temperature,  and  humidity 
measurements  during  the  experiments.  This  information  was  input  into  HITRAN-PC,  which 
was  then  used  to  calculate  the  expected  atmospheric  transmission  along  the  propagation  path. 

Photographs  of  the  target  sites  are  shown  in  Figures  11-15.  Figure  11  shows  blueboard  and 
flame-sprayed  aluminum  (FSA)  targets  at  the  2KS  site  used  for  pointing  alignment  and  system 
calibration.  Figure  12  shows  the  blueboard  and  berm  area  targets  at  the  4KS  site.  At  the  time 
the  photograph  in  Figure  12  was  taken,  only  a  single  blueboard  target  was  present.  In  early 
April  2003  another  blueboard  target  using  non-weathered  blueboard  panels  was  constructed 
and  placed  between  the  old  blueboard  target  and  berm  illumination  locations.  Aluminum 
alignment  targets  and  the  berm  illumination  area  for  the  7KS  site  are  shown  in  Figure  13.  The 
15KS  blueboard  target  is  shown  in  Figure  14.  Two  different  15KS  target  sites  were  used 
during  the  HD/DD  DC  experiments.  The  first  15KS  site  was  located  at  ~  14.9  km,  and  was  used 
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until  the  end  of  2002.  The  second  15KS  site  was  located  at  ~  15.3  km,  and  was  used  starting  in 
2003.  The  great  majority  of  the  experimental  measurements  were  conducted  using  the  second 
site,  and  all  15KS  information  given  in  this  report  pertains  to  this  location,  unless  specifically 
stated  otherwise.  A  photograph  of  the  propagation  path  from  B 770  to  the  15KS  target  board 
is  shown  in  Figure  15. 


Figure  11.  Photographs  of  left  and  right  sides  of  the  2KS  site,  showing  the  blueboard 


target  (left,  8’  x  8’)  and  flame-sprayed  aluminum  (right,  4’  x  4’)  calibration 
targets. 


Figure  12.  Photograph  of  4KS  target  site,  showing  blueboard  target  (right  of 
photograph)  and  approximate  berm  illumination  area  (yellow  circle,  not 
exactly  to  scale).  In  April  2003  a  second  target  using  non- weathered 
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blueboard  panels  was  constructed  and  placed  on  the  berm  (yellow 
rectangle). 


Figure  13.  Photograph  of  7KS  target  site,  showing  aluminum  (right)  and  plywood  (left) 
calibration/pointing  target  boards  and  approximate  berm  illumination  area 
(yellow  circle,  not  exactly  to  scale). 
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Figure  14.  Photograph  of  15KS  blueboard  target. 


Figure  15.  Photograph  of  propagation  path  from  B770  to  15KS  blueboard  target. 
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Panels  of  blueboard  insulation  (8’  x  4’  x  2”)  were  specially  prepared  for  use  at  the  target  sites. 
The  preparation  method  consisted  of  removing  the  external  extruded  surface  using  a  120-grit 
belt  sander,  thereby  exposing  the  diffuse  high  reflectivity  internal  material.  It  was  necessary 
to  remove  the  external  surface  on  all  surfaces  to  prevent  warping  of  the  blueboard  sheets. 
The  dimensions  of  the  blueboard  target  at  the  4KS  site  were  12’  long  by  8’  wide’,  and  the 
target  was  placed  at  approximately  a  60°  angle  of  incidence,  resulting  in  a  6’  long  x  8’  wide 
target  cross  section.  The  beam  diameter  at  the  4KS  target  is  approximately  4’  (285  prad  full- 
angle  divergence  in  output  space)  to  the  second  null  of  the  transmitted  beam  far-field 
distribution,  which  includes  ~  85%  of  the  total  beam  energy.  The  blueboard  target  at  the  15KS 
site  was  24’  length  x  16’  width  at  a  45°  angle  of  incidence,  resulting  in  a  17’  x  16’  target  cross 
section.  The  beam  diameter  at  the  15KS  target  is  approximately  9’  to  the  second  null  of  the 
far-field  distribution. 

Dirt  berms  were  also  used  as  targets  at  the  4KS  and  7KS  sites.  The  dirt  berm  reflectivity, 
calculated  from  calibrated  direct  detection  measurements  at  the  4KS  site,  was  determined  to 
approximately  match  a  Lambertian  reflector  with  a  total  reflectivity  of  3%,  for  angles  from 
normal  to  near  grazing  incidence.  For  the  approximate  incidence  angle  of  60°  for  the  berms,  the 
reflectivity  is  therefore  approximately  0.005  sr'1.  The  blueboard  and  berm  solid-angle 
reflectivities  are  given  in  Table  RR1  for  the  various  target  angles  used  in  the  HD/DD  DC 
experiments.  The  blueboard  provides  a  high  reflectivity  for  10  pm  radiation,  and  the 
reflectivity  does  not  decrease  as  quickly  with  increasing  incidence  angle  as  would  occur  for  a 
Lambertian  target.  It  should  also  be  noted  that  the  blueboard  reflectivities  were  measured  in 
the  laboratory  immediately  after  the  blueboards  were  prepared.  The  radiometric  results  which 
will  be  presented  later  in  this  report  suggest  that  the  blueboard  reflectivity  decreased 
significantly  as  the  boards  weathered,  since  they  were  continuously  exposed  to  the  outside 
elements  once  they  were  installed  at  the  target  sites.  Laboratory  reflectivity  measurements 
have  not  been  conducted  on  the  weathered  blueboards. 


Table  2.  Blueboard  and  dirt  berm  reflectivities 


Target 

Material 

Incidenc 

e 

Angle 

Reflectivity 

(sr1) 

Blueboard 

0° 

0.094 

Blueboard 

45° 

0.070 

Blueboard 

60° 

0.064 

Dirt  berm 

60° 

0.005 
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3.  Transmit  and  Receive  Temporal  and  Frequency  Characteristics 

A  detailed  understanding  of  the  temporal  and  frequency  characteristics  of  the  laser  transmit 
and  receive  signal  is  required  to  develop  accurate  analysis  procedures.  The  modelocked 
(spread  spectrum)  operation  implemented  to  provide  speckle  mitigation  in  the  heterodyne 
DIAL  measurements  differs  significantly  from  the  more  commonly  used  single-mode 
heterodyne  technique,  and  necessitates  more  demanding  system  requirements  (primarily  the 
ability  to  record  wideband  heterodyne  signals,  rather  than  narrowband  signals  centered  on 
the  heterodyne  intermediate  frequency).  The  signal-to-noise  ratio  (SNR)  advantages  of 
narrowband  operation  can  be  recovered  in  post-processing  of  the  wideband  signals,  because 
of  the  frequency-fence  structure  of  the  laser  modes.  As  will  be  shown,  the  laser  frequency 
modes  are  narrow,  and  occur  at  spacings  determined  by  the  laser  cavity  length,  allowing  the 
processing  bandwidth  to  be  reduced  to  only  those  regions  where  the  return  signal  occurs. 
Since  most  of  the  measurements  reported  in  this  paper  have  reasonably  good  heterodyne  SNR 
values,  this  processing  bandwidth  reduction  was  studied,  but  was  not  implemented  for  the 
analysis  results  presented. 

A.  Modelocked  Characteristics 

The  temporal  and  temporal  frequency  characteristics  for  a  single  outgoing  modelocked  laser 
pulse  are  shown  in  Figure  16.  The  outgoing  (monitor)  pulse  characteristics  were  captured 
using  direct  detection  of  a  small  portion  (<  1%)  of  the  outgoing  beam  energy.  Figure  16(a) 
shows  the  common  C02  laser  temporal  pulse  shape,  which  is  composed  of  a  gain-switch  spike 
and  a  relaxation  tail.  An  expanded  view  showing  the  modelocked  pulse  train  is  given  in 
Figure  16(b).  The  temporal  frequency  power  spectrum  showing  the  multiple  cavity 
longitudinal  modes  for  the  direct  detection  monitor  pulse  is  shown  in  Figure  16(c).  The 
modelocked  characteristics  are  achieved  by  placing  into  the  BBO  laser  cavity  an  acousto¬ 
optic  modulator  (AOM)  with  a  sinusoidal  variation  which  matches  the  cavity  round-trip  time. 
The  AOM  sinusoidal  variation  forces  the  laser  cavity  longitudinal  modes  to  be  in  phase  with 
each  other,  resulting  in  the  modelocked  temporal  pulse  train.  Modelocking  also  tends  to 
generate  more  longitudinal  cavity  modes  with  more  stable  amplitudes  than  allowing  the  laser 
to  free-run  with  multiple  longitudinal  modes  (as  will  be  seen  in  Figures  19  to  21).  The 
variations  from  perfect  theoretical  modelocking  operation  seen  in  the  BBO  laser  pulses  are 
believed  to  result  mainly  from  the  gradual  degradation  of  the  optical  coatings  on  the  intra¬ 
cavity  AOM.  The  AOM  was  coated  primarily  for  operation  at  11.15  jam  (13C1602  10P20),  but 
was  primarily  used  in  the  HD/DD  DC  experiments  over  the  wavelength  range  10.25  to  10.70 
jam.  In  addition,  the  AOM  is  an  element  in  an  open-air  laser  cavity,  and  is  therefore  in  a  high 
energy  density  region  which  is  not  completely  protected  from  dust  particles  and 
contaminants.  An  increase  in  the  visible  degradation  of  the  optical  coating  and  in  the 
difficulty  in  obtaining  optimal  modelocking  operation  was  noted  during  the  progress  of  the 
experiments,  starting  with  the  insertion  of  the  modelocker  in  January  2002,  and  continuing 
until  the  final  experiments  in  August  2003. 
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The  received  temporal  heterodyne  signal  for  modelocked  operation  is  shown  in  Figure  17(a), 
with  an  expanded  version  shown  in  Figure  17(b).  Since  the  transmitted  modelocked  laser 
beam  contains  multiple  longitudinal  modes,  the  transmit  laser  and  local  oscillator  were  not 
frequency  locked,  and  multiple  intermediate  frequencies  (up  to  100’s  of  MHz)  are  present  in 
the  heterodyne  signal.  The  heterodyne  modulation  of  the  modelocked  micropulses  is 
especially  evident  in  Figure  17(b). 


Monitor  Temporal  Pulse  Monitor  Temporal  Pulse 
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Frequency  (MHz) 


Figure  16.  Modelocked  laser  temporal  pulse  (direct  detection),  (a)  entire  pulse 
waveform,  (b)  expanded  time  scale,  (c)  temporal  frequency  spectrum. 
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Figure  17.  Heterodyne  receive  signal  from  a  modelocked  laser  pulse,  (a)  entire  pulse 
waveform,  (b)  expanded  time  scale. 


Heterodyne  temporal  frequency  spectra  for  modelocked  operation  are  shown  in  Figure  18. 
The  superimposed  heterodyne  spectra  of  the  signal  and  local  oscillator  [calculated  from  70-83 
(is  in  Figure  17(a)]  are  shown  in  Figure  18(a).  The  homodyne  spectra  for  the  local  oscillator 
[calculated  from  61-66  qs  in  Figure  17(a)],  with  the  DC  component  removed,  are  shown  in 
Figure  18(b).  The  local  oscillator  output  actually  consists  of  a  dominant  single  longitudinal 
cavity  mode,  with  two  smaller  features,  which  are  denoted  as  ‘spurs’.  The  spurs  occur  at  the 
single  and  double  AOM  frequency  shifts  (approximately  77  and  154  MHz).  They  are  believed 
to  result  from  slight  impurities  in  the  local  oscillator  dual  intra-cavity  AOMs,  which  causes 
some  of  the  unshifted  (in  wavelength)  beam  to  be  scattered  into  the  direction  of  the  shifted 
beam.  Since  this  can  happen  in  both  of  the  AOMs,  the  local  oscillator  output  can  consist  of 
three  frequencies,  the  desired  frequency,  a  single- shifted  frequency,  and  a  double-shifted 
frequency.  The  exact  frequency  shifts  are  slightly  different  for  each  local  oscillator 
wavelength,  since  the  AOM  frequency  shift  required  is  dependent  upon  the  output 
wavelength  desired. 
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Figure  18.  Heterodyne  receive  signal  temporal  spectra  from  a  modelocked  laser  pulse. 

(a)  signal  +  local  oscillator,  (b)  local  oscillator  only,  (c)  normalized  signal 
[(signal  +  local  oscillator)  /  (local  oscillator)]. 

The  optical  power  in  the  spurs  is  more  than  3  orders  of  magnitude  less  than  that  in  the  main 
local  oscillator  output,  and  would  not  be  evident  in  most  heterodyne  measurements.  In  many 
laboratory  or  short-range  heterodyne  measurements,  the  return  signal  level  would  be  much 
higher  than  that  in  the  spurs.  Also,  for  many  heterodyne  systems  a  single  longitudinal  mode 
transmitter  is  used,  and  the  intermediate  frequency  and  associated  receiver  passband  would 
be  significantly  lower  than  the  first  spur  frequency.  In  the  case  of  the  HD/DD  DC 
measurements,  however,  wideband  heterodyne  detection  is  required  because  of  the  multiple 
longitudinal  mode  (MLM)  nature  of  the  transmitter,  and  the  received  signal  optical  power 
levels  are  of  the  same  order  of  magnitude  as  those  of  the  local  oscillator  spurs.  A  number  of 
physical  methods  to  eliminate  or  significantly  reduce  the  spurs  were  studied,  but  it  was 
determined  that  a  simpler  and  less  costly  option  was  to  remove  the  spurs  in  the  processing 
algorithms.  The  normalization  procedure  used  for  spur  removal  is  described  in  more  detail  in 
the  data  processing  section.  It  should  also  be  noted  that  the  noise  level  evident  in 
Figure  17(a)  is  mainly  caused  by  the  154  MHz  beating  of  the  main  LO  frequency  component 
with  the  second  spur,  and  not  by  the  local  oscillator  shot  noise. 
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The  heterodyne  signal  temporal  frequency  power  spectrum  for  modelocked  operation  is 
shown  in  Figure  18(c),  after  the  normalization  procedure  has  been  applied.  As  seen  in  the 
Figure,  the  local  oscillator  spurs  at  approximately  77  and  154  MHz  are  removed  by  the 
normalization  procedure,  and  only  the  transmit  laser  cavity  modes  are  present.  The  laser 
cavity  modes  are  separated  by  approximately  44  MHz,  which  corresponds  to  the  cavity  round- 
trip  time  (~  23  ns  for  a  3.4  m  laser  cavity  length).  The  double  peak  nature  of  the  frequency 
modes  in  Figure  18(c)  is  believed  to  result  from  the  DC  wrap-around  inherent  in  the 
periodogram  computation  process,  with  the  negative  frequency  modes  happening  to  nearly 
fall  back  onto  the  positive  frequency  modes.  Multiple-peak  modes  can  also  result  from  the 
transmit  laser  containing  some  multiple  TEM  output,  but  this  is  not  believed  to  be  the  case  in 
this  instance. 

B.  Non-Modelocked  Characteristics 

The  temporal  and  temporal  frequency  characteristics  for  a  non-modelocked  output  pulse  are 
shown  in  Figure  19,  similar  to  the  presentation  in  Figure  16  for  a  modelocked  pulse.  The  gain- 
switch  spike  and  relaxation  tail  are  again  evident.  A  small  after  pulse  (at  39-40  qs  in  the  plot), 
probably  resulting  from  continued  RF  pumping  during  the  pulse  and  incomplete  cavity 
dumping,  is  also  present  for  this  specific  pulse.  The  expanded  scale  in  Figure  19(b)  shows  a 
significant  difference  between  non-modelocked  and  modelocked  operation.  The  effects  of 
MLM  operation  are  evident  in  the  partial  micropulse  temporal  structure,  but  the  pulses  are 
overlapped  rather  than  distinct  because  the  longitudinal  cavity  modes  are  not  phase  matched, 
as  is  the  case  in  modelocked  operation.  The  frequency  modes  are  shown  in  Figure  19(c). 
Multiple  modes  are  present,  but  as  seen  in  comparison  with  Figure  16(c),  there  are  fewer 
modes  for  non-modelocked  operation,  with  a  more  rapid  decrease  in  amplitude  away  from  the 
main  mode.  This  mode  structure  indicates  that  non-modelocked  operation  will  result  in  some 
speckle  mitigation,  but  not  as  much  as  for  modelocked  operation. 

The  heterodyne  receive  signal  from  the  non-modelocked  pulse  is  shown  in  Figure  20.  The 
heterodyne  modulation  is  evident  in  Figure  20(b),  and  follows  the  pulse  envelope  seen  in 
Figure  19(b). 
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Figure  19.  Non-modelocked  laser  temporal  pulse  (direct  detection),  (a)  entire  pulse 
waveform,  (b)  expanded  time  scale,  (c)  temporal  frequency  spectrum. 
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Figure  20.  Heterodyne  receive  signal  from  a  non-modelocked  laser  pulse,  (a)  entire 
pulse  waveform,  (b)  expanded  time  scale. 
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The  heterodyne  signal  spectra  for  non-modelocked  operation  are  shown  in  Figure  21,  in 
similar  fashion  to  that  for  modelocked  operation  shown  in  Figure  18.  In  this  case,  the 
normalized  spectrum  shown  in  Figure  21(c)  shows  the  effect  of  the  DC  wrap-around,  where  the 
set  of  frequency  modes  marked  as  (1)  can  be  considered  as  the  positive  frequency  modes,  and 
those  marked  as  (2)  can  be  considered  as  the  wrapped  negative  frequency  modes.  With  this 
unwrapping,  it  can  again  be  seen  that  the  frequency  mode  spacing  is  approximately  44  MHz, 
as  expected.  Comparing  Figure  21(c)  and  18(c),  it  can  be  seen  that  the  modelocked  spectrum 
contains  more  modes  with  higher  amplitudes  than  the  non-modelocked  spectrum,  indicating 
that  modelocked  operation  will  result  in  a  greater  degree  of  speckle  mitigation. 
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Figure  21.  Heterodyne  receive  signal  temporal  spectra  from  a  non-modelocked  laser 
pulse,  (a)  signal  +  local  oscillator,  (b)  local  oscillator  only,  (c)  normalized 
signal  [(signal  +  local  oscillator)  /  (local  oscillator)]. 
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4.  Data  Analysis  and  Results 

A.  Overview 

As  described  previously,  the  development  of  the  Coherent  Remote  Optical  Sensing  System 
(CROSS)  began  in  October  1998,  as  the  Laser  Airborne  Remote  Sensing  (LARS)  program  was 
nearing  completion.  Equipment  development,  laboratory  system  characterization  and  DIAL 
measurements,  and  preliminary  open-air  measurements  were  conducted  during  1999,  2000,  and 
2001.  Although  the  system  continued  to  be  upgraded  and  improved,  the  system  components 
remained  essentially  the  same  after  the  HD/DD  DC  measurements  began  in  February  2002. 
The  data  sets  included  in  the  radiometric,  DIAL,  and  system  characterization  analyses  for  this 
report  are  shown  in  reverse  chronological  order  in  Appendix  A,  Table  Al. 


B.  Radiometric  Theoretical  Analysis 

The  differential  absorption  lidar  (DIAL)  technique  utilizes  the  variations  in  absorption  at 
different  interrogation  wavelengths  to  determine  a  ‘spectral  fingerprint’  in  order  to  detect  and 
quantify  chemicals  in  the  laser  path.  In  the  simplest  scenario,  this  can  be  thought  of  as  a  two- 
wavelength  system,  with  one  wavelength  tuned  to  an  absorbing  feature  of  a  chemical,  and  the 
other  wavelength  used  as  a  reference  which  is  not  absorbed.  In  a  more  general  and  realistic 
scenario,  the  DIAL  problem  can  be  posed  in  terms  of  using  a  sequence  of  N  laser  wavelengths 
to  determine  the  concentrations  for  M  chemicals.  For  the  general  case,  the  commonly 
presented  2-wavelength  DIAL  approach  may  not  be  usable,  since  the  absorption  of  one 
chemical  may  interfere  with  either  the  on-  or  off-absorption  wavelength  of  another  chemical.  It 
should  be  noted  that  conceptually  there  is  no  difference  in  the  DIAL  technique  between 
measuring  the  concentrations  (see  Eq.  1)  of  the  IVf  atmospheric  and  Mn  target  gases, 
although  more  laser  wavelengths  are  required  to  measure  more  gases.  In  the  HD/DD  DC 
experiments,  the  Ma  contributions  were  determined  by  measuring  the  atmospheric  pressure, 
temperature,  and  humidity,  and  determining  the  expected  atmospheric  transmission  (Ta)  using 
an  atmospheric  model  (HITRAN-PC).  The  DIAL  equation  for  using  N  laser  wavelengths  to 
measure  M  (specifically  Mm)  target  gases  can  be  written  in  the  following  form 
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The  above  equation  assumes  a  hard  target  is  providing  the  return  signal,  and  therefore  a 
column-content  measurement  is  being  made.  If  a  distributed  target,  such  as  the  atmosphere  or 
aerosols,  provides  the  return  mechanism,  the  received  signal  will  be  spread  in  time,  and  a 
range-resolved  measurement  will  be  made.  The  chemometric  analysis  is  performed  similarly 
for  both  column-content  and  range-resolved  measurements,  but  must  be  done  for  each  range 
window  in  the  range-resolved  case.  Most  long  range  systems  utilize  column-content 
measurements  because  of  the  significantly  larger  return  signal  level  occurring  from  hard 
targets  as  compared  to  atmospheric  or  aerosol  backscatter. 


C.  Direct  Versus  Heterodyne  Detection  Signal  Processing  Issues 

The  differences  in  the  physics  of  the  detection  processes  between  direct  and  heterodyne 
detection  requires  a  complete  and  accurate  understanding  of  the  signal  processing  techniques 
that  are  used.  One  issue  that  needs  to  be  considered  in  the  heterodyne  detection  case,  is  that 
the  commonly  derived  heterodyne  SNR  relation 
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for  continuous  wave  heterodyne  operation  is  not  directly  applicable  to  the  case  of  a  pulsed 
laser,  where  the  signal  will  be  integrated  to  provide  a  result  proportional  to  the  received 
energy  Er¬ 
in  the  case  of  direct  detection,  the  relationships  between  the  current  out  of  the  detector  IDD, 
the  instantaneous  power  received  Pr,  and  the  received  energy  Er  are  given  by 

XDD  (' *)  «  |er  ( ' Of  =  |a  r(t)cos  ((0rt  +  cpr  )|2  =  ^|ar|2  =  Pr  (t ) 

2  •  (3) 

UDd=  I  Idd  (t)dtoc  J  Pr(t)dt=Er 

AT  AT 

The  measured  value  UDD  is  the  parameter  required  for  both  radiometric  and  DIAL  analysis. 
Since  it  is  proportional  to  E).,  it  is  then  also  directly  proportional  to  the  total  transmission  of 
the  target  gases  of  interest.  The  processing  is  actually  more  complicated  than  the  simple 
relation  given  in  Eq.  3,  but  the  key  relationship  between  UDD  and  Er  is  retained.  The  details  of 
the  processing  are  given  in  Appendix  C,  where  the  code  for  the  primary  direct  detection 
Matlab  calculation  program  DATREDUC.m  is  reproduced. 

Determining  an  equivalent  heterodyne  parameter  UHd  to  the  direct  detection  parameter  UDD  is 
significantly  more  complicated.  There  are  many  different  ways  to  calculate  a  UHD  parameter 
for  heterodyne  detection,  all  of  which  are  valid  (but  not  necessarily  optimal)  as  long  as  UHd  is 
directly  proportional  to  the  received  energy  In  the  case  of  heterodyne  detection,  the 
current  out  of  the  detector  is  given  by 

!hD  (0  06  K  (  0  +  £LO  ( t)|  =  |ar  (l)cos  (C0rt  +  cpr  )  +  aLO  (cos  ^LO1  +  9lO  )| 

=  |ar(t)|2  cos2  (®rt  +  <pr )+  2|ar  (t)aLO|cos(cort  +  (pr)cos((0LOt+  (pLO ) 

+ 1  aLO  P  cos2  (  “lo1  +  9lo  ) 

12  12  •  (4) 

=  J  lar  Wl  +  2  la  LO I  +2lar(t)aLo|cos[(ttV  _®Lo)t+((Pr_(PLo)] 

=  Pr(t)  +  PLO  +  2-JP,-  (t )  Plo  COS  [  (tO,-  -C0LO  )  t+  (cpr  -  CpLO  )] 

=  2-y/P,  (t )  Plo  COs[(cOr  -®LO)t  +  (cpr-(pLo)] 

The  above  derivation  uses  the  facts  that  the  heterodyne  detector  does  not  respond  at  optical 
frequencies,  that  the  detector  or  electronics  will  not  pass  baseband  (DC)  signals  (thereby 
removing  the  PL0  contribution),  and  that  the  local  oscillator  power  PL0  is  significantly  larger 
than  Pr  (thereby  removing  the  Pr  term).  The  above  relation  has  been  significantly  simplified, 
and  a  more  complete  relationship  would  need  to  include  the  contribution  of  the  LO  spurs,  as 
seen  previously  in  Figures  18(a)  and  21(a).  Since  the  optical  power  in  the  LO  spurs  is  of  the 
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same  order  of  magnitude  as  the  received  instantaneous  power  IJ,  the  beat  contribution 
between  the  spurs  and  the  main  LO  mode  is  also  present  in  the  detector  output,  and  occurs  at 
harmonics  of  the  AOM  shift  frequency  (approximately  77  and  154  MHz).  The  presence  of  the 
LO  spurs  led  to  the  construction  of  a  customized  processing  technique  designed  to  remove 
their  effect  in  the  calculation  of  UHd- 


The  calculation  of  UHD  is  primarily  based  upon  the  use  of  the  power  spectral  density  of  tD- 
The  power  spectral  density  Ss+N  of  IHd  is  computed  over  a  temporal  window  where  the  return 
and  local  oscillator  signals  are  present,  and  a  separate  power  spectral  density  SN  is  computed 
over  a  temporal  window  where  the  local  oscillator  but  no  signal  is  present.  UhD  is  then 
calculated  from  the  integration  over  temporal  frequency  space  of  the  point-by-point 
normalization  of  Ss+N  by  SN 
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Utilizing  Rayleigh’ s  energy  theorem,  it  can  be  shown  that  UHd  is  thereby  directly  proportional 
to  E,.,  as  desired.  The  factor  Nenf  is  termed  the  excess  noise  factor,  and  results  from  the  fact 
that  point-by-point  normalization  is  a  biased  processing  procedure,  due  to  the  nonlinearity  of 
the  division  operation.  It  can  be  shown  that  in  the  case  of  no  signal,  for  the  parameters  used 
in  the  heterodyne  processing  shown  in  this  report,  the  division  of  two  random  noise  signals 
results  in  a  value  of  1+Nenf  of  approximately  1.150,  as  compared  to  a  value  of  unity  for  an 
unbiased  processing  procedure.  Nenf  is  strictly  a  function  of  the  heterodyne  SNR,  with  Nenf 
approaching  zero  at  higher  SNRs.  Since  the  Nenf  correction  is  more  important  at  low  SNRs 
than  at  higher  SNRs,  a  constant  value  of  Nenf  =  0. 150  corresponding  to  the  no  signal  case  was 
used  in  the  processing.  The  details  of  the  heterodyne  processing  are  also  given  in  Appendix 
C,  where  the  code  for  the  primary  heterodyne  detection  Matlab  calculation  program 
DATREDUC_HD_X3.m  is  reproduced. 


D.  Chemometric  Theoretical  Analysis 

The  chemometric  processing  is  performed  in  exactly  the  same  manner  for  both  direct  and 
heterodyne  detection,  since  the  computed  parameters  UDD  and  UHd  are  both  directly 
proportional  to  E,.,  and  therefore  to  the  total  transmission  of  the  target  gases  of  interest.  The 
most  accurate  and  least  complicated  method  to  isolate  the  absorption  caused  by  M  target 
gases  is,  if  possible,  to  perform  the  same  measurement  in  the  absence  of  the  target  gases  as 
will  be  performed  with  the  gases  present.  The  initial  step  in  the  analysis  procedure  is  to 
normalize  the  receive  energy  by  the  transmit  energy.  The  data  are  then  background 
normalized  by  a  no-gas  spectrum,  leaving  the  spectral  transmission  from  the  target  gases  only. 
For  a  sequence  of  N  laser  lines  and  M  absorbing  gases,  the  set  of  equations  used  in  the 
chemometric  analysis  can  be  written  as 
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T(A-i)  =exp  (-2k(  ^  ,gi )  CL(g! ))  x  exp  (-2k  (A1;  g2  )CL  (g2 ) )  xL 
xexp(-2K(A,1,gM)CL(gM)) 

T(A^)  =exp(-2K(A,2,g1)CL(g1))xexp(-2K(A,2,g2)CL(g2))xL 

xexp(-2K(A,2,gM)CL(gM ))  (6) 

M 

t(^n  )  =  exP  (— 2k(^n  ’Si  )cl(Si))  x  exP  (— 2k  (A-n  >  g2  )CL(§2)) x 
xexp  (-2k  (  An  ,  gjyj  )CL  (gM )) 

where  T(Aq)  is  the  total  transmission  at  wavelength  A*  including  the  contributions  from  all  of 
the  target  gases. 

A  matrix  equation  which  is  linear  in  the  unknowns  CL  (gj )  is  formed  by  taking  the  logarithm 
of  each  equation 
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Applying  singular-value  decomposition  (SVD)  matrix  analysis  or  other  techniques  to  the 
above  equation  provides  a  solution  c  for  the  CL  (gj  .  Some  of  the  advantages  of  the  above 

chemometric  analysis  are  that  all  chemical  concentrations  are  solved  for  simultaneously,  and 
that  all  of  the  spectral  information  available  is  used  in  the  analysis.  These  characteristics 
provide  definite  improvements  over  the  conventional  2-line  DIAL  analysis  for  handling 
chemicals  with  overlapping  spectral  features.  Some  inaccuracies  have  been  noted  for 
chemicals  with  overlapping  spectral  features  when  using  the  basic  SVD  analysis,  but  it  is 
believed  that  applying  a  scheme  to  de-weight  these  wavelengths  is  feasible,  and  would 
improve  the  accuracy  of  the  chemometric  results. 

E.  Radiometric  Results 

The  direct  detection  radiometric  results,  presented  as  the  ratio  of  the  expected  return  signal  to 
the  actual  return  signal  (FdD  =  I?bD/UDD)  are  shown  in  Figure  22.  The  results  are  from 
experiments  conducted  on  different  days  using  both  modelocked  and  non-modelocked 
waveforms,  and  are  separated  into  groups  for  the  four  different  targets  used.  Target  #1  was 
the  earth  berm  located  at  the  4KS  site,  Target  #2  was  a  highly  weathered  (2-year  exposure) 
blueboard  target  at  the  4KS  site,  Target  #3  was  a  moderately  weathered  (2-month  exposure) 
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blueboard  target  at  the  4KS  site,  and  Target  #4  was  the  7KS  earth  berm.  Data  points  that  were 
obvious  outliers  have  been  removed  from  the  plot.  As  seen  in  the  Figure,  there  is  very  good 
agreement  between  the  measured  and  predicted  values  (within  a  factor  of  1  -  2.5)  for  the  4  km 
and  7  km  topographic  targets.  This  agreement  level  is  reasonable,  considering  that  the  DIAL 
measurements  had  the  most  emphasis,  and  the  radiometric  characterization  was  therefore  not 
as  extensive  as  had  been  conducted  in  past  experiments.  Previous  precise  radiometric 
measurements  on  the  direct  detection  system  conducted  in  1998  illustrated  that  the  level  of 
radiometric  agreement  achievable  with  the  system  fully  optimized  was  within  a  factor  of  1  to 
1.4.4  The  discrepancies  seen  in  Figure  22  for  the  blueboard  targets  are  most  likely  caused  by 
inaccuracies  in  the  assumed  blueboard  target  reflectivity.  The  simulation  results  used 
reflectivity  values  from  measurements  of  the  blueboards  immediately  after  they  were 
fabricated,  and  before  being  continuously  exposed  to  the  outside  elements.  The  blueboard 
reflectivity  is  expected  to  decrease  as  the  target  material  weathers  and  blueboard  dust 
particles  build  up  on  the  surface,  and  is  a  reasonable  explanation  for  the  results  from  the 
blueboard  targets.  The  previous  measurements  in  1998  were  conducted  immediately  after  the 
blueboard  targets  had  been  fabricated  and  characterized,  and  both  blueboard  and  berm  return 
signals  agreed  with  predictions  within  the  factor  of  1  -  1.4  mentioned  previously. 


Target 

Figure  22.  Direct  detection  radiometric  results. 
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The  heterodyne  radiometric  results  are  shown  in  Figure  23.  Similarly  to  the  direct  detection 
results,  the  ratio  of  expected  to  measured  return  signal  FHD  =  i3^hd  /  Urd  is  computed.  This 
ratio  is  further  normalized  by  the  direct  detection  ratio  FbD,  in  order  to  remove  any  errors 
caused  by  inaccuracies  in  the  assumed  reflectivities.  Data  points  that  were  obvious  outliers 
have  been  removed  from  the  plot.  As  seen  in  the  Figure,  the  heterodyne  return  signal  is 
within  a  factor  of  approximately  6  to  20  of  the  expected  value.  The  heterodyne  return  signal 
simulation  includes  heterodyne  efficiency  and  signal  reduction  factors  resulting  from  beam 
pattern  mismatch  between  the  LO  (uniform)  and  return  signal  (Airy)5,6;  distortion  of  the  return 
signal  phase  front  by  atmospheric  turbulence  (assuming  a  relatively  high  turbulence  level  of 
Cn  =  10  m  )  ;  and  loss  of  signal  caused  by  transmitter  beam  (pointing)  jitter.  Factors 
which  were  not  able  to  be  measured  and  are  not  included  in  the  simulation  include  non¬ 
uniformity  in  the  detector  quantum  efficiency  (probably  a  negligible  factor);  misalignment  of 
the  transmitter  footprint  and  receiver  field-of-view  (also  probably  a  negligible  factor);  phase 
front  curvature  from  misfocusing  of  the  heterodyne  receiver  telescope  (also  probably 
negligible);  and  angular  misalignment  of  the  LO  and  receiver  signal  on  the  detector  (possibly  a 
factor  of  ~2  reduction).  It  should  also  be  noted  that  the  results  of  Frehlich  and  Kavaya  [1991] 8 
suggest  that  the  coherence  loss  factor  caused  by  atmospheric  turbulence  may  be  significantly 
worse  than  that  given  by  Clifford  and  Wandzura  [198 1] 7  and  used  in  this  paper. 
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Figure  23.  Normalized  heterodyne  radiometric  results. 
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Overall,  the  difference  between  the  expected  and  measured  heterodyne  signal  levels  is  larger 
than  expected,  but  not  significantly,  and  could  most  likely  be  significantly  reduced  with  a 
more  precise  system  optimization.  The  best  agreement  noted  to  date  between  expected  and 
actual  return  heterodyne  signal  levels  for  comparable  systems  is  a  factor  of  ~49,  with  most 
systems  reporting  a  5  -  10  dB  (factor  of  3  -  10) 10,11,12,13  reduction.  With  the  possible  loss 
factors  not  included  in  the  analysis  of  the  CROSS  system,  and  the  added  complexity  inherent 
in  a  fully  wavelength-agile  heterodyne  system,  a  reduction  factor  of  6  -  20  compares  favorably 
with  the  previous  reduction  factors  from  other  researchers  of  3-  10. 

F.  Chemometric  Results  (Case  Studies) 

The  absorption  spectra  for  SF6  and  NH3  are  shown  in  Figure  24.  The  x’s  indicate  the 
absorption  coefficients  for  the  chemicals  at  the  specific  wavelengths  in  the  10R  and  10P  bands 
for  the  12C1602  laser  isotope  being  used  for  the  experiments.  SF6  has  a  broad  absorption 
feature,  which  will  affect  most  of  the  10P  laser  lines.  NH3  has  a  much  narrower  absorption 
spectrum,  and  has  its  only  significant  absorption  on  the  10R14  laser  line.  For  the  HD/DD  DC 
experiments,  only  3-4  of  the  10R  lines  shown  on  the  plot  (between  ~  969-982  cm1)  were  used  in 
the  laser  sequence. 
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Figure  24.  Absorption  spectra  overlap  with  12C1602  laser  wavelengths  for  SF6  and  NH3. 


The  results  from  the  HD/DD  DC  tests  using  the  spread  spectrum  (modelocked)  transmit 
waveform  and  the  diffuse  blueboard  target  at  the  4  km  site  are  shown  in  Figure  25.  The 
absorption  plots  for  both  the  direct  detection  and  heterodyne  return  signals  are  shown  on  the 
left  of  the  Figure.  The  horizontal  axis  in  these  plots  corresponds  to  time,  while  the  vertical  axis 
shows  the  laser  wavelengths  transmitted.  From  sequences  1  to  40,  there  was  no  SF6  present 
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in  the  absorption  cell.  Around  sequence  40,  a  first  SF6  insertion  was  made,  and  a  second  SF6 
insertion  was  made  near  sequence  70.  The  SF6  absorption  signature  is  seen  clearly  in  both  of 
the  absorption  plots,  although  the  direct  detection  plot  is  cleaner.  Rough  calculations  of  the 
expected  single  shot  SNR  give  a  direct  detection  SNRDD  ~  5,  resulting  from  speckle  spatial 
averaging,  and  a  heterodyne  single  shot  SNRHd  ~  2.5,  resulting  from  approximately  6  effective 
modes  in  the  transmitted  modelocked  spectrum.  After  the  DIAL  processing  techniques  are 
applied,  the  SNRs  shown  in  the  absorption  plots  become  SNRDD  ~  14  and  SNRHd  ~  7.  The 
graph  on  the  bottom  right  of  the  Figure  shows  the  computed  concentration-length  (CL) 
products  as  a  function  of  time  for  both  techniques.  The  table  on  the  top  right  summarizes  the 
measured  CLs  in  each  of  the  three  reference  /  insertion  regions.  The  results  from  the  two 
detection  techniques  differ  by  less  than  12%,  which  is  reasonably  good  agreement  given  the 
number  of  uncontrolled  variables  in  the  experiment. 

The  results  from  the  HD/DD  DC  tests  using  the  spread  spectrum  (modelocked)  transmit 
waveform  and  the  topographic  target  at  the  7.5  km  site  are  shown  in  Figure  26.  The  SF6 
insertion  timeline  is  essentially  the  same  as  that  described  for  Figure  25.  The  SF6  signature  is 
again  evident  in  both  of  the  absorption  plots.  The  expected  single  shot  SNRs  are  the  same  as 
for  the  4  km  measurements  (SNRDD  ~  5  and  SNRHd  ~  2.5).  After  the  processing  techniques  are 
applied,  the  SNRs  become  SNRDD  ~  14  and  SNRHd  ~  5,  indicating  that  the  direct  detection 
results  are  similar  at  both  ranges,  and  the  heterodyne  results  were  slightly  degraded  at  the 
longer  range.  In  this  case,  the  results  from  the  two  detection  techniques  differ  by  less  than 
20%,  which  is  still  reasonably  good  agreement. 


HD/DD  DIAL  Comparison  Test 
020926_000T1D  (Modelocking  On,  4KS  BB) 
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6.0 
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Figure  25.  Results  from  the  020926  experiment  using  the  blueboard  at  the  4  km  site. 
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Figure  26.  Results  from  020926  experiment  using  topographic  returns  at  the  7.5  km  site. 


Results  from  a  dual  insertion  of  SF6  and  NH3  are  shown  in  Figure  27.  The  SF6  was  inserted 
around  sequence  40,  and  the  NH3  was  inserted  near  sequence  70.  As  shown  in  Figure  24,  SF6 
exhibits  absorption  on  most  of  the  10P  laser  lines,  while  NH3  has  significant  absorption  only  at 
the  10R14  wavelength.  The  CL  values  for  SF6  and  NH3  measured  with  both  the  direct 
detection  and  heterodyne  receivers  are  shown  on  the  plot. 
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HD/DD  DIAL  Comparison  Test  03041 6_T2D 
(SF6  +  NH3,  Modelocking  On,  7KS  Berm) 


Figure  27.  Results  from  030416  dual  SF6  and  NH3  experiment  using  topographic  returns 
at  7.5  km. 


Heterodyne  and  direct  detection  DIAL  results  using  the  blueboard  target  at  a  range  of  15  km 
are  shown  in  Figures  28  and  29.  The  direct  detection  DIAL  results  are  significantly  noisier 
than  at  shorter  ranges,  which  is  expected  since  at  longer  horizontal  ranges  the  direct  detection 
system  is  signal-limited,  and  is  therefore  no  longer  in  the  speckle-limited  regime.  The  direct 
detection  signal  level  for  the  030612  data  set  is  extremely  low,  which  causes  significant  errors 
in  the  DIAL  intermediate  normalization  processing  algorithm,  and  appears  as  the  physically 
incorrect  result  of  the  measured  absorption  being  saturated.  The  heterodyne  DIAL  results 
also  show  some  degradation,  but  not  to  as  great  an  extent  as  the  direct  detection  results, 
which  indicates  that  the  heterodyne  system  is  still  operating  in  the  speckle-limited  regime,  as 
expected. 
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HD/DD  DIAL  Comparison  Test 
021 21 8T3D  (Modelocking  Off,  15KS  BB) 


Figure  28.  Results  from  the  021218  experiment  using  the  blueboard  at  the  15  km  site. 


i  i  HD/DD  DIAL  Comparison  Test  03061 2_T3D 
(NH3  +  SF6,  Modelocker  On,  15KS  Blueboard) 


“O 


1 


Sequences 

Measured  CL 

(PPm-m) 

Measured 

Ap(psi) 

Expected 
Ap  (psi) 

Reference 

5-35 
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#2  (NH3) 
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277  ±  305  (HDx3) 

36.6  (DD)* 

6.1  (HDx3) 
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*  Low  signs 
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Figure  29.  Results  from  the  030612  dual  SF6  and  NH3  experiment  using  the  blueboard  at 
15  km. 
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G.  Chemometric  Results  (Compilation  and  Comparisons) 

This  section  will  provide  a  summary  of  the  DIAL  chemometric  results  for  the  HD/DD  DC 
experiments.  Detailed  results  for  all  of  the  primary  HD/DD  DC  experiments  are  shown  in 
Appendix  B.  The  expected  and  measured  concentration-length  (CL)  products  from  all  of  the 
HD/DD  DC  data  sets  used  in  this  report  are  shown  in  Figure  30.  Prior  to  March  2004  (~  Day 
#60),  the  common  experiment  technique  was  to  begin  data  collection  with  no  gas  in  the  cell, 
then  add  SF6  approximately  1/3  of  the  way  through  the  data  set,  and  then  add  more  SF6 
approximately  2/3  of  the  way  through  the  data  set.  The  initial  SF6  insertion  CL  measurements 
are  shown  in  blue  in  the  Figure,  with  the  second  SF6  CL  measurements  shown  in  green.  After 
March  2004,  the  experiment  technique  was  changed  to  normally  have  the  first  insertion  with 
SF6  (shown  in  blue),  and  the  second  insertion  with  NH3  (shown  in  red).  The  dashed  lines  on 
the  plots  have  no  significance,  but  are  shown  to  more  clearly  indicate  the  different  gas 
insertions.  To  reduce  clutter  from  overlap  of  the  symbols  on  the  plot,  measurements  from  the 
4KS  blueboards  (4KS  BB)  were  shifted  left  by  Vi  day,  and  measurements  from  the  7KS  berm 
were  shifted  right  by  Vi  day. 


Day#  (2003) 

Figure  30.  Expected  and  measured  CL  values  from  HD/DD  DC  experiments. 
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Although  there  is  a  significant  amount  of  scatter  in  the  results  shown  in  Figure  30,  a  few 
general  conclusions  can  be  drawn.  During  the  time  period  of  these  measurements,  both  the 
heterodyne  and  direct  detection  DIAL  systems  were  being  continually  improved,  both  with 
respect  to  the  equipment,  and  to  the  alignment  and  experimental  procedures.  This  is 
evidenced  by  noting  that  prior  to  March  2004,  the  measured  CL  results  were  significantly 
lower  than  the  expected  values,  and  the  measured  CLs  never  significantly  exceeded  the 
expected  values.  After  March  2004,  the  measured  CLs  still  tended  to  be  lower  than  the 
expected  values,  but  were  closer  to  the  expected  value,  and  sometimes  exceeded  it.  This 
improvement  continued  until  the  final  measurements  made  in  August  of  2004. 

It  is  believed  that  the  measured  CL’s  being  consistently  lower  than  expected  probably  results 
from  a  systematic  error,  possibly  from  inaccuracies  in  the  gas  insertion  procedure  or  from 
errors  in  the  compensation  of  the  spectral  absorption  coefficients  for  the  experiment 
atmospheric  pressure  and  temperature.  Since  the  primary  purpose  of  the  HD/DD  DC 
experiment  series  was  to  compare  heterodyne  and  direct  detection  DIAL,  relative  results  were 
emphasized,  and  extraordinary  efforts  were  not  made  to  calibrate  the  system  for  absolute 
accuracy.  Previous  measurements  of  SF6,  NH3,  and  Freon  R-134A  made  during  ground  tests 
in  19984  with  only  the  LARS  direct  detection  DIAL  system  showed  slightly  more  accurate 
DIAL  results  than  achieved  with  the  HD/DD  DC  direct  detection  measurements.  The  more 
accurate  results  in  1998  were  possible  primarily  because  only  a  single  system  and  single  target 
were  being  used,  allowing  more  meticulous  alignment  and  system  calibration.  The  more 
complex  nature  of  the  HD/DD  DC  experiment  reduced  the  capability  to  optimize  all  of  the 
system  components  simultaneously,  and  would  be  expected  to  result  in  reduced  DIAL 
measurement  accuracy.  Even  with  the  slightly  reduced  accuracy,  the  comparisons  between 
the  heterodyne  and  direct  detection  DIAL  results  are  still  valid  and  provide  the  key 
information  required  to  analyze  and  compare  the  capabilities  of  the  techniques. 

To  better  study  the  effects  of  the  various  system  parameters  on  the  DIAL  accuracy,  the 
measurements  were  separated  into  groups  by  the  experiment  and  system  parameters.  A 
summary  of  the  DIAL  results  for  these  related  measurements  is  given  in  Table  3.  The  Table 
separates  the  measurements  by  the  detection  technique  used,  the  target  gas,  the  modelocking 
configuration,  the  time  period  when  the  measurements  were  made,  the  number  of  individual 
measurements  included  in  the  group,  and  the  average  and  standard  deviation  of  the  measured 
to  expected  concentration-length  (CL)  ratio.  It  should  be  noted  that  even  though  there  were 
more  than  158  individual  DIAL  measurements  made,  the  compiled  grouping  results  may  still 
not  be  statistically  significant,  since  there  are  many  parameters  which  may  affect  the  DIAL 
measurements,  resulting  in  a  small  number  of  similar  measurements  for  each  grouping.  The 
DIAL  results  have  a  significant  amount  of  scatter,  and  an  empirical  decision  was  made  to 
consider  any  individual  values  which  fell  below  50%  or  above  150%  of  the  expected  CL  to  be 
outliers,  and  to  remove  these  values  from  the  analysis. 
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Table  3.  Summary  of  DIAL  Results 


Target 

Detection 

Techniqu 

e 

Target 

Gas 

Mode¬ 

locking 

Time 

Segment 

Number  of 

Measurement 

s 

Measured  / 
Expected  CL 
Ratio 

4KSBB 

DD 

SF6 

On 

Post-Mar 

04 

5 

0.99  +  0.15 

4KSBB 

DD 

SF6 

Off 

Post-Mar 

04 

6 

0.93  +  0.10 

4KSBB 

HD 

SF6 

On 

Post-Mar 

04 

6 

0.77  +  0.17 

4KSBB 

HD 

SF6 

Off 

Post-Mar 

04 

5 

0.78  +  0.18 

7KS  Berm 

DD 

SF6 

On 

Post-Mar 

04 

7 

0.68  +  0.24 

7KS  Berm 

DD 

SF6 

Off 

Post-Mar 

04 

8 

0.78  +  0.28 

7KS  Berm 

HD 

SF6 

On 

Post-Mar 

04 

7 

0.87  +  0.16 

7KS  Berm 

HD 

SF6 

Off 

Post-Mar 

04 

8 

0.76  +  0.14 

4KSBB 

DD 

NH3 

On 

Post-Mar 

04 

5 

0.98  +  0.06 

4KSBB 

DD 

NH3 

Off 

Post-Mar 

04 

5 

0.93  +  0.10 

4KSBB 

HD 

NH3 

On 

Post-Mar 

04 

4 

0.78  +  0.26 

4KSBB 

HD 

NH3 

Off 

Post-Mar 

04 

3 

0.59  +  0.14 

7KS  Berm 

DD 

NH3 

On 

Post-Mar 

04 

5 

0.79  +  0.18 

7KS  Berm 

DD 

NH3 

Off 

Post-Mar 

04 

7 

0.82  +  0.16 

7KS  Berm 

HD 

NH3 

On 

Post-Mar 

04 

7 

0.83  +  0.20 

7KS  Berm 

HD 

NH3 

Off 

Post-Mar 

04 

6 

0.84  +  0.34 

4KSBB 

DD 

SF6  (Low) 

On 

Pre-Mar  04 

5 

0.82  +  0.08 

4KSBB 

DD 

SF6  (Low) 

Off 

Pre-Mar  04 

5 

0.80  +  0.12 

4KSBB 

HD 

SF6  (Low) 

On 

Pre-Mar  04 

5 

0.55  +  0.24 

4KSBB 

HD 

SF6  (Low) 

Off 

Pre-Mar  04 

5 

0.59  +  0.16 
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4KSBB 

DD 

SF6  (High) 

On 

Pre-Mar  04 

5 

0.80  +  0.09 

4KSBB 

DD 

SF6  (High) 

Off 

Pre-Mar  04 

5 

0.77  +  0.04 

4KSBB 

HD 

SF6  (High) 

On 

Pre-Mar  04 

5 

0.59  +  0.14 

4KSBB 

HD 

SF6  (High) 

Off 

Pre-Mar  04 

5 

0.61+0.11 

7KS  Berm 

DD 

SF6  (Low) 

On 

Pre-Mar  04 

3 

0.81+0.08 

7KS  Berm 

DD 

SF6  (Low) 

Off 

Pre-Mar  04 

3 

0.71+0.12 

7KS  Berm 

HD 

SF6  (Low) 

On 

Pre-Mar  04 

3 

0.84  +  0.15 

7KS  Berm 

HD 

SF6  (Low) 

Off 

Pre-Mar  04 

3 

0.67  +  0.25 

7KS  Berm 

DD 

SF6  (High) 

On 

Pre-Mar  04 

3 

0.70  +  0.04 

7KS  Berm 

DD 

SF6  (High) 

Off 

Pre-Mar  04 

3 

0.68  +  0.13 

7KS  Berm 

HD 

SF6  (High) 

On 

Pre-Mar  04 

3 

0.82  +  0.16 

7KS  Berm 

HD 

SF6  (High) 

Off 

Pre-Mar  04 

3 

0.64  +  0.17 

The  average  measured-to-expected  CL  ratios  are  shown  in  Figure  31  for  the  different  data 
groupings,  and  the  standard  deviation  of  the  ratios  are  shown  in  Figure  32.  The  data  points 
are  again  slightly  offset  on  the  x-axis  for  the  different  groupings  to  reduce  clutter  from  overlap 
of  the  symbols.  As  noted  previously,  definitive  conclusions  probably  should  not  be  drawn 
from  this  data  because  the  statistical  significance  is  questionable  because  of  the  small  number 
of  data  points  included  in  each  grouping,  but  some  probable  conclusions  can  still  be  made: 

1)  The  pre-March  measured  CL  values  are  generally  lower  than  those  from  post-March. 

As  stated  previously,  this  is  believed  to  result  from  improvements  in  the  system  hardware  and 
the  experimental  alignment  and  procedures.  Because  of  this  effect,  it  is  generally  better  to  rely 
more  on  the  post-March  data  to  draw  conclusions  about  the  effects  of  the  different 
parameters  on  the  DIAL  accuracy. 

2)  The  direct  detection  measurements  from  the  4KS  BB  appear  to  give  the  most  accurate 
DIAL  results,  with  the  closest  match  to  the  expected  CL  value,  and  the  least  variability. 
Referring  to  Figure  1,  this  is  expected  since  in  the  HD/DD  DC  experiments  the  direct  detection 
system  is  expected  to  outperform  the  heterodyne  system  in  the  speckle-limited  regime.  This  is 
a  result  of  direct  detection  spatial  averaging  providing  more  speckle  mitigation  than  is 
achievable  with  the  heterodyne  spread  spectrum  operation  in  the  system  configurations  used 
for  the  HD/DD  DC  experiments.  It  should  be  noted  that  the  amount  of  speckle  mitigation  for 
the  heterodyne  system  could  be  significantly  and  easily  increased  by  designing  the  laser 
transmitter  to  include  more  frequency  modes  than  was  possible  with  the  specific  BBO  laser 
configuration  used  in  these  experiments. 

3)  The  direct  detection  system  tends  to  be  in  the  transition  region  between  speckle-  and 
signal-limited  operation  at  7  km,  while  the  heterodyne  system  is  still  in  the  speckle- limited 
region.  The  DIAL  measurement  accuracies  and  variabilities  are  generally  similar  for  direct 
detection  operation  at  7  km,  and  heterodyne  operation  at  both  4  km  and  7  km.  This  is  seen  for 
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both  SF6  and  NH3,  in  the  post-March  data.  A  possible  explanation  for  this  behavior  is  that  the 
heterodyne  system  is  operating  in  the  speckle-limited  regime  at  both  4  km  and  7  km,  while  the 
direct  detection  system  is  in  the  speckle-limited  regime  at  4  km,  but  is  entering  the  transition 
region  between  speckle-  and  signal-limited  operation  at  7  km  (see  Figure  1). 

4)  Modelocking  improves  heterodyne  DIAL  accuracy.  This  is  most  easily  seen  by  noting 
that  the  modelocked  heterodyne  results  from  7  km  are  closer  to  the  expected  CL  value,  and 
have  lower  variability  than  the  non-modelocked  results.  The  same  results  are  expected  from 
the  4  km  measurements,  but  the  4  km  data  does  not  conclusively  support  or  contradict  the 
expected  result. 

5)  Direct  detection  results  are  better  with  modelocking  off.  Since  spatial  averaging  is 
the  primary  speckle  mitigation  technique  for  direct  detection,  it  is  expected  that  the  increased 
laser  power  available  with  non-modelocked  operation  will  provide  the  best  results.  The 
results  shown  in  Figures  31  and  32  generally  support  this  conclusion.  A  possible  reason  that 
not  all  of  the  data  points  demonstrate  this  result,  is  that  it  appears  that  the  system  transmit 
pointing  jitter  is  reduced  when  using  the  modelocked  pulse  train,  which  improves  the  DIAL 
results. 


Average  of  Ratio 


Data  Group 

Figure  31.  Average  measured  CL  /  expected  CL  ratio. 


39 


Standard  Deviation  of  Ratio 


Figure  32.  Standard  deviation  of  measured  CL  /  expected  CL  ratio. 
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5.  Summary  and  Conclusions 

The  HD/DD  DC  experiments  were  designed  to  provide  a  direct,  simultaneous  comparison  of 
the  radiometric  and  chemical  detection  sensitivities  of  the  heterodyne  and  direct  detection 
receiver  techniques.  Many  parameters  affect  radiometric  and  DIAL  accuracy,  and  it  is  often 
not  clear  why  measurements  differ  from  the  expected  results.  Numerous  unresolved  issues 
exist  about  the  capabilities  and  characterizations  of  DIAL  systems,  especially  with  respect  to 
systems  utilizing  heterodyne  detection.  In  order  to  address  some  of  the  major  issues,  and  to 
provide  an  empirical  comparison  of  performance  when  the  effects  of  all  individual  factors  can 
not  be  determined,  an  experiment  technique  was  developed  that  provided  the  first  direct 
simultaneous  comparisons  of  heterodyne  and  direct  detection  DIAL  measurements. 

The  radiometric  results  showed  that  the  direct  detection  return  signal  was  within  a  factor  of  1 
to  2.5  of  the  expected  value,  while  the  heterodyne  return  signal  was  within  a  factor  of  6  to  20. 
Previous  specialized  direct  detection  radiometric  efforts  at  AFRL/DE  had  resulted  in 
agreement  factors  of  approximately  1  to  1.4.  The  premier  heterodyne  radiometric  efforts  by 
other  research  groups  have  demonstrated  agreement  factors  of  3  to  10.  The  level  of 
radiometric  agreement  achieved  in  the  HD/DD  DC  experiments  is  reasonable  in  comparison 
with  these  previous  efforts,  considering  that  the  primary  focus  of  the  experiments  was  on 
multiple  wavelength  DIAL  measurements,  which  significantly  increases  the  system  complexity 
with  respect  to  single  wavelength  measurements. 

A  primary  issue  concerning  the  accuracy  of  heterodyne  DIAL  measurements  is  the  sensitivity 
of  the  heterodyne  detection  process  to  speckle,  which  directly  affects  chemical  detection 
capability.  For  the  HD/DD  DC  measurements  a  speckle  mitigation  technique  utilizing  a  spread 
spectrum  (modelocked)  transmit  waveform  was  employed  and  studied.  The  heterodyne  DIAL 
results  showed  slightly  higher  variability  (within  a  factor  of  2,  as  determined  from  the  minimum 
detectable  absorption  level)  than  the  direct  detection  results,  but  were  much  better  than  the 
results  that  would  be  expected  from  a  heterodyne  system  not  employing  a  speckle  mitigation 
concept. 

The  measurements  also  indicate  that  the  relative  DIAL  performance  of  heterodyne  detection 
improves  with  increasing  range  in  comparison  with  direct  detection,  as  would  be  expected. 
The  ability  of  heterodyne  systems  to  maintain  speckle  limited  performance  at  ranges  where 
direct  detection  systems  become  signal  limited  allows  smaller  systems  to  be  used  for  the  same 
range  requirements,  and  allows  a  longer  standoff  range  to  be  achieved  with  similar  size 
systems. 

The  primary  conclusion  of  the  HD/DD  DC  experiment  series  is  that  heterodyne  DIAL  is  a 
feasible  and  demonstrated  technique,  and  has  definite  utility,  especially  for  long  standoff 
range  operation. 
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Appendix  A.  HD/DD  DC  Datasets 


Table  Al.  HD/DD  DC  Datasets 


Date 

(yymmdd) 

Used 

in 

Figure 

Filename 

Target 

Modelocke 

d 

Gas  /  System  Configuration 

030828 

Test 

15KS  BB 

No 

No  gas  cell 

030828 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030828 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030828 

Testlc 

4KS  berm 

No 

No  gas  cell 

030828 

Testld 

4KS  berm 

Yes 

No  gas  cell 

030828 

Test2a 

4KS  BB  (new) 

No 

No  gas  cell 

030828 

Test2b 

4KS  BB  (new) 

No 

Empty  cell /SF6/NH3 

030828 

Test2c 

4KS  BB  (new) 

Yes 

No  gas  cell 

030828 

Test2d 

4KS  BB  (new) 

Yes 

Empty  cell /SF6/NH3 

030828 

Test3a 

7KS  berm 

No 

No  gas  cell 

030828 

Test3b 

7KS  berm 

No 

Empty  cell /SF6/NH3 

030828 

Test3c 

7KS  berm 

Yes 

No  gas  cell 

030828 

Test3d 

7KS  berm 

Yes 

Empty  cell /SF6/NH3 

030828 

Test4a 

None 

No 

System  characterization  test 

030826 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030826 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030826 

Testlc 

4KS  berm 

No 

No  gas  cell 

030826 

Testld 

4KS  berm 

Yes 

No  gas  cell 

030826 

Testle 

4KS  berm 

No 

No  gas  cell 

030826 

Test2a 

4KS  BB  (new) 

No 

No  gas  cell 

030826 

Test2b 

4KS  BB  (new) 

No 

Empty  cell /SF6/NH3 

030826 

Test2c 

4KS  BB  (new) 

Yes 

No  gas  cell 

030826 

Test2d 

4KS  BB  (new) 

Yes 

Empty  cell /SF6/NH3 

030826 

Test3a 

7KS  berm 

No 

No  gas  cell 

030826 

Test3b 

7KS  berm 

No 

Empty  cell /SF6/NH3 

030826 

Test3c 

7KS  berm 

Yes 

No  gas  cell 

030826 

Test3d 

7KS  berm 

Yes 

Empty  cell /SF6/NH3 

030820 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030820 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030820 

Testlc 

4KS  berm 

No 

No  gas  cell 

030820 

Testld 

4KS  berm 

Yes 

No  gas  cell 

030820 

Test5a 

None 

No 

System  characterization  test 
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030807 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030807 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030807 

Testlc 

4KS  berm 

No 

No  gas  cell 

030807 

Testld 

4KS  berm 

Yes 

No  gas  cell 

030807 

Test2a 

4KS  BB  (new) 

No 

No  gas  cell 

030807 

Test2b 

4KS  BB  (new) 

No 

Empty  cell /SF6/NH3 

030807 

Test2c 

4KS  BB  (new) 

Yes 

No  gas  cell 

030807 

Test2d 

4KS  BB  (new) 

Yes 

Empty  cell /SF6/NH3 

030807 

Test3a 

7KS  berm 

No 

No  gas  cell 

030807 

Test3b 

7KS  berm 

No 

Empty  cell /SF6/NH3 

030807 

Test3c 

7KS  berm 

Yes 

No  gas  cell 

030807 

Test3d 

7KS  berm 

Yes 

Empty  cell /SF6/NH3 

030724 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030724 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030724 

Testlc 

4KS  berm 

No 

No  gas  cell 

030724 

Testld 

4KS  berm 

Yes 

No  gas  cell 

030724 

Testle 

4KS  BB  (new) 

No 

No  gas  cell 

030724 

Test  If 

4KS  BB  (new) 

Yes 

No  gas  cell 

030724 

Test2a 

7KS  berm 

No 

No  gas  cell 

030724 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/NH3 

030724 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

030724 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/NH3 

030612 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030612 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030612 

Testlc 

4KS  berm 

No 

No  gas  cell 

030612 

Testld 

4KS  berm 

Yes 

No  gas  cell 

030612 

Testle 

4KS  BB  (new) 

No 

No  gas  cell 

030612 

Test  If 

4KS  BB  (new) 

Yes 

No  gas  cell 

030612 

Test2a 

7KS  berm 

No 

No  gas  cell 

030612 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/NH3 

030612 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

030612 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/NH3 

030612 

Test3a 

15KS  BB 

No 

No  gas  cell 

030612 

Test3b 

15KS  BB 

No 

Empty  cell /SF6/NH3 

030612 

Test3c 

15KS  BB 

Yes 

No  gas  cell 

030612 

29 

Test3d 

15KS  BB 

Yes 

Empty  cell /SF6/NH3 

030604 

Testla 

4KS  BB  (old) 

Out 

No  gas  cell 

030604 

Testlb 

4KS  berm 

Out 

No  gas  cell 

030604 

Testlc 

4KS  BB  (new) 

Out 

No  gas  cell 
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030604 

Testld 

4KS  BB  (new) 

Out 

Empty  cell /SF6/NH3 

030604 

Test2a 

7KS  berm 

Out 

No  gas  cell 

030604 

Test2b 

7KS  berm 

Out 

Empty  cell /SF6/NH3 

030529 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030529 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030529 

Test2a 

4KS  berm 

No 

No  gas  cell 

030529 

Test2b 

4KS  berm 

Yes 

No  gas  cell 

030529 

Test3a 

4KS  BB  (new) 

No 

No  gas  cell 

030529 

Test3b 

4KS  BB  (new) 

No 

Empty  cell /SF6/NH3 

030529 

Test3c 

4KS  BB  (new) 

Yes 

No  gas  cell 

030529 

Test3d 

4KS  BB  (new) 

Yes 

Empty  cell /SF6/NH3 

030529 

Test4a 

7KS  berm 

No 

No  gas  cell 

030529 

Test4b 

7KS  berm 

No 

Empty  cell /SF6/NH3 

030529 

Test4c 

7KS  berm 

Yes 

No  gas  cell 

030529 

Test4d 

7KS  berm 

Yes 

Empty  cell /SF6/NH3 

030527 

Testla 

None 

N/A 

System  char.  (HD  bias 
current) 

030527 

Testlb 

None 

N/A 

System  char.  (HD  bias 
current) 

030527 

Testlc 

None 

N/A 

System  char.  (HD  bias 
current) 

030527 

Testld 

None 

N/A 

System  char.  (HD  bias 
current) 

030522 

Testla 

4KS  BB  (new) 

No 

System  characterization 

030522 

Testlb 

4KS  BB  (new) 

Yes 

System  characterization 

030502 

Testla 

4KS  BB  (new) 

No 

System  characterization 

030502 

Testlb 

4KS  BB  (new) 

Yes 

System  characterization 

030416 

Testla 

4KS  BB  (new) 

No 

No  gas  cell 

030416 

Testlb 

4KS  BB  (new) 

No 

Empty  cell /SF6/NH3 

030416 

Testlc 

4KS  BB  (new) 

Yes 

No  gas  cell 

030416 

Testlaa 

4KS  BB  (new) 

No 

No  gas  cell 

030416 

Testld 

4KS  BB  (new) 

Yes 

Empty  cell /SF6/NH3 

030416 

Test2a 

7KS  berm 

No 

No  gas  cell 

030416 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/NH3 

030416 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

030416 

27 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/NH3 

030415 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030415 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

030415 

Testlc 

4KS  BB  (new) 

Yes 

No  gas  cell 
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030415 


030415 


030415 


030401 


030401 


030401 


030401 


030401 


030401 


030401 


030401 


Testld 


Testle 


Test  If 


Testla 


Testlb 


Testle 


Testld 


Testle 


Test  If 
Test2a 
Test2b 


Test2c 


Test2d 


Test3a 


Test3b 

Test3c 

Test3d 


4KS  BB  (new) 


4KS  BB  (new) 


4KS  BB  (new) 


4KS  BB  (old) 


4KS  BB  (old) 


4KS  BB  (old) 


4KS  BB  (old) 


4KS  berm 


4KS  berm 
7KS  berm 
7KS  berm 


7KS  berm 


7KS  berm 


6KSCC 


6KSCC 

6KSCC 

6KSCC 


Empty  cell /SF6/NH3 


_ No  gas  cell _ 


Empty  cell /SF6/NH3 


No  gas  cell 


Empty  cell /SF6/SF6 


No  gas  cell 


Empty  cell  /SF6/SF6 


No  gas  cell 


_ No  gas  cell _ 

_ No  gas  cell _ 

Empty  cell /SF6/SF6 


No  gas  cell 


Empty  cell  /SF6/SF6 


_ No  gas  cell _ 


Empty  cell  /SF6/SF6 

_ No  gas  cell _ 

Empty  cell /SF6/SF6 


No  gas  cell 


_ No  gas  cell _ 

_ No  gas  cell _ 

No  gas  cell 


_ No  gas  cell _ 

Empty  cell /SF6/SF6 

_ No  gas  cell _ 

Empty  cell /SF6/SF6 


No  gas  cell 


Empty  cell /SF6/SF6 


No  gas  cell 


Empty  cell  /SF6/SF6 

_ No  gas  cell _ 

Empty  cell /SF6/SF6 


No  gas  cell 


Empty  cell /SF6/SF6 


_ No  gas  cell _ 


Empty  cell  /SF6/SF6 

_ No  gas  cell _ 

Empty  cell /SF6/SF6 


Empty  cell /SF6/SF6 


Empty  cell  /SF6/SF6 


030321 


32. 

32: 


Testla 


Testlb 

Testle 

Test2a 


6KSCC 


6KSCC 
6KSCC 
4KS  BB  (old) 


4KS  BB  (old) 
4KS  BB  (old) 
4KS  BB  (old) 
4KS  BB  (old) 


4KS  BB  (old) 


4KS  BB  (old) 


4KS  BB  (old) 


4KS  BB  (old) 
4KS  BB  (old) 
4KS  BB  (old) 


4KS  BB  (old) 


4KS  BB  (old) 


7KS  berm 


7KS  berm 
7KS  berm 
7KS  berm 


030224 

030224 

030224 

030224 


030224 


030224 


030224 


030224 

030224 

030224 


030224 


030224 


030224 


030224 

030224 

030224 


Testla 

Testlb 

Testle 

Testld 


Test2a 


Test2b 


Test2c 


Test2d 

Test3a 

Test3b 


Test3c 


Test3d 


Test4a 


Test4b 

Test4c 

Test4d 


030204 


030204 


Testla 


Testlb 


4KS  BB  (old) 


4KS  BB  (old) 
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030204 

Testlc 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

030204 

Testld 

4KS  BB  (old) 

Yes 

Empty  cell  /SF6/SF6 

030204 

Testle 

4KS  BB  (old) 

No 

No  gas  cell 

030204 

Test  If 

4KS  BB  (old) 

Yes 

No  gas  cell 

030204 

Test2a 

7KS  berm 

No 

No  gas  cell 

030204 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/SF6 

030204 

Test2c 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

030204 

Test2d 

7KS  berm 

Yes 

No  gas  cell 

030129 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030129 

Testlb 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

030129 

Testlc 

4KS  BB  (old) 

Yes 

No  gas  cell 

030129 

Testld 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

030129 

Test2a 

7KS  berm 

No 

No  gas  cell 

030129 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/SF6 

030129 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

030129 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

030123 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

030123 

Testlb 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

030123 

Testlc 

4KS  BB  (old) 

Yes 

No  gas  cell 

030123 

Testld 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

030123 

Test2a 

7KS  berm 

No 

No  gas  cell 

030123 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/SF6 

030123 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

030123 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

021218 

Testla 

4KS  BB  (old) 

No 

No  gas  cell 

021218 

Testlb 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

021218 

Testlc 

4KS  BB  (old) 

Yes 

No  gas  cell 

021218 

Testld 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

021218 

Test2a 

7KS  berm 

No 

No  gas  cell 

021218 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/SF6 

021218 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

021218 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

021218 

Test3a 

15KS  BB 

No 

No  gas  cell 

021218 

Test3b 

15KS  BB 

No 

Empty  cell /SF6/SF6 

021218 

Test3c 

15KS  BB 

No 

No  gas  cell 

021218 

28  Test3d 

15KS  BB 

No 

Empty  cell /SF6/SF6 

021218 

Test3e 

15KS  BB 

Yes 

No  gas  cell 

021218 

Test3f 

15KS  BB 

Yes 

Empty  cell /SF6/SF6 

021217 

Testla 

15KS  BB 

No 

No  gas  cell 
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021217 

Test2a 

4KS  BB  (old) 

No 

No  gas  cell 

021217 

Test2b 

4KS  BB  (old) 

No 

Empty  cell  /SF6/SF6 

021217 

Test2c 

4KS  BB  (old) 

Yes 

No  gas  cell 

021217 

Test2d 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020926 

Testla 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

020926 

Testlb 

4KS  BB  (old) 

No 

No  gas  cell 

020926 

Testlc 

4KS  BB  (old) 

Yes 

No  gas  cell 

020926 

25 

Testld 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020926 

Testle 

4KS  BB  (old) 

Yes 

No  gas  cell 

020926 

Test  If 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020926 

Test2a 

7KS  berm 

Yes 

No  gas  cell 

020926 

Test2b 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

020926 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

020926 

26 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

020926 

Test2e 

7KS  berm 

No 

No  gas  cell 

020926 

Test2f 

7KS  berm 

No 

Empty  cell /SF6/SF6 

020916 

Testla 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

020916 

Testlb 

4KS  BB  (old) 

No 

No  gas  cell 

020916 

Testlc 

4KS  BB  (old) 

Yes 

No  gas  cell 

020916 

Testld 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020916 

Test2a 

7KS  berm 

No 

No  gas  cell 

020916 

Test2b 

7KS  berm 

No 

Empty  cell /SF6/SF6 

020916 

Test2c 

7KS  berm 

Yes 

No  gas  cell 

020916 

Test2d 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

020821 

Testla 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020821 

Testlb 

4KS  BB  (old) 

Yes 

No  gas  cell 

020821 

Testlc 

4KS  BB  (old) 

No 

No  gas  cell 

020821 

Testld 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

Test2a 

7KS  berm 

Yes 

No  gas  cell 

mmmM 

Test2b 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

020821 

Test2c 

7KS  berm 

No 

No  gas  cell 

020821 

Test2d 

7KS  berm 

No 

Empty  cell /SF6/SF6 

020806 

Testla 

4KS  BB  (old) 

Yes 

No  gas  cell 

020806 

Testlb 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020806 

Testlc 

4KS  BB  (old) 

No 

No  gas  cell 

020806 

Testld 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

020717 

Testla 

4KS  BB  (old) 

No 

No  gas  cell,  10P20  only 

020717 

Testlc 

4KS  BB  (old) 

No 

No  gas  cell,  10P20  only 

020717 

Testld 

4KS  BB  (old) 

Yes 

No  gas  cell,  10P20  only 
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020717 


020717 


020717 


020717 

020717 

020717 


020717 


020717 


020711 


020711 


020711 


020711 


020711 


020711 


020710 

020710 

020710 


020710 


020710 


020710 


020710 

020710 

020710 


020710 


020703 

020703 

020703 


020703 


020703 


020703 


020703 

020703 

020703 


020703 


020703 


020703 


020529 


Test2a 


Test2b 


Test2c 


Test2d 

Test3a 

Test3b 


Test4a 


Test4b 


Testla 


Testlb 


Test2a 


Test2b 


Test2c 

Test2d 

Test3a 


Test3b 


4KS  BB  (old) 

No 

Empty  cell  /SF6/SF6 

4KS  BB  (old) 

No 

No  gas  cell 

4KS  BB  (old) 

Yes 

No  gas  cell 

4KS  BB  (old) 

Yes 

Empty  cell  /SF6/SF6 

7KS  berm 

No 

No  gas  cell,  10P20  only 

7KS  berm 

Yes 

No  gas  cell,  10P20  only 

7KS  berm 

No 

No  gas  cell 

7KS  berm 

Yes 

No  gas  cell 

4KS  BB  (old) 

No 

No  gas  cell,  10P20  only 

4KS  BB  (old) 

Yes 

No  gas  cell,  10P20  only 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

4KS  BB  (old) 

No 

No  gas  cell 

4KS  BB  (old) 

Yes 

No  gas  cell 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

None 

No 

System  characterization  test 

None 

Yes 

System  characterization  test 

4KS  BB  (old) 

No 

No  gas  cell,  10P20  only 

4KS  BB  (old) 

Yes 

No  gas  cell,  10P20  only 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

4KS  BB  (old) 

No 

No  gas  cell 

4KS  BB  (old) 

Yes 

No  gas  cell 

7KS  berm 

No 

No  gas  cell 

7KS  berm 

Yes 

No  gas  cell 

7KS  berm 

No 

Empty  cell /SF6/SF6 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

4KS  BB  (old) 

No 

No  gas  cell,  10P20  only 

4KS  BB  (old) 

Yes 

No  gas  cell,  10P20  only 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

4KS  BB  (old) 

No 

Empty  cell /SF6/SF6 

4KS  BB  (old) 

No 

No  gas  cell 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

7KS  berm 

No 

No  gas  cell 

7KS  berm 

Yes 

No  gas  cell 

7KS  berm 

No 

Empty  cell /SF6/SF6 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

7KS  berm 

No 

No  gas  cell 

7KS  berm 

Yes 

No  gas  cell 

4KS  berm 

Yes 

Empty  cell /SF6/SF6 

Test3b 


Test3c 


Test3d 


Test4a 

Test4b 

Test4c 


Test4d 


Testla 

Testlb 

Test2a 


Test2b 


Test2c 


Test2d 


Test3a 

Test3b 

Test3c 


Test3d 


Test4a 


Test4b 


Testla 


49 


020529 

Testlb 

4KS  berm 

Yes 

Empty  cell /SF6/SF6 

020513 

Testl 

4KS  berm 

Yes 

Empty  cell  /SF6/SF6 

020513 

Test2 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

020513 

Test3 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

020507 

Testl 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020507 

Test2 

7KS  berm 

Yes 

Empty  cell /SF6/SF6 

020506 

Testl 

4KS  BB  (old) 

Yes 

No  gas  cell 

020506 

Test2 

4KS  BB  (old) 

Yes 

Empty  cell /SF6/SF6 

020506 

Test3 

7KS  berm 

Yes 

SF6  (steady  state) 

020506 

Test4 

4KS  BB  (old) 

Yes 

SF6  (steady  state) 

020506 

Test5 

4KS  BB  (old) 

Yes 

No  gas  cell 

020506 

Test7 

4KS  BB  (old) 

Yes 

No  gas  cell 
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Appendix  B.  HD/DD  DC  DIAL  Results 


VHD/DD  DIAL  Comparison  Test  030828_T3D 
_ (SF6  +  NH3,  ML  On,  7KS  Berm,  12”  HD  Rx)  3^ 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.79  ±  1 .82  (DD  SF6) 
0.22  ±  2.15  (HDX3SF6) 

-20  ±96(DD  NH3) 

6  ±140  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

2.00  ±1.64  (DD) 

7.82  ±2.53  (HDx3) 

0.9  (DD) 
3.4  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

232  ±  82  (DD) 

319  ±254  (HDx3) 

5.1  (DD) 

7.0  (HDx3) 

10.0 

030828_000T3D_CL05s  (12"  HD  Rx) 
Modelocker  On;  7KS  Berm;  SF6/NH3 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  030828_T3B 
_ (SF6  +  NH3,  ML  Off,  7KS  Berm,  12”  HD  Rx) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.1 9  ±  1 .25  (DD  SF6) 
0.25  ±  3.36  (HDx3SF6) 

-11  ±  101  (DD  NH3) 

26  ±  233  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

2.38  ±0.95  (DD) 

9.55  ±1.48  (HDx3) 

1.0  (DD) 

4.2  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

81  ±  108  (DD) 

178  ±244  (HDx3) 

1 .8  (DD) 

3.9  (HDx3) 

10.0 

030828_000T3B_CL05s  (12"  HD  Rx) 
Modelocker  Oil;  7KS  Berm;  SF6/NH3 


Sequence  # 
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IH 
I  CP  U 


HD/DD  DIAL  Comparison  Test  030828_T2D 
(SF6  +  NH3,  ML  On,  4KS  New  BB,  12”  HD  Rx) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-1 .74  ±  3.67  (DD  SF6) 
-0.28  ±  3.67  (HDX3SF6) 

-82  ±  185  (DD  NH3) 

9  ±257  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

12.89  ±  3.29  (DD) 

7.53  ±2.41  (HDx3) 

5.6  (DD) 

3.3  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

211  ±205  (DD) 

***  ±***  (HDx3) 

4.6  (DD) 

***  (HDx3) 

10.0 

***  Signal  drop  -out  on  heterodyne 

030828_000T2D_CL05s  (12"  HD  Rx) 


Modelocker  On;  4KS  BB  (New);  SF6/NH3 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test  030828_T2B 
(SF6  +  NH3,  ML  Off,  4KS  New  BB,  12”  HD  Rx) 
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VHD/DD  DIAL  Comparison  Test  030826_T3D 
(SF6  +  NH3,  ML  On,  7KS  Berm,  12”  HD  Rx) 


t5 

CD 

5 


P 


| 

I 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.15  ±0.94  (DD  SF6) 
-0.01  ±  2.65  (HDx3  SF6) 

-7  ±106  (DD  NH3) 

13  ±  160  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

5.95  ±1.52  (DD) 

9.49  ±1.27  (HDx3) 

2.6  (DD) 

4.2  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

4  ±  83  (DD) 

375  ±295  (HDx3) 

0.1  (DD) 

8.2  (HDx3) 

10.0 

030826_000T3D_CL05s  (12"  HD  Rx) 
Modelocker  On;  7KS  Berm;  SF6/NH3 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  030826_T3B 
(SF6  +  NH3,  ML  Off,  7KS  Berm,  12”  HD  Rx) 
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VHD/DD  DIAL  Comparison  Test  030826_T2D 
(SF6  +  NH3,  ML  On,  4KS  New  BB,  12”  HD  Rx) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.11  ±0.82  (DD  SF6) 
-0.05  ±  1.92  (HDx3SF6) 

-4  ±  50  (DD  NH3) 

26  ±  224  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

9.04  ±0.57  (DD) 

1 0.88  ±  1.15  (HDx3) 

4.0  (DD) 
4.8  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

426  ±  69  (DD) 

234  ±247  (HDx3) 

9.3  (DD) 

5.1  (HDx3) 

10.0 

030826_000T2D_CL05s  (12"  HD  Rx) 
Modelocker  On;  4KS  BB  (New);  SF6/NH3 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  030826_T2B 
(SF6  +  NH3,  ML  Off,  4KS  New  BB,  12”  HD  Rx) 
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VHD/DD  DIAL  Comparison  Test  030807_T3D 
(SF6  +  NH3,  Modelocking  Off,  7KS  Berm) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.04  ±  1 .67  (DD  SF6) 
-0.07  ±  1 .69  (HDx3  SF6) 

-6  ±129  (DD  NH3) 

9  ±105  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

7.50  ±1.27  (DD) 

10.50  ±  1.17  (HDx3) 

3.3  (DD) 

4.6  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

332  ±138  (DD) 

346  ±153  (HDx3) 

7.3  (DD) 

7.6  (HDx3) 

10.0 

030807_000T3D_CL05s 
Modelocker  On;  7KS  Berm;  SF6/NH3 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  030807_T3B 
(SF6  +  NH3,  Modelocking  Off,  7KS  Berm) 
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HD/DD  DIAL  Comparison  Test  030807_T2D 
(SF6  +  NH3,  Modelocking  On,  4KS  New  BB) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.11  ±0.80  (DD  SF6) 
0.07  ±  1.49(HDx3  SF6) 

-4  ±  58  (DD  NH3) 

4  ±152  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

10.27  ±  0.53  (DD) 

7.90  ±1.63  (HDx3) 

4.5  (DD) 

3.5  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

457  ±  75  (DD) 

21 9  ±91  (HDx3) 

10.0  (DD) 
4.8  (HDx3) 

10.0 

030807_000T2D_CL05s 
Modelocker  On;  4KS  BB  (New);  SF6/NH3 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test  030807_T2B 
(SF6  +  NH3,  Modelocking  Off,  4KS  New  BB) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

0.05  ±  1.08  (DDSF6) 
0.02  ±  1.55(HDx3  SF6) 

-1  ±  66  (DD  NH3) 

18  ±  169  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

9.75  ±0.68  (DD) 

7.34  ±1.79  (HDx3) 

4.3  (DD) 

3.2  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

503  ±  64  (DD) 

105  ±130  (HDx3) 

11.0  (DD) 
2.3  (HDx3) 

10.0 

030807_000T2B_CL05s 
Modelocker  Off;  4KS  BB  (New);  SF6/NH3 
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HD/DD  DIAL  Comparison  Test  030724_T2D 
(SF6  +  NH3,  Modelocking  On,  7KS  Berm) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

0.1 7  ±  1.02(DDSF6) 
0.07  ±  1.67  (HDx3  SF6) 

5  ±  52  (DD  NH3) 

10  ±  119  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

7.79  +1.38  (DD) 

9.42  +1.48  (HDx3) 

3.4  (DD) 

4.1  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

425  ±182  (DD) 

231  ±210  (HDx3) 

9.3  (DD) 

5.0  (HDx3) 

10.0 

030724_000T2D_CL05s 
Modelocker  On;  7KS  Berm;  SF6/NH3 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test  030724_T2B 
(SF6  +  NH3,  Modelocking  Off,  7KS  Berm) 
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VHD/DD  DIAL  Comparison  Test  03061 2_T3D 
(NH3  +  SF6,  ModelockerOn,  15KS  Blueboard) 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  03061 2_T3B 
„  (NH3  +  SF6,  Modelocker  Off,  15KS  Blueboard) 


i-swr 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5  -35 

1.01  ±  2.33  (DD  SF6) 
0.22  ±  2.21  (HDx3SF6) 

28  ±  129  (DD  NH3) 

6  ±205  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

26.73  ±  3.53  (DD) 

6.47  ±1.62  (HDx3) 

11.7  (DD)* 
2.8  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

1801  ±28  (DD) 

631  ±581  (HDx3) 

39.4  (DD)* 
13.8  (HDx3) 

10.0 

*  Low  signal  level,  error-induced  saturated  absorption  , 
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VHD/DD  DIAL  Comparison  Test  03061 2_T2D 
(SF6  +  NH3,  Modelocking  On,  7KS  Berm) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

0.05  ±  0.68  (DDSF6) 
0.01  ±  1.24(HDx3  SF6) 

1  ±  60  (DD  NH3) 

-3  ±132  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

9.10  ±0.61  (DD) 

6.77  ±1.72  (HDx3) 

4.0  (DD) 

3.0  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

41 0±  45  (DD) 

397  ±  90  (HDx3) 

9.0  (DD) 
8.7  (HDx3) 

10.0 

03061 2_000T2D_CL05s 
Modelocker  On;  7KS  Berm;  SF6/NH3 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  03061 2_T2B 
(SF6  +  NH3,  Modelocking  Off,  7KS  Berm) 
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VHD/DD  DIAL  Comparison  Test  030604_T2B 
(NH3  +  SF6,  Modelocker  Out,  7KS  Berm) 


Sequences 

Measured  CL 

(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5  -35 

0.03  ±  0.60  (DDSF6) 
0.16  ±  2.01  (HDx3  SF6) 

-1  ±  55  (DD  NH3) 

2  ±155  (HDx3  NH3) 

Insertion 
#1  (NH3) 

50-65 

426  ±  67  (DD) 

462  ±149  (HDx3) 

9.3  (DD) 
10.1  (HDx3) 

10.0 

Insertion 
#2  (SF6) 

80-95 

8.00  ±0.86  (DD) 

6.19  ±4.13  (HDx3) 

3.5  (DD) 

2.7  (HDx3) 

4.5 

030604_000T2B_CL05s 


Modelocker  Out;  7KS  Berm;  SF6/NH3 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  030604_T1D 
(SF6  +  NH3,  Modelocker  Out,  4KS  New  BB) 


62 


HD/DD  DIAL  Comparison  Test  030529_T4D 
(SF6  +  NH3,  Modelocking  On,  7KS  Berm) 


Sequences 

Measured  CL 

(ppm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

0.05  ±  1.72  (DDSF6) 
-0.16  ±  2.01  (HDx3  SF6) 

0  ±125  (DD  NH3) 

9  ±148  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

8.96  ±2.70  (DD) 

1 1 .28  ±  0.96  (HDx3) 

3.9  (DD) 

4.9  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

203  ±  95  (DD) 

457  ±119  (HDx3) 

4.4  (DD) 
10.0  (HDx3) 

10.0 

030529_000T4D_CL05s 
Modelocker  On;  7KS  Berm;  SF6/NH3 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test  030529_T4B 
(SF6  +  NH3,  Modelocking  Off,  7KS  Berm) 


tfcnrti-^fcKtra- 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

0.05  ±  0.76  (DDSF6) 
0.07  ±  2.36  (HDx3SF6) 

5  ±  88  (DD  NH3) 

-2  ±  238  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

10.31  ±  0.97  (DD) 

9.32  ±1.84  (HDx3) 

4.5  (DD) 

4.1  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

392  ±120  (DD) 

251  ±257  (HDx3) 

8.6  (DD) 

5.5  (HDx3) 

10.0 

030529_000T4B_CL05s 
Modelocker  Off;  7KS  Berm;  SF6/NH3 
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HD/DD  DIAL  Comparison  Test  030529_T3D 
(SF6  +  NH3,  Modelocking  On,  4KS  New  BB) 


■  Iet.-H-  **+.+LS-y- 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.12  ±0.78  (DD  SF6) 
0.03  ±  1.69(HDx3  SF6) 

-7  ±  78  (DD  NH3) 

13+  219  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-65 

9.56  ±1.03  (DD) 

7.66  ±1.43  (HDx3) 

4.2  (DD) 

3.3  (HDx3) 

4.5 

Insertion 
#2  (NH3) 

80-95 

418  ±128  (DD) 

339  ±196  (HDx3) 

9.1  (DD) 
7.4  (HDx3) 

10.0 

030529_000T3D_CL05s 
Modelocker  On;  4KS  New  BB;  SF6/NH3 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test  030529_T3B 
_ (SF6  +  NH3,  Modelocking  Off,  4KS  New  BB) 
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VHD/DD  DIAL  Comparison  Test  03041 6_T2Dm 
(SF6  +  NH3,  Modelocking  On,  7KS  Berm) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.09  ±  1 .22  (DD  SF6) 
-0.18  ±  1 .48  (HDx3  SF6) 

-2  ±  88  (DD  NH3) 

-6  ±145  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-95 

7.53  ±2.18  (DD) 

7.25  ±1.28  (HDx3) 

3.3  (DD) 

3.2  (HDx3) 

4.3 

Insertion 
#2  (NH3) 

80-95 

41 6  ±  51  (DD) 

516  ±160  (HDx3) 

9.1  (DD) 
11.3  (HDx3) 

10.0 

03041 6_000T2D_CL05s 
Modelocker  On;  7KS  Berm;  SF6/NH3 


Sequence  # 


VHD/DD  DIAL  Comparison  Test  03041 6_T2B 
(SF6  +  NH3,  Modelocking  Off,  7KS  Berm) 
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HD/DD  DIAL  Comparison  Test  03041 6_T1D 
(SF6  +  NH3,  Modelocking  On,  4KS  New  BB) 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.1 4  ±1.41  (DDSF6) 
0.02  ±  1.19  (HDx3  SF6) 

-0.14  ±1.41  (DD  NH3) 
0.02  ±1.19  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-95 

8.68  ±1.11  (DD) 

7.69  ±1.33  (HDx3) 

3.8  (DD) 

3.4  (HDx3) 

4.3 

Insertion 
#2  (NH3) 

80-95 

490  ±  68  (DD) 

51 7  ±  93  (HDx3) 

10.7  (DD) 
11.3  (HDx3) 

10.0 

03041 6_000T1  DCL05AS 
Modelocker  On;  4KS  New  BB;  SF6/NH3 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test  03041 6_T1  B 
(SF6  +  NH3,  Modelocking  Off,  4KS  New  BB) 
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HD/DD  DIAL  Comparison  Test  03041 5_T1F 
(SF6  +  NH3,  Modelocking  Off,  4KS  New  BB) 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test  03041 5_T1D 
(SF6  +  NH3,  Modelocking  On,  4KS  New  BB) 


■ 


Sequences 

Measured  CL 

(PPm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5  -35 

-0.06  ±  0.59  (DDSF6) 
0.13  ±  1.13(HDx3  SF6) 

2±  91  (DDNH3) 

0±  77  (HDx3  NH3) 

Insertion 
#1  (SF6) 

50-95 

4.38  ±0.53  (DD) 

5.18  ±1.67  (HDx3) 

1 .9  (DD) 

2.3  (HDx3) 

4.3 

Insertion 
#2  (NH3) 

80-95 

442  ±  96  (DD) 

342  ±138  (HDx3) 

9.7  (DD) 

7.5  (HDx3) 

10.0 

03041 5_000T1D_CL05s 
Modelocker  On;  4KS  New  BB;  SF6/NH3 
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HD/DD  DIAL  Comparison  Test  030224_T4D 
(Modelocking  On,  7KS  Berm,  0ND/5X) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

0.03  ±0.49  (DD) 
0.15  ±0.98  (HDx3) 

Insertion 

#1 

50-65 

4.92  ±  1 .00  (DD) 
6.96  ±  1 .40  (HDx3) 

2.2  (DD) 

3.0  (HDx3) 

3.0 

Insertion 

#2 

80-95 

9.23  ±0.76 (DD) 
13.71  ±1.23  (HDx3) 

4.0  (DD) 

6.0  (HDx3) 

6.0 

030224_000T4D_CL05s 
Modelocking  On;  7KS  Berm  (0  ND  /  5X) 


Sequence  # 


HD/DD  DIAL  Comparison  Test  030224_T4B 
(Modelocking  Off,  7KS  Berm,  0ND/5X) 
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V 


HD/DD  DIAL  Comparison  Test 
030224  T3D  (Modelocking  On,  4KS  BB,  1ND/5X) 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
030224 T3B  (Modelocking  Off,  4KS  BB,  1ND/5X) 
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HD/DD  DIAL  Comparison  Test 
030224T2D  (Modelocking  On,  4KS  BB,  1ND/10X) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

0.14  ±  1 .23  (DD) 
0.10  ±  1 .09  (HDx3) 

Insertion 

#1 

50-65 

6.09  ±0.72 (DD) 
5.13  ±0.91  (HDx3) 

2.7  (DD) 

2.2  (HDx3) 

3.0 

Insertion 

#2 

80-95 

10.63  ±0.96  (DD) 
10.38  ±1.59  (HDx3) 

4.7  (DD) 

4.5  (HDx3) 

6.0 

030224_000T2D_CL05s 
Modelocking  On;  4KS  Blueboard  (1  ND  /  10X) 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
030224 T2B  (Modelocking  Off,  4KS  BB,  1  ND/1  OX) 


r^kta  P 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

0.00  ±0.62 (DD) 
0.02  ±  1 .64  (HDx3) 

Insertion 

#1 

50-65 

5.60  ±0.60 (DD) 
2.65  ±  1 .96  (HDx3) 

2.5  (DD) 

1 .2  (HDx3) 

3.0 

Insertion 

#2 

80-95 

10.60  ±1.07  (DD) 
6.03  ±  1 .13  (HDx3) 

4.6  (DD) 

2.6  (HDx3) 

6.0 

030224_000T2B_CL05s 
Modelocking  Off;  4KS  Blueboard  (1  ND  /  10X) 
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Vo 


HD/DD  DIAL  Comparison  Test 
030224  TIDj [Modelocking  On,  4KS  BB,  3ND/10X) 


Sequences 

Measured  CL 

(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

0.03  ±0.53 (DD) 
2.02  ±  1 .46  (HDx3) 

Insertion 

#1 

50-65 

3.98  ±0.69 (DD) 
5.15  ±  1 .19  (HDx3) 

1 .7  (DD) 

2.3  (HDx3) 

3.0 

Insertion 

#2 

80-95 

9.68  ±0.89  (DD) 
9.31  ±1.55  (HDx3) 

4.2  (DD) 

4.1  (HDx3) 

6.0 

030224_000T1D_CL05s 
Modelocking  On;  4KS  Blueboard  (3  NDs  /  10X) 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
030224 T1  B  (Modelocking  Off,  4KS  BB,  3ND/10X) 
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VHD/DD  DIAL  Comparison  Test 

030129 T2D  (Modelocking  On,  7KS  Berm) 


Sequences 

Measured  CL 
(pprrt-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5-35 

-0.03  ±1.28  (DD) 
-0.06  ±1.01  (HDx3) 

Insertion 

#1 

50-65 

5.64  ±0.51  (DD) 
5.52  ±  1 .05  (HDx3) 

2.5  (DD) 

2.4  (HDx3) 

3.0 

Insertion 

#2 

80-95 

9.52  ±1.84  (DD) 
10.41  ±1.36  (HDx3) 

4.2  (DD) 

4.6  (HDx3) 

6.0 

0301 29_000T2D_CL05s 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
030129 T2B  (Modelocking  Off,  7KS  Berm) 
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HD/DD  DIAL  Comparison  Test 
0301 29_T1D  (Modelocking  On,  4KS  BB) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5-35 

0.03  ±0.75 (DD) 
-0.10  ±1.25(HDx3) 

Insertion 

#1 

50-65 

6.48  ±0.61  (DD) 
2.55  ±1.17  (HDx3) 

2.8  (DD) 

1.1  (HDx3) 

3.0 

Insertion 

#2 

80-95 

11.88  ±0.81  (DD) 
6.45  ±0.71  (HDx3) 

5.2  (DD) 

2.8  (HDx3) 

6.0 

0301 29_000T1  D_CL05s 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
030129_T1B  (Modelocking  Off,  4KS  BB) 
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VHD/DD  DIAL  Comparison  Test 

030123  T2D  (Modelocking  On,  7KS  Berm) 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
030123 T2B  (Modelocking  Off,  7KS  Berm) 
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HD/DD  DIAL  Comparison  Test 
030123  T1D  (Modelocking  On,  4KS  BB) 


Sequences 

Measured  CL 
(pprrt-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

0.07  ±0.92 (DD) 
0.03  ±  1 .12  (HDx3) 

Insertion 

#1 

50-65 

5.52  ±0.80  (DD) 
2.64  ±  1 .15  (HDx3) 

2.4  (DD) 

1 .2  (HDx3) 

3.0 

Insertion 

#2 

80-95 

12.29  ±0.53  (DD) 
8.46  ±  0.83  (HDx3) 

5.4  (DD) 
3.7  (HDx3) 

6.0 

0301 23_000T1 DCL05S 
Modelocking  On;  4KS  Blueboard 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
030123 T1B  (Modelocking  Off,  4KS  BB) 
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HD/DD  DIAL  Comparison  Test 
021218T3F  (Modelocking  On,  15KS  BB) 


Sequence  # 


V  HD/DD  DIAL  Comparison  Test 

021 21 8_T3D  (Modelocking  Off,  15KS  BB) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5  -35 

0.04  ±  1 .00  (DD) 
0.33  ±  1 .26  (HDx3) 

Insertion 

#1 

50-65 

4.37  ±  1 .36  (DD) 
5.92  ±  1 .24  (HDx3) 

1 .9  (DD) 

2.6  (HDx3) 

3.0 

Insertion 

#2 

80-95 

7.61  ±  1 .98  (DD) 
9.68  ±  1 .98  (HDx3) 

3.3  (DD) 

4.2  (HDx3) 

6.0 

021 21 8_000T3D_CL05s 
Modelocking  Off;  15KS  Blueboard 
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VHD/DD  DIAL  Comparison  Test 

020926 000T2F  (Modelocking  Off,  7KS  Berm) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 

Ap(psi) 

Expected 
Ap  (psi) 

Reference 

5-  35 

0.03  ±0.56  (DD) 
1.41  ±  3.65  (HDtl) 
0.34  ±2.11  (HDxl) 

" 

Insertion 

#1 

50-65 

4.75  ±0.47  (DD) 
9.14  ±  4.65  (HDtl) 
9.02  ±2.50  (HDxl) 

2.1  (DD) 
4.0(HDt1) 
3.9  (HDxl) 

3.0 

Insertion 

#2 

80-95 

9.28  ±0.92  (DD) 
14.57  ±2.22  (HDtl) 
13.02  ±  2.11  (HDxl) 

4.1  (DD) 
6.4  (HDtl) 
5.7  (HDxl) 

6.0 

020926_000T2F_CL05s 
Modelocking  Off;  7KS  Berm 


JL 


Sequence  # 


VHD/DD  DIAL  Comparison  Test 

020926 000T2D  (Modelocking  On,  7KS  Berm) 
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VHD/DD  DIAL  Comparison  Test 

020926 000T1D  (Modelocking  On,  4KS  BB) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

-0.01  ±0.51  (DD) 
0.13  ±1.19  (HDtl) 
-0.06  ±0.80(HDx1) 

Insertion 

#1 

50-65 

6.35  ±0.85 (DD) 
8.99  ±2.02  (HDtl) 
5.82  ±0.77  (HDxl) 

2.8  (DD) 

3.9  (HDtl) 
2.5  (HDxl) 

3.0 

Insertion 

#2 

80-95 

12.28  ±0.54  (DD) 
17.18±  2.88 (HDtl) 
11.73  ±1.55  (HDxl) 

5.4  (DD) 

7.5  (HDtl) 
5.1  (HDxl) 

6.0 

020926_000T1  D_CL05s 


Sequence  # 
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VHD/DD  DIAL  Comparison  Test 

02091 6 000T2D  (Modelocking  On,  7KS  Berm) 


Sequences 

Measured  CL 
(pprrt-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5-35 

0.00  ±  1 .15  (DD) 
1.17  ±3.44  (HDtl) 
-0.09  ±  1 .99  (HDxl) 

Insertion 

#1 

50-65 

4.79  ±  1 .32  (DD) 
2.37  ±2.53  (HDtl) 
5.03  ±1.22  (HDxl) 

2.1  (DD) 

1.0  (HDtl) 
2.2  (HDxl) 

3.0 

Insertion 

#2 

80-95 

8.59  ±0.91  (DD) 
6.30  ±4.90  (HDtl) 
9.96  ±4.13  (HDxl) 

3.6  (DD) 

2.8  (HDtl) 
4.3  (HDxl) 

6.0 

02091 6_000T2D_CL05s 
Modelocking  On,  7KS  Berm 


0  1  0  20  3  0  40  50  6  0  70  8  0  90  1  00 

Sequence  # 


VHD/DD  DIAL  Comparison  Test 

02091 6 000T2B  (Modelocking  Off,  7KS  Berm) 


Sequences 

Measured  CL 
(pprrt-m) 

Measured 
Ap  (psi) 

Expected 

Ap(psi) 

Reference 

5  -35 

-1.98  ±0.62 (DD) 
0.22  ±2.79  (HDtl) 
0.43  ±2.20  (HDxl) 

Insertion 

#1 

50-65 

2.73  ±  1 .23  (DD) 
3.84  ±2.26  (HDtl) 
3.11  ±3.1 2  (HDxl) 

1 .2  (DD) 

1.7  (HDtl) 
1.4  (HDxl) 

3.0 

Insertion 

#2 

80-95 

7.28  ±  1 .20  (DD) 

1.65  ±2.90  (HDtl) 
5.72  ±1.98  (HDxl) 

3.2  (DD) 

0.7  (HDtl) 
2.5  (HDxl) 

6.0 

02091 6_000T2B_CL05s 
Modelocking  Off,  7KS  Berm 
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VHD/DD  DIAL  Comparison  Test 

02091 6 000T1D  (Modelocking  On,  4KS  BB) 


Sequences 

Measured  CL 
(pprrt-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

-0.03  ±0.94  (DD) 
0.09  ±1.35  (HDtl) 
0.00  ±  1 .43  (HDxl) 

Insertion 

#1 

50-65 

5.16  ±0.68 (DD) 
2.94  ±1.48  (HDtl) 
2.72  ±0.82  (HDxl) 

2.3  (DD) 

1.3  (HDtl) 
1.2  (HDxl) 

3.0 

Insertion 

#2 

80-95 

10.10  ±0.92  (DD) 
10.40  ±  2.15(HDt1) 
7.47  ±1.59  (HDxl) 

4.4  (DD) 

4.6  (HDtl) 
3.3  (HDxl) 

6.0 

02091 6_000T1D_CL05s 
Modelocking  On 


Sequence  # 


VHD/DD  DIAL  Comparison  Test 

02091 6 000T1  A  (Modelocking  Off,  4KS  BB) 
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VHD/DD  DIAL  Comparison  Test 

020806  000T1D  (Non-Modelocked,  4KS  BB) 


Sequences 

Measured  CL 

(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

0.03  ±0.77  (DD) 
0.33  ±2.52  (HDtl) 
0.22  ±  1.36  (HDxx) 

Insertion 

#1 

50-65 

5.80  ±0.39 (DD) 
4.36  ±1.86  (HDtl) 
3.43  ±1.97  (HDxx) 

2.5  (DD) 

1.9  (HDtl) 
1.5  (HDxx) 

3.0 

Insertion 

#2 

80-95 

10.56  ±0.89  (DD) 
8.67  ±2.54  (HDtl) 
6.65  ±2.36  (HDxx) 

4.6  (DD) 

3.8  (HDtl) 

2.9  (HDxx) 

6.0 

020806_000T1  DCL05S 


Sequence  # 


HD/DD  DIAL  Comparison  Test 
020806 000T1B  (Modelocked,  4KS  BB) 
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VHD/DD  DIAL  Comparison  Test 

02071 7_000T2D  (Modelocked,  4KS  BB) 


Sequences 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psD 

Expected 

Ap(psi) 

Reference 

5-35 

0.08  ±0.71  (DD) 

0.11  ±1.82  (HDtl) 
0.04  ±1.37  ( HDxx) 

Insertion 

#1 

50-65 

5.77  ±0.55 (DD) 
5.62  ±1.67  (HDtl) 
6.08  ±1.71  (HDxx) 

2.5  (DD) 

2.5  (HDtl) 
2.7  (HDxx) 

3.0 

Insertion 

#2 

80-95 

9.92  ±0.59 (DD) 
10.40  ±  1.10 (HDtl) 
10.90  ±1.32  (HDxx) 

4.3  (DD) 

4.6  (HDtl) 
4.8  (HDxx) 

6.0 

02071 7_000T2D_CL05s 


Sequence  # 


VHD/DD  DIAL  Comparison  Test 
_ 020717_000T2A  (Non-Modelocked,  4KS  BB) 
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VHD/DD  DIAL  Comparison  Test 

„  020529  000001  (HD  Tx-Normalized,  4KS  BB) 


Sequences 
[Time  (hh:mm:ss)] 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 
Ap  (psi) 

Reference 

10  -25 

[1 :00:13  -1:00:32] 

0.26  ±  0.98  (DD) 
■0.34  ±  1.1 3  (HD) 

- 

" 

Insertion 

#1 

40  -55 

[1 :00:52  -1:01:11] 

7.95  ±  1.24  (DD) 
2.25  ±  1.42  (HD) 

3.5  (DD) 
1.0  (HD) 

30 

Insertion 

#2 

70  -85 

[1:01:31  -1:01:50] 

15.46  ±1.47  (DD) 
2.54  ±  0.81  (HD) 

6.8  (DD) 
1.1  (HD) 

60 

Steady 

State 

100-115 

[1 :07:10  -1:08:29] 

16.91  ±1.41  (DD) 
0.97  ±  0.89  (HD) 

7.4  (DD) 
0.4  (HD) 

60 

Insertion 
#3  (*sat) 

130-145 

[1 :07:49  -1:08:08] 

20.14  ±2.53  (DD) 
0.53  ±  1.16(HD) 

8.8  (DD) 
0.2  (HD) 

90 

Insertion 
#4  (*sat) 

160-175 

[1 :08:28 -1:08:47] 

17.43  ±2.74  (DD) 
1.35  ±  0.99  (HD) 

7.6  (DD) 
0.6  (HD) 

12.0 

25 

_  15 

- ah* 

0 

jP*'  ■*--  . . . — 

HD/DD  DIAL  Comparison  Test 
020529,000001  (HD  Not  Tx-Normalized,  4KS  BB) 


Sequences 
[Time  (hh:mm:ss)] 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 
Ap  (psi) 

Reference 

[1 

10  -25 

00:13  -1:00:32] 

0.26  ±  0.98  (DD) 
-0.01  ±  0.55  (HD) 

- 

" 

Insertion 

#1 

[1 

40  -55 

00:52  -1:01:11] 

7.95  ±  1.24  (DD) 
2.28  ±  0.70  (HD) 

3.5  (DD) 
1.0  (HD) 

38 

Insertion 

#2 

[1 

70  -85 

01:31  -1:01:50] 

15.46  ±1.47  (DD) 
2.98  ±  0.46  (HD) 

6.8  (DD) 
1.3  (HD) 

6.0 

Steady 

State 

[1 

100-115 

07:10  -1:08:29] 

16.91  ±1.41  (DD) 
3.62  ±  0.48  (HD) 

7.4  (DD) 
1.6  (HD) 

68 

Insertion 
#3  (*sat) 

[1 

130-145 

07:49  -1:08:08] 

20.14  ±2.53  (DD) 
3.98  ±  0.26  (HD) 

8.8  (DD) 
1.7  (HD) 

98 

Insertion 
#4  (*sat) 

[1 

160-175 

08:28-1:08:47] 

17.43  ±2.74  (DD) 
4.1 4  ±  0.28  (HD) 

7.6  (DD) 
1.8  (HD) 

12.0 
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HD/DD  DIAL  Comparison  Test 
02051 3 000001  (4KS  BB) 


Sequences 
[Time  (hh:mm:ss)] 

Measured  CL 
(ppm-m) 

Measured 
Ap  (psi) 

Expected 
Ap  (psi) 

Reference 

10  -25 

[1 :00:13  -1:00:32] 

0.29  ±  1.12  (DD) 
■0.02  ±  1.85  (HD) 

Insertion 

40  -55 

9.63  ±  1.72  (DD) 

4.2  (DD) 

38 

#1 

[1 :00:52  -1:01:11] 

5.56  ±  1.92  (HD) 

2.4  (HD) 

Insertion 

70  -85 

18.38  ±1.15  (DD) 
8.07  ±  1.40  (HD) 

8.0  (DD) 

68 

#2 

[1:01:31  -1:01:50] 

3.5  (HD) 
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Appendix  C.  MATLAB  Programs 


Section  Cl  -  DATREDUC.m  -  Direct  Detection  Data  Reduction 

Section  C2  -  DATREDUC_HD_X3.m  -  Heterodyne  Detection  Data  Reduction 

Section  C3  -  HDPulseSim.m  -  Heterodyne  Pulse  Simulation  (Optical  Field  Based) 
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Cl.  DATREDUC.m  -  Direct  Detection  Data  Reduction 


1  Hat  lab  program  ummc.M 

■Created  Cron  ENERGY,  P.  30  Sep  by  D.C.  Eenrt 

■  c<Ht#MTea  r*diiv«d  4<C«  ptrAHAttift  frvm  M-AbLt  digit iletf  Mavfrfotm  -dar* 
i  (Cl ,  Vpk,  tpk,  dte-£E,  Verna ,  VceEl 


Hvdi  fi  cat  lorn  - 

2  4  Oct  96  -  Transferred  iron  Eenfc  Ftfiac  tc  GLArs 

il  1  Jnh  Q7  |  lu-! '  i*  ;l  *  1  =  j r> - 1  *n-:!  tuff =*004  nlmJiJh-H  oh  t  L  ny  p]atp  rni’ii-HT 


than  on  average  plot 

G-7  Way  97  Corrected  Vpk  cals  jlarlcji  to  bo  redo  tranced  f  tooi  Vref  baseline t 
rather  than  0  V 

D7  Gun  97  -  Corrected  tpk,  cslru.lax.l-n  to  be  referenced.  Erom  start  of 
date  window  rather  than  start  of  nignal  window 


A  trays  work  in  h«-  i Ah  -  |jnif«»r**d.  coiuauj  r^Jhui.,  n*t  tow  rewat 

.  n  Id  pulae  I  ivarj.ee  iiomcoitally) 

.  L  la  sample  I  i varies  vertically) 

■  date  Li-r^-rn-o^l-jnnl 


ran  in  era  puls*  {Aanpl  I Lng|  time 
t  pulse  time 

lasibda  wavelength,  laser  line 
Npu.l*r  pulse  number 
Nseg  sequence  nmtber 


Processing 


Cprlntf  | 
fprintf I 
Cprlntr | 
f  P  i  i  n  r.  f  | 
tprintt I 


V,n> ) 

Kaeq  -  %d\n'rNwsi) 

Jflinea  =  U\rMalinia] 

Waanplea  ■  id'-.n k  H  Haample  a  |> 


-  taamplr 
run  -  riinCt) 
mem  ~  men  i  L  S 


foi  p-1  !fclppldoa 
blp*n)-n 

end 


for  n-1 :Mseq 
end 


i  n (invert  frnn  digitlSGt  umlta  tn  >>•> 
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t  pf 


s 


fprinlf 

Cprint? 

('  \fs 1 J 

(  ’■Digitizer  units  tc-  voltage 

conversion 

CprLnt? 

«  ' 

M  HMS  '  *  ( 1  f4  QS£$  0.1  □  a  | 

m 

Ms1  L 

Sprint £ 

i 1 

LftRE  (*<2/4OW>-Ii&0fl) 

LL3 

\r 1  h 

Ip Tint ; 

(  ' 

HP-DJL3  (1-1  | 

[H] 

\r.«  L 

ferine? 

<  1 

tfdlmar  (+11 

r  h  i 

\h'  L 

datalvolt  -  inputs  ' 

Brimroae  [-15 

LB]: 

:  i  ' 

iC 

data-- 

*l*fcll 

(.  1  detail: volt- 
data.  *  [1/4D9GI 

i  (d*tft2volt—  ! 

'HM  r  daetsivolt—  'n')} 

-D . 100 

1 1  ‘  1  t  d<S*r.*fu*lT..--  1 1  1 J  J 

dats-daza.  *  I 4 D 9« I  -1 . -U00 
a  lfisrl  I  (  (daaa.2volT.=- '  :H  ‘  |  :  d*£ata£’xolt==  ■  h  '  J  > 
date-data. 

els Eli  (  (daialvaJi.—  '  ?!  ”  |  1da.ta2volt==  '  k.'  J  f 
date-data 

e lsei  f  (  (dat  aZ volt—  '  3  ■  |  :  4data2voltP-  '  ,b '  >  > 

dAT-H-  diCj 

end 


J 


1  - - «*  Average  pulses  lor  eacii  laser  ws™  length  >>*- - 

avgdsta  -  i-ros-IFaenpleai-Niirceal 


lor  L=l: Klines  j 

for  w-l  =  Mse-g  ; 

*Bflinea+l  i 

avqdi"  a  ( \  r  l  j  -  jvqdn  ca  ( : ,  1 )  tdAT.ii  i :  r  p  |  p 

end 

avgdata  | :  r  1 1-  =o.vgdaia  (i ,  I]  /  K  seq  j 

4iif|id  ; 

I  - - Plot  OVftXlA^ftd  pjJiflfl  »? - 


M  TONS  -  4  p 

Mfrdw^-  NttM«  ; 

Heels  -  £1jc(  [Hiines-I.l  /SfrcrwsS +1  ; 

H Cools  =  HoolJi  p 

i£  rHiuies— 1}  ; 

Nfr«ra=2  p 

Hfc-sle-?  ; 

end  p 

If  [flwl«-c2]  p 

Hfrd-la-2,  p 

end  j 


•  ^.-figured1  tee Lt ii*n “ .  (1LM 1 1  ■  3 DO* Nfrale..  2 2 5 “ N £  z v«e  I 


NllneaCh  =  HlJLnos  t 

fdt  IdtrlfMl  ifuaeCh  ; 
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ChLinc  |LcJi5=lcJi 
find 

if  (Nlina.fi>£P) 
clear  C'nlirjt 
while  < IS 

lisp  I  1  '  | 
disp  |  1  *  *  * 

■disp  I  1 

l  31 sp  I '  "  | 

LineUsage-zercs  I  Lj  9*9  S 

LijieUi  j-gsi-LlnerjijisM  - 1 
LineErrcieifO  “  cer?  I  Line Jnde.Kr  IOC'S 

for  l-l;Hlincs 

LineUsagcILlnel  nrfexO  |1H  1M  =1? 

Li  lit  Choice 
pause 

minoa^h  m  G 
for  l-l:Hlines 

II  (LiDftUsage  ILiJialEdcKC  LI*  m=l> 

NLineaCh  -  NiineeCb  *  1 
ChLine  I  HI  inesCh.?  ”1 
and: 
end 

f  WUtf&Zh—Q} 

•licpf  1  1  1 

diop{  '  ***  Ho  lined  chosen  Cor  lisp  Lay  ■  *  *  1  S 

d.isp(  1  ‘■T'  cratl'  t-o  continue  mMmt  S 

cLia  p  i  '  1  l 

P*U*+ 

el ae l £  tHlln*  aCh>2  0  S 
dispC  '  1 ] 

disp{  1  *‘*J  HaKirrun  oi  10  lines  may  he  chosen  h  m  1  1  ] 

■lisp-;  '  ***  ^fat>  to  continue  ***'  j 

disp i  1  1  I 

pause 

Ml  -2.il 

break 
end 

end 
And 


Choose  lines  T.u  be  diploya-d  *k+l  ) 
■cret>  to  continue  *’■'  ) 


t 


J 


i 


i 


} 

t 


3  of 


3 


1C  CbattMIg=0|.  j 

Coe  1  ill- 1  =  flli  heia^li 

i  igure  <  ■  Foa  1 1  ion 1  ,  1 3Q*-«  U  ■  |  icii  -1  h  3  DG-fl  D*  (lcb-1  ?  r  iOC ,  AGO  j )  r 
and  : 

end  j 


*igr*l-  ; 

iC  [bat-rtifiij— -3 1-  j 

HUiUluZ.  f 

end  j 


Hyir-ai^fiflj  -  n' 
dacp  r  1 ) 

VJ 1/  E'qilC  If y-l.npUL 4  '  S-U r  f ..  .  P  ID.  ■  j" ,!  [HI  ;■  ■  1  f  '  £  ”  | 

if  Naur  Egwer-y—  1  f 1 5  I  I  eurfquery—  a  y '  1 1 
SURFDISP_B 

ipcj.ntf(’  Sti  "  5 

fptrinc.r(  nrf*e4  plot  [S]  Knr  I 

fpcintf(  ‘Print  (E4WI  [F/Pi»1.  \ti  1  J 

lpt  JLnif(  'Print  (‘-oloij  [£]  '  |> 

query-ljiput'j 'G£iiE.Lbue  I  i  =  1  i's'J 


;  :■  n i  r,nH  9 

J 


t 


f  H  (quajFy^  1  S  ! )  |  (qba^Y31-4  1  a  1  >  I1 
50RFDI5E_B 

else  i  E  ( (queEy==  '  t: '  1  |  quory==  '  p '  |  i 

EJ-rihC: 

el  ae  l  C  ( tq^lery,,^  1  PU r  |  ;  <qjuery— ' '  bv '  .1 ) 


print  t 

el  ae  i£  ( Lguery—  ‘ '  i  I  query—  1  c “  1 1  7 

print  -dpac  -PcoLdi  ,■ 

■1  i4  ? 

bze-aic  ? 

and  t 

end  j 

end  i 


i  E  (batctii  lp— fl  y  j 

IE  |*iiLCdlipilg==  1  V  '  J  s 

delete-JHlinea  +  U  : 

end  j 

end  i 


^  ViCiu'-at'liyplat  ; 

SigRefQuery  t 

far  Ldr= 1 :  HI  1 ncsCb  j 

fl-TUte  llcli}  f 

PI ofcSigRe;  : 

and  4 

Query  t 


nh  j  in  lyigreJT— 1 J 


;  Pl&i  U*;?; 
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if  |fiatehIlg=Q|  j 

■tilePldt,  r 

1  PlatSiglleE  4 

and  : 

Qwry  ; 

end  j 

if  fbatctrgXij""^  y  j 

Cor  lc*=l  EHlinaadh  4 

delete^lchl  r 

end  i 

end  ; 

elaa.r  NllnaaCh.  Ch.Llna  t 

1C  [hatcnri5=fl|  4 

fi-TUte  I  1  Position  ^  \  IQ .  ■■I 0  r  567  *  305 1  I  ‘ 

end  j 


5  pf 


3 


r  — - «<  Begin  integrated  signal  oa.lcoia.t itms  ■>»■ 

Pt_sl  -  l_ffl/dslt 
Ht_s2  '  t_aZ/delt 

dt_ri  -  t_rl/delt 
Ht_r2  =  t_il/delt 

lC  Nt_al  <-  U 
H”  1=  1 
end 

if  f H_<t2  >  Hi&dfrpL+e 
Nt_ai  "  Ivasnpiea 

and 

1C  Ht_rl  <-  □ 

Nt_£l-3 

end 

it  Mt_c2  >  Ks-anplia 
Hz  r2  -  HaanpLas 

end 

Ht_ft  -  fltrt-NtjHl-H 
Nt_r  -  Ht_r2-Ht_rl+i 


Cfaig  *  iaroa  i;Wa*qrH]  Lnaal-  4 

Uref  “  itroa  ^Hse^Hlinrea^  t 

Vref  =  zeros  i;Hseq,,trilttfca|  4 

V  -  SeroG  i  i 
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6  of 


s 


Sprint'  I  1  \Si "  }■  r 

tax  1*1  jHlinea  j 

Cprintn  1  Riivluc I pilj-  laaax  line  Id  of  Id.  |3Lap  1/3)  \:r 1  J 1. Wlirtsa)  ; 
£or  ai-liMse-g 

n=(n-Ll  *bfl^rreH-^l  j 

^.MndnS-^un  {-Si t*  I  NC_E  ] :  Nt_f£,  hJ )  ; 

foi  1  =Ht_il;bli  1 3  | 

Vref  Cmr  1 1  -Pee  £  lnH  1 )  -+dnta  [I,  n3  Fdel  t  r 

U  var=Uva  rt-  ^data  1  1  .  n  |  -awandat)  *  (dat  a  (1 ,  n)  -rreandat  [  *  del  t  ) 

efad  I 

VseC  lnH  1)-Prel  in,  1)  /  (Nt^r'deit)  i 

Veim  in,  l)E*qrt  (D«r/ £Ht  z*dal£))  j 

end 

end  j 


\  fprintf | y  \n' ! 

Mf  i-i  iHAlrm 

fprintf 4  1  Reducing  laser  lane  Id  of  Id  I  Step 
Co r  n=l:Hf»q 

a-  (r?-l  J  ""d  line ■&  +1 

L Vp-Jc  I  n ..  1 ) ,  piei  ndex  |  '-max  |  data  I  Ht_s  1 :  Ht_aZ ,  n] ) 
Vpk |Hj 1 ) -Vplc 3 i  Vra  f |nP 1 ) 

1  tpk  In..  1 3 -tsninple  Ipkindeit  j 

~pk  iBj  l)=tcaaiplE  Ipk-indc-ic+HL  il-li 

fn:  lH^t_H;NT._« 

rJsi  g  Crr.r  1 1  -Paig  In,  l^data  Uj  a)  'delt 
and 
end 
end 


t 

i 

_.'J>  r'  jl,Efl.inra) 


I  fprintf  4 1  Xu*  |  : 

Cor  i»l :  Wilrnijs  j 

fprintf -I1  Redmcirrg  laser  lane  %d  of  lid  I  Step  3/3)  \r‘  .  lHtfllneB) 
Cor  n=i:H±eq  j 

lfc^(IIrl)-lKlillM-H  1 

P  fmr  il  ■* Ifsag  imf  1 1  -  -J  fit_s  ‘deit  l  ‘  Vie  C  I  n,  1  ?  ; 

n-1  Cuad  {mP  1  |  =11  (re,,  1]  /  Vpl  In.  1  !■  f 

end 

and  j 

Sprint £ I  1 \n 1 h  ; 

daug^run  |U]  /Haag  ; 

dsdevatd  IPJ  t 

i  f  Lb%&dh  f  Jl*j— ri  t  ; 

first  -  1  ; 

uhl La  ill)  r 

s.f  <  first-— 1| 

cloi<3ftrP0*l[]'"9|  t 

IfiHUhe*--  Jl 
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tide  CERtype.,  '  File  ■  ,riL*,  1  .  aut  '  ,  nimZjtr<L^ivg*lC-'9|  , 

1  ^hV.Bl  1  +  'SD*v  -  1  rn^etrUflO^Utotov/Uh.1^  r  1  1  *■  I'll  ; 


elae 

rltla  [EEfttypa.,  1  =11.-  ’  rftl*]h 

end 

xlahel C 'P-Lsc  Mm  e  r '  i 
iO*b#l  t*>'  iSnV-d  1  | 

s-v™a}ii.=  r  Hmn-fiv<l|t  r  wra-av^  r  ymn-av ( 3 ^  r  VTTX'-'av  ( t  h 
11  (Tfttypa--  ■TJt*  > 

*ki  d  [  |  mn  xju  min  I  '.■ .  ynn  I  ynui  j  I 
end. 

tlr*t  -  Q 
end 

xaplac=-lnpnr.{  'Rascal*  plo:  fY  -  :■  ■’  |  M  3 1  ;  1  ,  ; '  ]i 

if  i  i  replot”  1  F 1 J  I  I  replct—  '  p  ’  H 
prJLDt 

eLoeif  ( Ireplo-t—  ‘ I  I  <pppiot—  '  i  ■*'  3  ? 
print 

!* rspplci t“  F  I  1  nc “  "  '  I ) 
print  -dpac  -  Pooler 
nla*  1C  ( (raplc-r-=  'V  l-'Hicplot--  '  y r  j ) 
b  leak 
eLae 

vri  l  p.^liipnt  i[  1  FjiLut  nau  Ji  m  ■ .  i .  '] 
unax-input  $  'Enter  new  K  max-  1  I 
yni.n=lnpue{ 'Enter  new  V  min.  1 l 
yriait-itfput  <  1  EjiL  et  njw  f  ‘I 

sxiiM  |jaj-r.  kit -ax  '/min  ymaxjl 
and 
end 
end 


t 


t 

J 


I 


J 


r 


I 


%Epint  *  input  |  hCcmvert  to  energy  (Y/jKIJ  $  t  sj  1  **  l  : 

IlCHEplnt^'Y-l  |  CEpInt='y'  J  >  t 

<  Rv- input ('Enter  conversion  factor  IV/W  or  V.a/Jlr  'f  ; 

1  Cirat  =  1  t 

I  uMl#  1 1  >  ; 

I  If  (first |  r 

I  pint  (Hn.,  D/Svh  f 


title f | ?Rt  ype , 

'  File  - 

Hililer'  Average  ■  '  ^nunZatr  IUflvq/RvJ 

"  fJP 

j  'Jbav  = 

“,  nmnlstr^ laa*1 Dsdev/ Davgl  i  r  |4JP])  j 

slfl.be  1  L'Fulse 

Number' ) 

i 

ylabel 1 'Energy 

') 

i 

CLrat  ■»■  0 

i 

'  end  t 

1  raplnt-lnput  CHaacila  plot  I  3F  A[6TS  ]  i  1  ► ps ' )  ) 

4  if  l  Ireplot'— 1 Y1  H  i Eepl^C"1  y 1 )  ^  break;  end  i 

i  Kmin-lnpijt i  ‘Enter  new  K  min:  1  i  ; 

1  jfflhu^lnput  \  'Entail,  nau  K  ifjxl  'J  } 

4  /nun-input  | 1  Enter  new  E  mm:  !>  ; 

l  ymax=input 4 'Enter  new  5f  max:  '  >  j 

*  4Xl*Ujflilh  ymifi  y#wi>!ll  i 


T  Df 
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t  and! 
'end 


tract.  =  '.SCO,  CO,  2 67 , 4C-3 ] 
t  HftJHi  L#t -liiy  r«  i 1  i^el  t  i«sh h 


ir*<:T  | 


p 

rJ  I \i  ■  tog r^LTi  ) 19 ,  r  mr.Zst  r  (Hpo-Lee  b  i- , 


IhlKtCCr/EUvq,  lOD'i 
-  -title  £  iTfitype,  1  Fixe  -  \fileP' 
l  "  pjlsca  total)  ' j j 
iKlfcb#  1 1 p  NdlftMi  1 1  led  3 1-ii'fi^l r  J 
*y  libel  I  '  Hunter  cl  Occuttencea 1 | 


ItrbiLedi  ; 

t  aacaut-Lnput rwratc  A3CEE  112c  lor  bis  Lour  an,  [5f/ [H] ) £ % ^**5  p 
I  if  ( (aBC-ou.L-’-"1 1 1  i*  laatcut’-—  'y 1 1  I  break;  and  t 

1  [iiequ  uina  ]  “hi  a  1 I UV  L~avq ,  1 U  Q )  p 

I  ri -wqu^  CT-yq'j Npu  1  i,hfSi  ; 

I  cave  ijuii.aac  imu  -arcii  -<fcrabie  ; 

l  fprlntfC  \ji  1 )  p 

I  rpcinLfd  File-s  fxequ-B&S  and  yyu-BBc  uilu.cn  ')  ; 

1  break  r  cud  p 

land  p 

it  IbatdiCLEp^h  j 

d*l-utu|l]  p 

t  delete  1 21  ; 

end  p 
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C2.  DATREDUC_HD_X3..m  -  Heterodyne  Detection  Data  Reduction 


M2HS04  J:l? 


I 

I  DATPE0\y:_UZ_X.3 
1 

Created  by  DCS  cron  DiegoFtacecsing_DC5r  origj.nn.lly  written  by 
i  Diego  ?l4rro:KL,  Din  Bckelkazip-lokei^  and  Dan  S-cn'c 

I 

'  14  Wot  iJ-2  -  Changed  conpotation  at  CRFjcI  to  use  Welch  per-iodoqraiis  with  equal  signal 

1  and  noln«  windows  in  or  das  to  match  raxoLutloBX  lor  pormall  ratio  n 

*  Idiviei-onl .  Fenane-  CNRk3  . 

1 

I — - - - - - - — 


\  DotiuI  tipi  Ming;  tirio  raw  data  Into  Its  lospcctlvo  wavelength . 
'  Ho  header  Cile  ia  uced  £or  thas  vecei-on. 

i  Have  length  .sequence  assumed  to  be  standard  13  wavelengths, 

■  slatting  with  lUFiS,  IOpIO,.  10P]£,.  .  ,  .  ( I0F3Q 


ibS  eloj»  All 
•TS  clear  all 
lb3  pack 


Tstart  «  cputine 


diep4 ' 


so  =  Input  I 

del?  “  L/Ga 

FPI  -  input  I 

Hseq  -  input  | 

Wlrnea  “  input  I 

dispr  1  > 

■SigRvfUOHryl 

Ht_cl  -  round  I  t_sl/deltJ-»-l 
Ml  s  1  -  pound,  [t  ai/doit) 
Ht_rl  “  round  I t_rl /dolt > +1 
Ht_r2  -  round  [t_r2 /del  L) 


inter  sampling  Itequenry  I h.2 ;■ 

inter  FRF  (He I  [IQ] 
inter  Nseq  1 101 1 
inlet  Nil re  a  : Li | 


;  r._Hiri._r.?,i_r.j 


TestFlq  “  G 
iC  LTestrlg— Of 
dlapC  1  f 

diept1  Feeding  in  Test  H  P""1 ..  txt  1  I 
dlspf  'f 

load  Te  stHDDD . tit 
IDS  plckfll* 

'■•DS  load  (together! 

? *  1  Fu  L  n  [.  R-  L  enq  (-h  !  t'-G  ff tMD&Eif 


;  Hornal  (Uf  or  tea-  (If  read- ia 
i  I  HncTToI  reed-ln 
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Pol  n  ta  I  &  rflfi  u  t=Tota  IB* i  n  t  s  /  Ha  e  q /H  1  ines;  j 

APnJji(,ePer5iio,t.-2(iQtH  i 

1  fct/n*Atd  ; 

■■  JoriTi  daza  array:  data  (paint ,,  sbnb,,  1  int  i 
daC-i-Ser^e  (folfltafar&hat,  rfflidr  HLliftiil  ; 

Cor  im  -  I;Jf=eq 

tOl  11  =  1  ifcflLMfl  ? 

Hpt  L-round ( ( [ran- 1 ) " Hli  ne  r^IL- 1 1 w Poin  teFerShot  + 1 }  ? 

MptI -sound  ( ( (nn- 1 J '  Mil nea-i-Ll  J  *PoiataPerS-liat  |  ? 

( ;  r  in4 1 1  i  -Tea  cADDD  <Hpt  L :  1 )  ; 

1 1 1  Hpt  1 --round  I  Rpt  1 II  ;  BE Ciiy-1  r  keyboard,.-  end,  ; 
l  IE | Hpt2- =r ound ( Hpt2 1 | i EEC lq=2 ;  keyfr::  u rd  r  and  ) 

end  ; 

end  f 


while  <  (I )  r  ■■----■™  -  ---  - 

MsBKpHaifcq-1 

dfl!  H-SBr^o  (PPlritePerSJlol.,  Nneq.r  Milner  I 
Cor  am  -  Lilfaeq 
fo:  11  -  liMUftW 
nn--l:iin-il  ’■Mlinee+11 
il(nn.?=  1 3  3iHS  |  j  nn=nntl|'  end 
Apt  L  -townd  i  Cjirt— 1  ]  T  tfei  r  l  v  PerShdC  ■  L I 
MptL-saund  (nn  *  PointaPeraho'c  J 
data  ( :  r  ms,  1 1  *  -7*a  t.HDDD  (Mpt  L :  1 J 

end 
end 
ehd 


f 


4  •■ 


O3-O0-2'8>  1'2A  nririjel  cv*  cection 


jfcip  bat  pulse  i  E 1 D 3  E  i- 


cltar  rertHDOL 
p*G=k 

elsft 

load  I-a  stHHDDb 
d eta  -  ada  ta 
clear  /data 
Nu^[  ■  1 V/ 

end 

t  lne’-l  /sa :  1  /  sa :  PnmtsPtcahot  /  sa 


f 


I  Past  read-in 
J 

t 


Hirpi! '  "  !• 
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■lisp-;'  starting  Noise  ?£D  c aLculati i>n  loop'  | 
diep<‘  H) 

.'■it':T.:i  -  naxCHt  #J?-Mt  iiltlPMl  K.'l-PJt  E 1-1-1  | 

NTTTJb*ar  -  rail  Ucg£  ITSpta  M 
HFIT  =  Z'NFET'baSE 

ior  nr  "  1:1 ;Hjeq 
Cor  11  =  L 1. 1 1  k-ll  1  rwr 

n.c-i-ae-™,to  I Nt_rl :  Nt_r2 ,  roi,  11 1 
l  [  Fun  i  r  r  |  "puel  Jii  |  nni be ,  Cl#  i  i ,  NFFI ,  i,i  i 
r Pnfi .  rftl-pwil$degr*mtneiLH*r  ;  |  ,hfft4  a-0 
P:  ih?ED  r : ,  roiH  ll  I  "Fun 
and 
end 

diep*'  Pone  TTcj.ee  PS"  calculation  l^op-1  J 

dispr  ■> 


t 


f 


! 


iiH+«  £:iiijmri  f  cat  iTitirF  H:'rc':Un  H-:\i  r  T'-hIl]! 


I  Rc-aErango  daia  into  corract  uavrElengcli  .sequence  OEdar 

dj.sp<"  Starting  spur  freqneDcj  loop'  I 
diapr  ■*> 


de:I_fn  =  fr.(2|  p  %  f  ir.l  spur  "requenclB-Ei 

H_ffil  -  roO  nd,  ( 1 75 .  Oet-  fn  L  ]  I  l  ; 

H_CnZ  -  round  I  I  E  l .  Jee-ln  [  1 1 J  /dEl_fjiJ  ? 


for  ran  *  IrNneq  ; 

Cor  11  =  li.KlJ.nE3  j 

[  ICY  r  H  |  n4«fi  I  rolp  1 1 1  1  -  n^fi  isJel  i  &+*  PSD  I  N _ 1 5"l  ]  =  H_l«2  ,  «nr  111  i  ; 

PpurFreqi'ninj,  1 1  ^  "  f n  IHxo±ek  Irai,  ll| +-W_Cni-»-l)  ? 

3-F  Tnda-H  (mrij  1 1  !■  =  nun  t  <11-1]  /If  lines  p 

; 

end  } 


qqepu.r  "  m  ed.  is  n  <  Epu r E r-eq ?  r 

tiqure-C 1  Posi-iT  j  '  H  |  50  IDO  400  300  |J  ; 

■plul  <  3F  XndSM ,  SpLiEfTC-q/  lofi ,  '  .  '  J  j 

tltl*  ( K  SgMt  Frequencies  <CEi£l  1  )  r 

av-axls  r  Hmn-au-<l>  ?  amK^avClfi  r  ym«av  <2  )  :  ynuc^av [ 4 p  j 

anS.a<rD  Hini-q-i- ]  ywm  I) 

figure ( 1 >: :n [300  4D0  4DC  300|]  f 

plot  <-Tq.spLir/le£ !  ? 

jiold  on, 

t-lul.  -IqqipDj;  / 1*6  .  '  to: 1  J  l 

title {  '’led jar.  Spur  Fcequency  IMJic!  r  | 

diSp<  '  tr  r  h  r  y  inn  ini  .  1  n  I  :  r;i»7r--:!  rfHi'ij  .m  i  !.h  :-  Inr"  \  r 
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■lisp  < '  ■>  r 

i  Rearranqercenx  assumes  first  wavelength  in  sequence  in  always  lowest  vavelengih.,  and 
%  fri-iSlW-Ht  IpU  f  fX4«q  JApOy 

■  Th.io“  out  pulse  a  prior  to  first  acrurcenct  «?£  rn  sin-.m  spur  frequency,,  “ill  cause  Ids  a 


l  of  one  sequence 

LYTij  JHiKlndeH  ]  -  ms* iqqacjrl 

for  rani- 1 ;  Hseq- 1 
Cor  h/wl=l  .mines 

h!ifl-  |  rail-  i  |  *  W]  j  n*:<  '  »  H  .  mAjilndH-x  -i 

mil- f  Loc-c  I  I  nnfl  -0 .  5)  /Hi  uvea }  -*i 
wwD=nnC-  [tuiu-1  J  -mines 
•Sate  I  r j mr.l .  nw  1  s -date  L - . jutO r ww&l 
end 
end 

dsxa< :  rHafr^r  il  -  FI 
Haeo  =  Nseq-1 

query-input  (  ’  Enter  <  set.*  tc  continue 
dalanMl:?] 


J  Ci-HJd  L  fHlti'iJ  jhidhlt  of  it  jiu  IrfHt 

f 

'  H  '  =  '  C ; 


clear  HoiseEE-D  j 

f-acir 


I  Calculate  FSDb 
diSp  i '  *  'f 

dacpC 1  Storting  FED  col  ruler Lon  loops '  J 

dispr  ■>> 

Nt_ff  -  Nt_s2-Nt_al-H 

WC-_r  - 

1C  [Nr  avHx  i) 

diep('  iS  Ht_s  «■  Ht_r  * 1  1  j 

fceylcard 

end 

Npx.ii  “  Nx_r 

HFTTPaae  -  ce  iHlcg£  IPptsh-«-l]i 
HEFT  -  2  "  HFF  TLrs  l: 

1  MLen  -  HFET 
Ht*n  *  NFFT/J4 

HorseFai  =  zeros  (NLen/Z-!-l,'Nseqi  mines | 
iighilPSd  -  BGfOB  (NLeh/2-O  .  HiMqrtfllj|i$:!sl 


t 


■ 


;t  lor  peiiodogienr 
■  I  foi  Ha:  eh.  oae lode gram  a 
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MecnrSD  =  i^ros  (HLer,/ Z+l.  Hscqi  MliJi-ifil 


F 


for  cm  *  l:Waeq 

fprlntr^  1  £*;  qua  cl  r  I.  1  \n'  r  r«Jij 

for  11  -  LrWlinta 

iLoi5e=fiata  |  Ht_c  1 :  im„  11J 

\  Pfin  r  fh  1  -pu^l-sh  I  nGi  ,  £  1 ,  H  >  N&j?ii  r  ok! 

i  L  Pnn  1 1  nj  -ptricdoqrair,  (noise  ■  !  I  .>  Wien ,  a-a ) 

JtoiEflp3fD£:  rmril|  =Pnr 
if  ilHviz  tPS  Dll,  imr  i  1 1  — D  i 

WalsePaJJi  1 1 1  [tvt.,  11 )  =Hd±«raD  1 2 ,  mn .,  1 1 } 

Li  iWG  i^tPE  D  (Hltn  /2  J-l ,  mr:r  il  1  — O I 

Mol;:  &P-SD  |  HLan/2i-l  r  run..,  1 1 J  =Hul  saPSD  |HLan/  2 ,  imi,  1 1  > 
end 

lapdijnj— ceil  lMt_o/Nt_r  l 

S-ig  PSD-  =e  cos  £WLeny  Z  + 1 ,  1-lplgiTr  I 

fP:  lipd  -  1  JNpiirjlUfl 

ejttra-counral(  EHpdqTr.s^  ttpta-Ht_a‘)  f  lUpdqms-l?  ^ 
yi=Ht_al^<Wpd-Lj  *  (Hp-LS-CKtra} 

Prf-Nl  ttipt-5-1 
3  lgnttlTdat  a  (H  1 :  H2  h  rn? ,  1 1 J 
|  Pi  ih,  f  :'i  j  -pw-H  leh  1  s  tgnal  B  r  ] ,  [  I ,  Mian,,  sa  j 
lPnnH  in  j-peciodogcan  Icignai,  [] ,  NL*nr  sal 
SigP3D(: ,  Hpi/  =?r.n 
eft<3 

E-ignalPEDi : j  inn,  ll]  ~mean  £ BigPSD'  J 


F 

J 


F 

;  k3 
r»  xz 

i  *  eliminate  diwi-Sc  by  ceiap 
F 

h  1  cLiminait  divide  by  zecos 
F 


F 

] 

1 

J 

r*  X2 

F 


Jtomtf'SD  £ ; ,  1 1  h  -SiqnalF&D  i : ,  imr  i  i  I  .  /EfOi-EtFSr- 1  r  ,  on,  LI  ] 


*  f-*t  qcj-2  !N^P172  S*  H.2  3-.pt.  thins  ffrr  hfcetfc 

’  nz  3-  (HaiseFa  D  |  qq- 1  r  ran,,  11  S'  -*-Ndi  aeFED  £qq,  miH  ll  |  +WaaaeP5Ii  i  qq+1  r  rm ,  1 1  f  h  /  J  ; 

1  WuniF,3biqqrnfrf.PLi;i-3ign.alPSD<:i;qr  ran,  1 1 K  /n=-3  / 

I  e-nd  ? 

1  m  3  -  | BnlaePED  ( 1 ,  cm,  LI  |  +JtolacFED  iZ ,  cm,  1 1 1  +-Ho1*ePBD  [3 ,  raiH  11)1/3  / 

\  *^nSJD<lr«F  lU-sl^^tPSCs(]  ,sfflP  1 J  i /n*3  t 

I  r.z 3-  I  Hoi aeFBD  (HrFT/5-lH  wwnr  11 !  +-Kd istPE D rHFFT/2j im,  11 1  +-SoiHePE  D LHFFT/2+1  _  nnir  11] )  /3  t 

I  HoniPSD(irFFr/2414!im.  HjaaigBaLPaZHUFFT/24-l^nii,  1.11 /nz3  ] 


end  j 

end  ? 

d 1 ap i '  ' > 

d.isp£'  C'ont  B'SZ  r-alc'j.atlen  1  c-aps 1  I 

dlapr  J 


J'Li-t  tlra  jiailwj,  signal  and  nolAv  spvztra..  and  rL&fnall sad  signal  ipaerrun 

YO  =  JDO  j 

aol  -  1  f 
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SG2  =  1 

iC  iHlJJ1*o>-3Dfr 

YD  -  TOD 

flci»l/3 

so2-£/-3 

end 

f 

t 

tmn  "■  tineili 

trix  -  Gall {tinai 

f 

f  ar  nrrl  :N;e  j 

f|'M  ntf  !  1  ■  ”,‘ :  ^ .  i  «■!  r . !  :h  f:  :i  \h '  rWH] 

tor  1 1-1: Hi lues 

p 

fi^uce  I  1 :  ■  \  inn '  ,  [£5-i-2  5*  l  II- 1  > "col  r  10-2  0  *  I  11-  i  J '  sol  f  GAO'  Ki ,  TOO  *ac?J  i 

Clg-uce  |  “BaaitiDJi^ ,  ?Z1+Z1'  ill-3  |  *acl,  YD-21*  |11-1|  *9  =  1,  ec-D^cl,  3Da*sc2]  * 


1  P-3  nr p ,  DineNuine  Hire  Index  ( 1  f  111 •  ;  1  1 

1  text  (nnx- □  -  IS*  ^i?y-  xmn  $ ,  yiu-] .  OS*  (ynx-  yim  |  j  ... 

i  j  F  :t  r  [a  i  m 

1 

subplot \4 r 1 P 1)  f  1 

plot  ( tine*le£.. date  1 :  H nTF.r  1 1  !< )  ; 

aet(gcflj  ' Fonts  ire ' ,  |3  | ]  p 

1  *at,(gdi,  '¥T^CKt4t»lrP  FI  1  ■ 

title  1 'Tine  Senes' ; 
uVrj>:  is  j  >:nri=avi;3|  j  >:n:t=av  (21  j  ynn=av  [21  t  ymx=avr4] 
■aaisUfl  20  -IrO  ]  - 0 1  >  r 

hold  on  p 

riots  igRail  t  % 

Tino  ao  noil 

F 

display  Jig-  and  uLndous 

subplot  M,.  1,1:] 

plot  (fh^L^firMral^*i  PSD  ;  :  rimJ  11^  ,  1  :  '  [ 
bold  on 

plot  (fu/lofi,,31g-nalPSDC:  Trm,  11 1 J  ’  .  '  j 

1  plot  t  Cn,  |i-  cnPED  ^  i  j  rmH  111  ,  '  b-  1 1 
set(gcflj  '  For,i  =  iLer  H  [?l  3 
.  *t,c  t/u*,  '  VTiSKUfc*! 1 «  ri  1 
title  !.  'Rk  and  Hz  PSPs' | 

j  KoriCAU  [  1 1  i  mi=*v[2]  j  ™nj*~1'"'  P'l  J  vnx-ai.'Nj 
axia(|D  iU-J  niji^ir  juaJ  ynx] } 
hold,  on. 

p  l  Signal  ;  F:x  ;■  and  noise  P3Da 

f 

t 

jr 

subplot i 4,  lr 3| 

p  L  oc  ( f  n/ 1*  f.  r  31^-na  3  J'ShC  j  ,  m...  \  J  |  Hoj.*a?&h  ( : , »|) 

■  plot  C  f  n ,  PzrnPED  ( i ,  rmH  111  !  1  b-  1 1 
seLCqoaj  '  For.-slz*-'  ,,  [P]  ] 

'  set  Cgc*. 1 TTiokhobel 1 .  M 1 
title  1. '  Signal  PSD '  J 

av-axis  j  xnn-avni  i  *rn5c-av(21  j  yran^aii  [3]  j  ynx^aiu  |-4 1 
axia-UG  £00  niji (Or>nn?  ynx]  ? 
hold  on. 

r  Signal  FED 

t 

I 

99 


2 3:1?  EH  C:SRQS\DataAnfll.yaj.g\Ri0.5JrDira\aD-DD  SC\ EftTREDia:  HD  K3  vEep.n  7  d i  11 


subplot \4 r 1 r i )  t  Normalized;  PSD 

Ipl^t  (iB/l«3hr rtSifthPSDf  ;  r tiJii  LLI* 1  ^  ; 

set  (qcaj  1  Fonts  lie r ,  |9|J  *' 

sat  (g^i ,  1  V7 IckUb*  L 1 .  r  ]  I 

title  [ 'Hemal ized  F-5T  I 

av-aiis  j  xnn=av[l|  t  hu=iv(2]  j  ynn=av  [21  j  yni=av  N)  ? 

2  DO  plh  <0  r  i'nrt  !■  tfFiftl  i  ; 

hold  on  7 

end  i 


rpiltltf  tr  IntOt  cjf#t>  to  ftMKfciBU*  ^h'  I  r 

query- input  I  1  Enter  <a>  to  stop  displays  r  1  P  “  3= 1 1  ? 

dillOtf:  [1  itfliM*)  f 

i f  1 l-query-™ "  G '  I  I  Cquery"™  1  ?  1 1  > ;  breahj  end  i 


end 


. . . 

iiap I'  ')  f 

dJLep  i 1  5  t^rtlng  CPP  ^ilculati-on  loops1  5  r 

dlapr  ->  r 

l  ***  CHRtl  *** 

FH  -  zeros INseq, Hiinesh  ? 

HU  -  z£-i  ±L  i,Miu  i|;r  Hlln.* ji|^  r 

CHEtl  "  zeros  INeeqr Hlinesl  r 

1  1 1  f-  l>::  fn|  |;  j  1  |  J  si  -|:j:ira'  hi:. I  rnijH-n  U[r.',ie:U»< 

Cor  tui  - 1 :  M j-L-q  ; 

for  ll-l;Nlines  ; 

DCs  iq  ionH  11]  “  nean  (data.  (Pt_jsi :  Nt_s2  a  ttji ,  11]  >  f 

Kvwiwulll  -  Ihwn.wnwwt_risnt_ra.iw,  11) >  ? 

end  ; 

on  d  ] 

l  Cutout:*  return  pouts,  nolac  powo c ,  |  signal  powar)  P  and  CbfBtl 

lor  tm-l  :Hseq  f 

for  IWiNllnoa  ; 

RFtl  lim.  111  - 

nun, |  |dato  |Ht_sl:Nt_s2jiinrll|  -DCsluijovr.,^!]  j  .  *  |daca  [Ht_s  1 :  Hc._s2 ,1131,11 I  -pCilgtnfi,  11 1 J )  / 
mil  Innp  111  - 

sun  I  |  da  to  I  Ht_r  1 :  Ht_c2j  im,  11J  -  DCnz  <mT,  11] )  .  ■  I  da  ta  I  Hc_r  1 :  Ht_c2 ,  rai,  1 1 1  -  DCnz  ]amH  11  1  J  ? ; 
■CWJU .  jam,  1 1 1  -  |Rkcl  irsn,  11J/H  <&TL_a2  Mr.  iltlM  ^’c_f2- -Kt  rl>l )  ^  *|fflt  1  fra*,  U M  1  1  ] 

end  r 

end 
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J:l7  EH  C:SRQS\DataAnfll.yaig\R05grDira\aD-DD  3C\ DATFiEDVC  HD  K3  vRep.rr  a  D i  11 


4  .,r  *** 

ZOI  CiJi- 1 J  Wifir-i; 

Cor  11-1: Alines 

RHf.l|m#  11[  =  jun  Signal  rani  i  rnn.,  11 JJ 

KH  C  3  |n%  III  -  Siffli  P  m  r  1 1  S  > 

‘CWH£Hnn,lI.I  -  PRi  1  frrn  ,  11  >  /HH±  1  (nn..  1  i )  -1 
and 

end 

i 

r 

i 

**  CMPk  J  -  P-oint-b'/-pni.Tie  Crequency  nacnttlisat ion 


*  ExcesoNi  Factor  -  Q.JJ05 

l  ExcesaHzFactor  -  O-TIQ 

EpcaaflNiFftet&T  -  fl  ISL< 

SHE  “  surH'PrciiFED-l-EHces.sHi  Factor! 

5uE  cm  =  2 :  Hseq 

Cor  11  -  IrXlinea 

UHRk.  3  lain,  11 1  -  (1/  1  HFFT/Z^-l  J } k  EUR  ( 1 ,  mi,  1 1 1 
and 

end 

t  appro jl  CHExI  random  noine  di'i  factor 

Ft  approx  CHExZ  random  noise  d±v  Sector 
r  '  ,jp|J  :  :>:■!  CMP.!-!  3  t  |  vising  djU  frt-r.:.! 

? 

F  *  scaie^  Carter  Is  inly  approximate 

1  ***  3L,splay  CKRa  "* 

— - - — - — - - — ■ — — - - - — - - — - 

Snoot h  CXXe  -o  ceduoe  effort?  of  a-ingle  point  nonlineanes  rn  -scale  Carters 


TiW»t*  -10 

r  '  iin  ii  i"  1  :  ■..  |  hxI  -■  i !  •  H--  1 

Cor  EITi-1  J  Hitf  -if 

Cor  ll-l;Alines 

iui=  [hh-IJ  *NIinea-*-ll 

c^i. )  ^rtfr.P  u )  -  cnn/  P m  3-e-n.n 

end 

and 

;  1  craata  pulsa  zljiia  for  Hanning  TSt^-ilIi 

F 

t  \  WHKiMI#  PJLW  UM  Mill.  1 1 

t 

Rda-a  =  CH Bttl 

Re  n  n  i  Jig  _7  Rnoc  t  h 
qq_tl  ■a  Hdata 

i 

f 

Rdaza.  -  CHRfl 

Hanoi  nq  7  Smooth 
q.q._Tl  -  LI  data 

t 

Rdaca  ■  CHRx3 

Hanni  nq_7Enooth 
qq_0  =  Hdrta 

t 

101 


J:l7  FH  C:\j™\Dttt^^alyaiB!VBOaPrmra\H[]-DD  ncYD&TKEDIK 

HD  S3  vEep.m 

S  of  ll 

PcIIltl  =  meonfi3q_Ct  .^qq_tl) 

7 

5s;  1 _ C*t  1  “  nwfln-S-5^_0-/-n_ti ! 

' 

-o  r  11=1  j  Nil  lies 

t 

I: i  IT ure  |  'Position1  ,  L£^-*-2  5*  [11-1  ? - 

sol , YO-2 D  *  Ui- 1 ?  *  e Cl H SOD* ac2 , 

1-gO*Bc2| ) 

? 

subplot  H,  ir  11 

1  Scaled  DTCRs 

over  le/ 

plot  (CHRcl  (: ,  LI)  H  '  V  1 

I 

hdld  :>h 

plot  i'CNRI  It;,  IIJ  ,  /Ecl_£  ]U(lr  I  L.I 

,  ‘fl-'J  F 

plot  (cmrm.3  (j.115  -  /sci  att  <  l  r  in  , 

T")  f 

set  (go*.  'Jdr.tSize1  F  L 5 1 1 

; 

It  set  (gcflj  1  K7.icle  Label'  t  []  1 

r 

title  t  'Qve.rLa.jied  Scaled  CHPeF  ) 

av-axis  ;  manr-av  ( 1 1  f  >:n>:“av  (21 

j  yim-av  [3]  ;  yn.T"*v  [4-1  p 

™m.i o  c  f acxirt  jhi*  mtarQjymnt  yniwll 

? 

1  hold  on 

subplot  HP 1  r 2~  1 

; 

*  DNhtl 

plot  (CHE:t  1  [ :  j  11) ) 

r 

£wr  1  S-  ■:  ■  '  r  T^l  5 

t 

‘i  set  (qca.,  1  YTictLabei1 ,  U  1 

p 

title | 1 CNRtl ' | 

F 

j  Jinn-wvlll  ;  ;rtraS-4V  [21 

;  Vnn-SL'  1 3 1  ■=  j'pjv-sjI  \4  1  r 

axis  (Iran  *nx  min  1 D.,  yuci  |  vtt-ik  J 1 

7 

1  hold  on 

r 

subplot  M,.  1,3] 

l 

1  CKR21 

pl-Gt  (CNMl  C  i  p  1 1  :> !« 

set  (qcaj  'Font Size 11 P  |S|.I 

i 

1  sot  (go  j,  1  YTlckLab*! r ,  [  ]  | 

t 

title  1 'CNEfl 1 1 

r 

a.-axia  i  >:nn-avi;i|  i  xnx=a,v  [2] 

j  yriQ*a.v  [3]  t  ynx=-avi4J  p 

axia(fK»n  >i»x  ro,  yrwi  l  yMMli 

r 

1  hold  on 

? 

subplot 1 1 , 1 r  4  1 

1  CNRk j 

plot  (CHF:>3  ( : i  11) ) 

7 

set  (go*,.  'ronf;  .-■?l  F  [  P|  1 

; 

t  set  {geo ,  1  YTlokLabel  %  [  ]  j- 

7 

tLtld  [ '  tmxl '  | 

i 

flwni!  j  ann-avl II  j  ^nx-avlll 

;  YTan“avl3l  ;  yns-avlfl  r 

o:ti  s  i;  |  urn  v.TLf.  min  [Oj  yuan  |  ynu<]| 

F 

1  hold  on 

P 

i 
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■lisp i "  Done  MS4  calculation  Loop"  |  j 
dispT  ^>; 


I  '**  Save  CHEr  **■ 

*«ve  iioet  cti|i%1  i“KHfj  r.m*.n:\ 

■lisps'  KIDD  .  :lLi  L  ^  JMdil'  [  r 
dlap  i  *  '  ;i ; 

T*hd  *  cputln»-T*t*it  ; 

1  7  tad 

f  print  £  ( *  1  Enter  meq  I  to  dlnp-lo?  \n ^  k  ; 

q’jeryrinpnt  ( "  Enter  **[]>  to  atop  2  '  |  ? 

delete  i!l;  Klines  I  r 

i  delete  1 1]  p 


1  Flc-ts  fee  acLe.Lod  sequences 

wtiilelquery-Tti  j  * 

nn  -  query  ; 

TfO  -  400  ; 

scl  w  i  j 

SC2  »  %  ■ 

if IKlines>-30|  p 

YD  =  TDD  i 

KJl-1/2  ? 

SC2-2/1  r 

and  j 

tna  -  time  f 1 1  p 

T.lfifi  -  4*11  (t^S  P 


%  fpiiiitfi1  Stgueiice  a;  j  Xn'jimJ 
for  11-1= Klines 


EL:] ut'a  I  1  p  [251-25+ rkl  l!-*aclrVCi.^D*  CU  1J 'aoL,  BOO^iie:  ,  TDfl+flcSn 

fig-uce  rPceition  1 ,  f £&+25*  <11-31  mazl,  10-25'  <11-11  " ncl,  @DG“5c2  ,  30fl'5c2  j  ^  . ,  . 

1  Karre '  j  LI no Name  rLlae  Index  ■;  1 ,  ll]  ri  :  |  I  t 

I*XG  [rtjifc- [i  ,  IS*  <xirji-fimn  1* ,  ,  OS*  (ynii-ymnl  ,  ... 

LlnelilBirie[L±BeIndex<lPil|l ,  i  |  r  ►Color1'*  [Q  1  1J*  p 


subplot lr  n 

plot  (tine*  lefij  -lata  1 1  ,  mri,  1 1 )  ^ 

«Hif  (qtri  ,  1  I  :  '  ,  I  0 1  ) 


Tine  aeries 
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sot  Cgc^j  '  tt lck Label '  d  [  ]  |  i 

tit-1*  L  'T  in?  5er  ey  1  i  t 

airjuin  ;  yin-av f  1 1  ;  xnx"-av(2l  j  ynn-av  r 3.1  ;  ynx-av  I 4  i 
asclaUQ  3D  -1-a  1-01)  i 

axiaUO  100  -L.0  1.0] S 

hDld  DJI  t 

1>ldG.Sifl&*ll 


display  sifl1  arid  r*f  uind-yu* 


^■Lbjl-lul  | 4  r  1  r  £j 

plot  ll  En/  le  6  r  5i-jn-a  jlFED  ( i ,  n3iH  111.  "  ~  ’  I 
hold  on 

pi  Cat  (  f  h/  l*  S  r  NpI  PS[]  (  :  r  nn,,  1 1  K  1  :  '[ 

plot  ( tn ,  l*c  mP5D  { ,  rniH  11 1  ■  1 1—  1 1 

soz  (gca,  '  r  ..•  .  -  "  ,,  [  ■&  5 ) 

set  (gee,,  1  rXickL-afcei1  j  []  I 

title!  'A*  and  Hz  P EL>sm  \ 

fly-axis  j  xnn-BYllI  ;  Knx^av  [$-J  ;  ynj-i^avl3l  ;  ynx-a'v  I A I 

axis  MO  -  LI  >v  nln  (0j  yum  |  ymn]l 
**1  a  ( [0  500  b  1-  n  (0  r i*d )  i'bkl ) 
hold:  on 


SLy  [rctfj  arid  noli#  P£i^& 


) 


subrplo  t  M  P  1  r  3J 

plot  ( fn ■' le E ,  Eigne  1FBD r : , nnH  11 1  -HoiseP&D r : , xdih  1 1 1  I 

P  L  ot  i  f  1 1 ,  tii  rwF3D  C  !  r  hu ,  1 1 1  ,  1  b  '  | 

set  (gca.,  1  Font  Eire r ,  [9]) 
set  Cgc  a ,  '  i"T  ickLabe  1 \  l]\ 
t  Ltle  I  lrqhtl  K-U T  > 

ay-ra*:i5  r  >nn-fl.v(ll  j  xrix-av  1.2 1  i  ynn-av  [3j  ;  ynx-av  Ml 

axla  i;  1 0  IDTih*  mlii  1 0.  ypn  I  yrriM  1 1 
axis  MG  SOQ  numCOryJin)  ynx] ) 
hold  on 


?  1  Signal  FETP 

F 


I 


subplot  id, 1,  4  I 

pi  ot  ( fn/  U  E  p  HoiijsIhSD  M  P  nn ,  1 1 )  !■ 
set  (go*.  'rPr.tSise'  ,  LB]  3 
!■  aet  Cgca  ,  1  liTickLatel'  #  []  | 
tifc.If  £  r(feXMllI*d  MG1  } 

a  van  is  ;  xnn-av  ( 1 1  :  xrrx-av  1.2 1  ;  ynn-av  1.3.1  ;  ynx-av  I  -1 1 

1  axlaMO  -Erio-»  min  1 0,  yum  |  ymx]| 
axi  s  MO  tO-J  nlji (0 r ynn)  ynx | ;■ 
hold,  on 


r  ■  Normalized  FED 

? 


J 


end 

di  &p  t '  r  b  ; 

fprlutr  ['  Enter  soq  1  to  display  \p*  |  ) 

quecy-ihpUt  f  1  Enter  '0>  sl-dp  'I  ; 

delete  11 iKlinea)  ? 
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C3.  HDPulseSim..m  -  Heterodyne  Pulse  Simulation  (Optical  Field  Based) 


| 

<  Hat  lab  program  HD?uL  selim.  rr 

1 

Frograji  to  B-inrulate  rrodelcche-d  laser  pul  s  e  in  both  tng^ral  and  fceguer.cy  spare 

1 

*  Cleared  14  hijc  Q3  by  d-C.  -iSHfifG 
i 

1  .Hod  l  fleer  ion  * 

*  dd  nun  yy  -  hihk. 


clear 

pao!t 

warning 

t  Ua*tl:&  -  'UWS  1 

I  DserlE  -  '&L*PF  1 

I  Din  IS  =  1 "?  S53CI  F_PC  1 
I  UeezlZ  -  'SenftJ^FC' 

□serlS  -  1 3er.it_LT  1 

jtKT  “  conptter 
platfcm  =  jew:  [li«4 ) 

<:le*r  fi 

Land-;'  rsojiJ lOO^cloclj  j 

candnC 1  e;  ■  ■  : !  r  s™ {-  lOO1'  clott]  > 


’  Constanta  I  physical  and  system 


if 


?  I  Uba  wily  cut  5.  a 


r 

r  1  Initialize  random  (  ganaratir 

?  *  initialize  Gaussian  rand, on  I  generator 


qe_det 

- 

0.6^ 

r  * 

detector  quant  urT  effios.en.cy 

qq 

» 

1.6e-13 

-  ? 

charge  -t  electron 

hh 

<* 

i  ■% 

Plesficif1^  cnn*t4i!nt 

capeed 

- 

2.&99eS 

r 

speed  of  light  in/sj 

wvl 

- 

IQ. 5yie-b 

r 

wavelength  IC12-01C-  10PZQ1 

Tran 

= 

2e-S 

r 

DAG  sampling  Interval 

PW_DPj“ 

- 

1/  ^TSoinJ 

r 

data  collection  bandwidth. 

12. 

V 

3-3 

r 

cavity  length  fn| 

d.Mir_ir 

3  3. 7rt-6 

: 

rieq  oftflHt  tmeuoun.  Pjc  end.  |.u 

iotuaLly  rancM&h 


■  Create  m^eloohed  temporal  pulae 


IT  ^  3  Da  6 

dt  -  i>.£'le-9- 

Nt  =  TT/dt 

l>t  -  le-S 


■  din*  window 
?  I  true  point  spacing 
?  -  nunber  of  tenporal  points 

r  ■•  .yyr««  H  rirrz:  ]  i  hrj 
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BfiS/ 04 

3:lfl  FH  C : \ LUOS (Mode is\ He  ter oivne 

Fuise  EiirjulationiUD.  . -.  'iUDF-jiseSin ti  £  of  10 

N_D-_ 

=  TT/Jft 

j  l  actual  nunl-ar  of  sacple  points 

tt 

-  zecos-jln^tl 

i 

ttl 

«  z*  toi!lr  Mt,  | 

; 

tt_Dt 

-  zecoc4lrM_Dtl 

; 

99 

=  zeros! lrHtl 

J 

aa 

-  EOCO*!lrHtl 

; 

FF 

-  zecosjlpHtl 

p 

si 

a  1 

J 

aigi 

-  0.4e-6 

; 

wrdtSti  para rr.  of  farct  F.ici,an 

G2 

■“i 

m 

EJ 

II 

J 

#i%2 

-  2.fl#-6 

; 

parafi  or  Pi  elan 

aiqe 

-  i,30e-9 

j 

Hldtli  of  Gaussian  reaporal  modes  1  -  3  na  mm? 

dr.  a 

=  Z+LXycepoad 

J  1 

spacing  batuaon,  Cause lan  temporal  node  a 

M  t,_a 

-  dt_a/dt+I 

' 

number  of  point  a  an  a  ingle-  interva  l 

for  ii- 

i  i  H_Dt 

j 

Actual  time  o ampler 

tc_Dt  [ilj-  (li-lj  *Dt 

- 

shift  to  stare  at  "  na  to  match  real  An ta 


Ol'id 

-or  Ll=l:Hr 

tel  (HI  ^tt  (ii.J  -7e-6 

■iq  tii l "0 . 0  ; 

else-  r 

7-3  Eii!l  -31*  ltd  Eli  I  /le-*l  t*jkp  I  (-(U  Ul^Kl  Eiil  |  /  (2-aiql'sLglJ  > 
GZ k  1  tt  1  l  it 3  /  le-6)  '■eip  ( ( -ttl  <  i  i  l  ■  tt  1  ( li  1 .1  /  (2  -  aiq2  * slgZ  ^  | 
and  | 

end 

lot  ii-liWt 

tta-nwd  1 1 1  [aaj  ,dt_a]  -dt_aVZ 
aa(ii  h=a*pi;  [-tia.*tta]  /  (2 k e-lga ■  sl-ga \  j 
end 


r:iH4ln  h  |  .ir»i.:  :n  ti  f  Nh  .■  ■!» 


E-f  -  -as- -ha 


pul**  ttJ.in  vitft  nar.|ieKfi]  *iw*lw 

pp  ia  Inetantamecue  -p^wer*  square  of  field  uu 


£  igur*(  ffMitlflll'  J  [50  SDD  130  4-9Q  | )  j 

plot <tC/ie-£rBG!  r 

.hold  or.  ; 

fi  lur  j  1 1  /  j  a  .Er»d  '  ■  1  J  ; 

plot <tt/le-firppj  ' t  ' ? 

ylabel  ( r Optical  rower  Canpcnents  iHodelocfccd;  INi  '  )  j 

ylflibel  (  (>ih:  ■“ )  ; 

ylabel (' Optical  Power  I  M  l  1  ]  r 

ai'-as-is  ;  KinmaavCi^  ;  hit* -air (2|-  j  yrw^J-'U' ( 3 1>  |  yiuc-av(4|-  ■ 
anaE([0  ?r/le-f  ymn  ynn]  J 

Hold  off  j 
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i^£/0^^^1^^J^^£^LHGfiVHDdeis\HeLerDd^ri^^ii^^^iinjLdtion\KD^^\FD^i5je5ir?-T^^^^^^3^C^ 

'  GcGgHdl-e  IccT-enc^  optiPtiitiD  fnr  mod  bItcH  *4  tnriiorol  pul^f* 

Nirnujoa  ■  vmin  Ufgg  mt>  k  ; 

H ITT  -  S'llFFTbaM  ; 

Fp  =  uecisU#  NFfT)  ? 

Fpp  -  E4-E4»11PHZn)  ; 

ire q  -  zeros'll*  NF=T3  ; 

ppl  c  SSfuaUrHfFTj  ;' 

ttl  -  EecocU.NrFr'!  ; 

r&  -  mrppMuriFm  ; 

Fpp  -  f£t^bift(Ipj  r 

df  =  w  (feirpr*dt|  ? 

EH  -  l/[£-dtS 

for  j ]™I ; MFFT  ; 

Treq  j  ]  =  (jj  -HFFT  /2-1 J  *dt  ? 

ebd  j 

ClgunTPultlon'i  [IQO  5DiI  33D  S001 3  p 

AfcibpLftt:  £3,,  L,  1^  ; 

plot.  ■!£  ceq/leE,- obs  (Fppj  i  ; 

tstl*( 1  :^  r.  .  :Al  F\iu&:  PiK^iuncy  S£*ot  \t  iim  Fp/.  f  ’  1  ; 
xlabel  £ 1  Frequency  iFIH  z  I  r  I  ; 

’/label  i  1  Ft'eqycfiziy1  Bpectmn1  ]  p 

4ubpl-&E  l3#  lP£t  ; 

plot  |  1  req/le  0 ,  real  [Fpp3 J  ; 

tit Ji# ( rapt4eq.l  row* '  Fi^queocy  Spec: i ait.  in.a-1  iFppJ  t m  J  r 

xlabel  £ 1  Frequency  iMH  e  !  '  I  r 

ylobel £' Frequency  &pectrLm' 3  p 

asbpLot  i 

plot  |  £  req/le E ,  map  (Fpp'i  J  r 

title  ( 1  dpt  seal  Pout:  Fi  .equaocy  ^eactrurc.  |Imaq  IFppi  p 1  }  t 


jilabel  £ 1  Frequency  <|E>9Kzf  I  '  I  ? 

ylabel  ( 1  Frequency  Epeotnai1  }  p 

ppl  -  if  ft I  Fp'  HFFTH  df  I  ; 

ior  11-IjHFFT  } 

ttl (ill  - < i i-ij / £HFFT*di  |  ; 

■jJid  f 

figure- £ 1  Position"  t  3G  S10  4  D3]|  p 

plot  i  t  tl/l,e-  E  r  ppl  J  ; 

hold  or. 

plot ^tt/l*-&r ppj  'g-s'J  i 

tl-tl-fS  pW:ieii  i:rim  ihi  1  j  ; 

xlabel  £ 'Tiree  t,ust ■  J  , 

y label  ( 'Opr. : oil  Pot-. .  .-:  |H|  1  |  p 


sv-axia  r  Knui-av  «j  1  p  r  i-jra-i-av  ■;  2  }  r  /nn-dv  ■;  $  h  r  yntf-tr?  ( ■i  h  r 
axis|[D  "E/le-fi  ymm  ymx] ^  r 
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1  CGnpute  Btfal*  I-adCor  Lo  trattj.4liiq  to  creeled  energy  4  Integra ted  inctah1 Eaneowe  power) 


F!fc':|.i  -  fiUrtlpft}  *dc 


I - - — - - - — 

l  CUMlO  :.3HH.>K-iTn  Ln  w.jivi"  r  :  -rn-. 


ppjbfi 

-  Zaras  ral  za  i pp'i  j 

pp_L^Jac 

-  zero-a  Isize  Ippl  J 

Fp_LQ 

=  zerc-s  [size  (Fp) ) 

rpp_w 

-  E*  !■:>*■  £#t 1*  £fpp)  J 

ppI_LO 

■  zeros  1  sue  Ipplj  ;■ 

Fp  LOae 

-  zeros  [ai za  ( Fpj  i 

rpp_LjV-ao 

"  zer^a Isize llppj i 

ppI_tjCiac 

-  zero  s | si se (ppl ) ) 

f 

7 


i 


di=p  r  '  J  ? 

All  -  Input*]  1  Srith  :  ,„t.  si}Vj|-  a-m  \v  Ta  ::T.ui  [  Ll ,  Sft-*-  3 1  ;  1  J  • 

iliA  -  input 4  *  Enter  2nd  apur  scale  factor  ll.Se-JJ:  1  1  f 


i^hiT?  "  a^pi'rpnd 
ptul"J  '  2*p.ihrajid 

1  s_p.?  ■-  1,3a  .1 3 

delf_LO  -  77e6 
l  aig-LO  =  E^j 


;  *  arbitrarily  chbeen  phase 
i  arbitrarily  chosen  phase 
jl  no  la*  }*v*  1  For  l.u  ^pa^tinri 
;  I  spacing  between  Lb  main  mode  ami  apura  (Hz  I 
i  \  spectral  width  oC  ID  Frequency  redes  ■. H z i 


for  1L”  1 ,:  Mt  r 

pp_LQ£Ll,l  3  5  +  S7T+ tx;.£  1 2*pl  *dal  1  J^l*  Ct  ( II  |  t^h  ITT)  +  ... 

e-154'cob  £2'pi"  (2Tdede_WJJ  'tt  Iiij4pbii54l  ; 

pp_LOac  1 1 1  >  =  sll*  cos  [  2*pl  *del  l_LO*  tt  ( ii]  *-p!n  iTT)  - 

# 1 54- gm  tz*Pi *  iz* m l e_K3) #  tt  4 1 1 j  +ph ii5+i  f 

end  r 


Fp_l*5  “  f  ft  (pp_2/3^dt  i  NFET I  F 

Fpp_LO  =  EftS^-lFt  4Fp_L-DJ  r 

^2_E^  “  irrc 

Fp.Lbac  =  Eft  (pp  LOac*dc.PSfFTri  ? 

Fpp_LQac  "  fftahrft ^Elp_LOac?  r 

ppZLOac  =  iff  t  (Fp  LjCCiC  ■  HFF7  ■’dfj  i 

figure ( '  Position"  ,  LTOO  3D4I  3SD  E0D|1  f 

*ubpl*t £2,  1, iy  i 

plot  4  f  ceq/le  £  r  abe  ( Fpp_  L/T4  I  f 

>iu  Id  on  j 

j;lo"  ( I  ceq/le6rabe  ,.  'r-'j  f 


t  ltl  e  ( 1  LjO  Frequency  Spectzure  [Fpp_LO|  '  J  ; 

xlabal  ( 1  Praquancy  43iKi  I-  ‘  |  t 

Ylabel { 1  Frequency  Spectrin1  J  ? 

subplot (2j 1,2  £  ? 

Plot  4fteg/l*firr&el  £1tep_[£?ii  f 
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i^£/0^^^1^^J^^£^LHafihMDdeis\HeLerDd^ri^^ii^^^iinjLation\tiD^^\ED^isje5irf-T^^^^^^3^C^ 


Lriuld  on  i 

pin'  ’ifceq/leGr  resl  Lfpp_LCfld  ,  't-'f  ^ 

titles1^  Frequency  Speccrurt  IReal  |FppJi  [ '  f 
Ulibui  (  '  i  :\:Lji.,  ::....  «J*IH .. :  ‘  |  ! 

y  label  ( 1  Frequency  Epectr™1  1  ; 

S  Jbpll-t  [3j  1,.  3h  f 

pldT.  (fteq/lsGr  ineg  Lfpp_Ld  I  i 

hold  on, 

senl  1  -sq ■-.'  |  f  i«q/ 1  c  €  r  in  as  >\  Fpp  ~Jj*  i  |  ,  1  l  -  1  ] 

titlet 1 i£>  Fre-suency  Spectrum  llmag  IFppv  !  '  \ 

jc  label  (' Frequency  (MU  z  |  '  |  r 

X  lib*  l  (  1  PrmijUH  ncy  Ufi*n  1  ]  t 

siguifi  ( ■  r.:.  :-i.^n ,  [70Q  sn  aso  flaon  j 

plot  i  t 1 1/ le-  6  r  PPl_LO  I  .■ 

hold  on  j 


plot<tt/le>-Srpp_LOr  1  1  h 

plo^ttl/le-Ei  ppI_LOacH  'r-'J  j 

placftt/l*  6rpp_U^e,  'y- '  i  f 

titlefLO  Optical  Bower  HIFF7I  (tfh fe  1 
jc  label  ( 'Tine  ■!uj['j  : 

ylabel  ( 1  L£?  DpL  i-v-s  1  Fs^cr  "  J  ; 

av^ascis  r  xmei-av^lji  r  xjnx^av(E>  r  yn*i*av(3|  t  ymK*iv(4|'  r 
7tVle*-£  ymn  yiiH.1  >  i 


•  Create  H >:  end  L£?  npf.  iP«l  field  ¥«r sidles  fiPm  upti?*i  p^nr  i  erns 


■uu  i 

a 

aqrx  |pp) 

■uu_l/j 

“ 

aqct IppJLOl 

tv_r 

- 

nr,mu_rHi  J  Nfc'rn 

FUu_r 

" 

f  ftahift  (ru_rj 

Pj_1a3 

- 

fft  iuu_LC-dt,yo7 

F\ju_lsQ 

= 

1  f  ratal  i  t  (Fu_U^ 

1  Display  field  /  power  trapo-i/al  f  spc-cciraL  values 


diipT  ')  ; 

disp  i !  If#  Fee  S3  return  to  continue  ill'  )  ; 

d lap  r  1 3  / 

pOUiHi  ; 

disp^1  Displaying  field  ■■'  power  temporal  /  spent  cal  valuer  1  S 

dlapr  1 J  4 

de2et,e(l:SJ 


figure  ij 1  Fos  :■:  i-E^n 11 H  I  -55  60*0  S5G  -800 1 1 

ploHtt./ie-e,pp3 

Jtald  nr 


s  in  pp  and  llu_c  teu^orel  display 
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i/2£/0^^^1^^J^^£^LHGShMDdeis\HeLerDd^ri^^ii^^^iinjiation\tiD^^\FD^lse5irf-T^^^^^^3^C^ 


plot  >[  t  t/lm- &,  uu_s ,  '  c : "  h  J 

titli#5  1  v-r  I.  Ji-yi  VoU*Z  illl  1  i  ; 

k label  ( 'Tine  ^iisi'J  ; 

ylab#l  £  'Optical  E'Cik'.  .  [U  |  1  |  j 


figure  ( ■  1  .... -_ izn- ,  [  10  30  S5C  400 1  ■ 

yybtil-S-C  ££,  L,  ]> 
plot  >j  f  ceq/le  <•  r  aba  ( FppJ  J 

Jiuld  OP: 

plot  ( fceqV  lefirabn  ( ruu_r?  r  ‘r;  ') 
tatl*( ' Frequency  Spectrum  (Power  and  Fields 
Kl^t^-1  (  1  |  r  s*i jijh  in:  -j  |tf!!*|  ‘  j 
ylabel  ( '  Frequency  Spectrum'  I 
subplot  [2 j  1  r  2 1 
plot i f ceq/le  6  r angle I Fppi i 
no  Id  on, 

plot- <fceq/lefir  angle  iruit_rl-  f  '  :  > 

t  itie  ( 1  Frequency  Spectrum  Fhase  i.Powe  and 

Slab* l  £ 1  FrHLjLm .-i.±y  ^Hxi  ■  | 

ylabel ( 1  Frequency  Spectrum  Fhae-e  1  \ 


,;  FZ  Fp  and.  Fu_z  spectral  display  ilppjju n 


J 

Fp  .  i  p  j  ' r  y  r 


? 


:1  dj  Ft  i  Fli  1  | 


HgU¥4-£ 1  Jfoftilij&ft1  n '[70Q  601)  550  4001)  ;  'F3  pp_Lo  *ftd  UU_i£J  fceflpfttftl  display 

plot  i|tt/ie-Erpp_tJO|i  f 

ho  I  d  op  i 

plot  (  tt/le-tiruu_LC?r  ‘  r  ■ '  I  r 

titic-C  LC  Cp:i:al  Power  and  Field  |Maaqrt<W>f  |  f 

tfldbtil  £ 1  I  <:  ii  h  ■  !  ; 

ylabel ( ‘Or  tecs  1  Fowez  :\nd  Fie,.  '  y  j 


figuie  ( '  Fisizi  zn’ j  [700  30  5bC  flOOl)  ;  MA  Rp_L0  suid  Fu_LG  spectr  diap  | Ppp_I0,  Fuu_LO) 

eMbplOt  L2- 1F ; 

plotf  ireq/lefi  r  abs  (TJip_UOSl  I-  j 

ho  1  d  or.  f 

plot  <fceq/lefiraba  (Fuu_[/?*  ,  1  r :  1  > 

title C'LjC1  Frequency  Spectrim  IPiiwer  and  Field,  j  Fp_La(d  |Fa_3XS|  1  J  ? 

£  1  ^  ^Sl  1  I  r 

ylabel ( 1  Frequency  Spectrum'  J  ; 

suhpl-zt  |2,  lr2^  } 

plot  <  f  ceqy  le  fi  T  angle  i  Fpp_L0 1  1  ? 

hold  on  j 

plot  ^ fceq/le6rangle  IF«u._L0l  *  ' :  ;  ’  )  ; 

t itle( '  LO  Frequency  Spectrun  Phase  ■; Power  and  Field i  Fp_ljC-|  Fu_ljCr  ^  ; 

Klabfrl  ( '■  Fr-ij  .  :: ..  y  <rfH.  ■  ‘  \  \ 

ylabel  (“  Frequency  Spertrmi1  I  ; 

i  Check  thet  spectrum  at  received  optical  pover  is  canvnlutaon  of  P.3t  opt  field  apectruiTs 

icnvflq  -  0  ;■  %  execute  ,£inq 

iC [lcnvClg==l>  r 

iJleipt1  III  fUi  rijiisvej  lUvi<:fl  tOUT  me  #11  1  1  r 
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EH  C:\LTlCffl  hMcdeiB\HeLe:rc-lyrie  Fulze  Eiir^jlatignXHD.  ■ /'■KDF-jjg-e-Eiff-Ti  7  d£  IP 


diL.pi;'  1(1  PECis  icmrn  ca  begin  liJ'  J 
dispE'  ') 

pause 

dlapi;’  Rx  cj-vfivuliit Inn  cljagl  irvjtLrvi  stiXEad'  ] 
di  c-p  l  1  1  > 

Ui  bfslze]  =  sl2«[Puu_r| 
tluw  - 

xFll  »  Fou_e  I  If  <Jf5i.zeA2-Huat/2^I3  ;  C Fa i.at/2  +H1111  e /2 ]. 3 
xFp  -  ffpp(lr  |Hsl«/2-lftiM/a+U  s  |»ilza/a+HDia/an 
xFc  ”  canv[xFuH  IfFUl  ^-zf 
xFc  =■  xFc  ( 1 F  Hose  /2-*-l ;  ^-Jse+Kujse/Z ) 

rUfttctl'i  I-  .  :  r  raoo  ssn  40&1  3 

plot  Eabe-lHFpl  I 
hold  tn. 

pl-ot  iflbslyiFol  .  '  ;  1  J 

end 


:  1  ■  IT  j  Stf| 

t 

r 

;  -iVHrlHy  pi>h*rt t  «*.r<:=  CdnV  **:hi:-  rn 


t 


1  Ct'-eefe  j l  ji|>u * c. hull  4i f  Lii  optical  yOhS/t  1%  ^envcil  unui:  if  L'j  optical  fi*14  opue c^un 


ienvflq  =  G 
i£  [IcJWflg— 1^ 

dls-pf'  111  LG-  convex nil an  ch.-ecfe  routine  '3 

dj.ia.pf1  III  pta* >  -iHr.nri:  co  U-iin  Ilf  '  3 
diepf1  '> 
pause 

dlS-pf1  LP  BHIWLilE  ibh  "hec*.  r<=--i r. L r •-  1 

diapC  '3 

>h  Mil  Zb]  -  £  J  24  |  Flili  LU| 

Huae  -  TSOflO 

xFu  -  Fuu_I£  1 1  f  (HjlaE/’Z-Nuse/i'-l  J  J  EHslie/2+-3tuia/2J  J 
KVp  -  fpp_Lol'l,  EH*  !  (N*iS*/2tHUJn/&)> 

xFr  ■  con.v[xIoF  kFuI  'd!: 
xFc  -  xFeElrK1us«/2^1  i^Jsa+'KIufA.f 2) 

Eig-uce  I  1  S'ositi';:;'  .  L-900  400  55-0  400]  3 
plot  (ohafxFpl  1 

hal 4 

pl  ot  (abs  jxTc| ,  '  r:  ' ) 
end 


r  %  execuce  flag 


t 


I 

r.  *  size  -Di  convoiuEi-En 


t 

r  ;v*rlH7  prnec  and  conv  npec^ra 


r 


1 — - - — — — - — 

*  Clear  Uunaade^  ^arlafe-lee 


clajix  ia  gg  ppl  ppIJLfl  pp[_[£dc  Fpp  Ff>P_W  Fpp_U^e  ju_j  uu_I.d  Fu_t  Fiin_f  fu_U>  Fu-j  bO  ; 


■  Calculate  optical  po*t«r  signal  vn  heterodyne  de'wtar 

diflp  i 1  ill  iidleLiL  jcion  a!  ope  leal  pAuar  signal  hot* radyno  dor*,-,  to?  ip|  '  1  j 

diapf '  ■  ■  ■  Fie 5?!  refjrr  to  begin  1  (*  '  )  s 

dispT  '  J  i 

SfaUs*  j 
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■iisp-;'  Calculation  started"  )  r 

de.leWi-i;.l;4)  i 

I  Revalue  optical  |j f.-ww i 

dispr  ■  5  t 

Eppin  -  inpUtl ‘  K:n»-r  PpCi-^el  «R  dw^tfi'lor  {pdf  E2.QJ:  1  J  ; 

diapT  1  )  ; 

RppLr,  -  EpplJi^Ia-lJ  t  convert  fj.  Dili  pJ  to  J 

pp  =  pp-  |  Eppin /Epp^  r  scale  cacti1. red  optical  power 


I - uy 

optical  poorer  acallnq 

F_L0in  = 
d lap < 1 
F_LDin  - 

■  Input  I '  Enter  LC  optical  fewer  on  detector  « rH l  ■  J  r 

')  ; 

-  P_LGin*le-3  r  convert  from  nH  to  H 

pp_LO 

ppLGae 

-  F_LOj.nl'pp_M  r  ”  acalt  LG  power 

=  P_ljQln*pp-_LOac  i  scale  A£-couplcd  LG  power 

-  Calculate  IjQ  abot  ncj.it  optical  poweo 

(”5 1  iji- 1  h tu  foe  pull  tanpopaj  a-sh  iiii.  thirf.  fl?ata  and  band]  LJijt  to  inatcl.  r.-Ac  h» 


alg_Hi 

=  aqrt  H2*  rhh*c.:;peed/wvl|  *EW_0RC7qo_det> -nean  (pp_LGH  r 

-  Sig_ah+tandn^il5e,!j?p-_l0]  J  r 

Fp  pn 
Fpp_sn 

<=  fit  L'pp  an’dt.MFTTj  t 

-  fftahif  t  rip_wO  ,- 

-  M.Mirid  I  c.v.j^c/d  f  | 
atp_im  -  lerosil^HIFTl 
atp^f-an  tl,  l  iN^BfHQfcCh  =  1 


otp_foi 

lii  Pm-N_lKTOftC+ 1 1 HJFTJ  -  1 

Fp_^rt 

Fpp_an 

-  Fp-J&n - ■  « cp_f M C  (PK/EW_®AC)  i  *  r.u  ir.fi  1  HI- <Hf*  p*u*- 

-  f  ftahift  [Fp_an^  ; 

ppl_aa 

-  i  f  St  1  Fp_an  ‘  VlFn  ■  dl  H  WITT) 

|AA  PLM 
pp_ap. 

-  iicolpp_atil  ; 

«  ppI_EnislHl:BBJ  r 

—  Calculate  heterodyne  optical  p:?wer  on  detector 

r-crforiL  AC-ccupIlng  now  to  avoid  ciracs  Icon  sine  diatortlon  in  miner!  cal  calc 


pp_hd 

"  pp  ■+  pp_LOac  +  2*coa  [2*pi*  rdell_IF'tt >  ^  . *aqrt  |pp.  *pp_LQ^  +■  seal fpp_sn)  ? 

Fp_hd 

Fpp_bd 

-  fft  rpp_JidBdtH^rm  r 

=  rrtsbirt.  (Fpjuij  ? 
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figure ( '  Fuse  t  ion  ’  t  [ICO  51U  400 ( !|  r  Fiq  ppjhd  tensors!  p-lot 

tluT.i|  tt.'j  t-rhrpp  hd./  | 

tit  It  r  KD  Opt  i  “-.  1  Pore  r  IifJ’3  r 

x libel t ' Tint  i  u s t 1 )  j 

ylabftLf  rVFtLe*l  Pfln*#t  iwi  '  i 

\  a'f”-a>:.i3  f  xnntv^ll  ;  xmwav(.2l  ;  ynn-avISI  ;  yrnc^av  ri]  p 
1  ascl  ■  ( j  ]mjL  Jinx  D  ynx  |  j  | 


ligure( '  Pus:. lien"  j  [100  10i]  350  1001)  r  Fiu  pp_Jid  aptctrnn 

#Bbpl*t  12. tr Ik  t 

plot  ^  £  ceq/led  r  aba  £,FJip_JidJ  j’ 

title- ( 1  Frequency  3pe-.it run  (3D  Optical  tauex)  1  I  p 

Jtlabtl  E 1  T  r  . . ■■;  :,  !,i  ■  I  r 

Ylabel (' Frequency  Bpectri^n 1  }  j 

eubp Lot  jC-S j  iLp  £  h 

plot  ij  1  ceq/le  b ,  angle  i  Fpp_hd  |  J  r 

r.  i  T.  l  ir  P  ■  14  H‘.|UH  li  e  M  tJfrt-fi!  L  jrLiii  Pfi4&*  I  HO  -Spt  3.dll  £  t 

a label ( 1  Frequency  !MH  e I ‘  I  j 

ylabc I  ( 1  Frequency  f-pcctruffi  Fhaae  '  !  f 


%  CaiquLaba  eurxauT.  ;A£  coupled  |  *ut  of  douotof  find!  hi  a  a  eunfr  api  fc-.q.  *Ysot  n&i&iri 

I_hd  =  <qa_det.*q;q]  ^  (Jib*cspeed/w¥\lk*pp_Iid  r  ■  includes  10  shoe  noise 

I - Curapute  biaa  current  end  biaa  current  shot  noiae 

diap i *  1 )  r 

I_bii5  =  input  1  1  Intel  bias  current  inA)  ;  .103]:  '  \  j 
di  sp<‘  ')  ; 

I  blaa  -  I  bias*le-3 


siq_Ib  "  aqrt 1 2 T  qq ■ BTW_IWiE ■ I _biae I 

-  $ lg_*b* t ^rji-ji. $ h  i  .'.h ipp_hfl) J 

Fi  b  =  fit U  bias* dtPNfPD 

Fiz_b  -  fitshif  t  i,U-_b) 


convert  from  nl  to  1 
at  dev  of  bi-a-f  cure  shot  noiae 
cytrtitt  *hac  -noi** 


ri_b  -  Fi_b  r '  a  tp_f  ort/adrC-  <BK/M_PA.C| 
Fr  i_b  -  ut3bl.lt  i.E-_b) 


scale  to  maintain  no  tee  power 


nZ  b 


ilftlFi  b ' NFFTJ  df . FFFT ! 


\M  lvb]  -  aiB*  lljbidisi 
l_bl aa  -  i  ±z_b  .j  1 H 1 :  BE) 
1  bias  -  x«aljl  jviaii 


;  luc'Qu  i  bias  !  ■  i  -  real 


l - Corepute  total  AC- rc-jpX  ed  current  out  of  detector 
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I_dec  =  l_hd  +■  l_bias  r 

I - - — 

1  CJlcuiatA  vulijyw  it  Input  till  A/b  radno*  di~*  to  jotual  a  yjittui  aanpllJig  Iritarval 
Detector  - >  K?*&  anp  ->  1  HHn  hiqiti  pass  ->  x£2  aTip  ->  5-G-ohn  Kf 3? 

V_rtDC-  -  l_d*t  *  74  +  1.0  •  i|(  '  5fl 

V  AD  -  intarpl  (tt,  V  ADO,e.e  St  | 

HTFVbaaaV  =  ceil (log2 (N_Dt) ^ 

HTIYV  -  2*HFIT1MI**V 

dfV  -  l/i|CTFTY¥*Dth 

for  zj]-llHFFTV 

rreqV  [}  ]  J  -  ( j  J  -HFET.7  S  - 1 1  h df V 
end 

Fjy_aP  -  ffr.  |y_AP*0t.l*PFTU| 

FW_AD  -  iitshiLl  t  l[F¥_AD? 


f  igoii  [  700  7  D  D  5Sfl  400 |) 

plot i t t_0i / le -f, V_A3 I 
title ( 1  Signal  at  A/D  rv I  1  i 
;•:  label  ( '  i  ,r-  "  h  !  ‘  3 

ylabei (' Signal  at  A/D  iVI 1 J 

I  av-axIs  4  xnri-av  i;]|  4  unsc^av  £21  4  ywn*#T  [ 3]  4  yni-  jv  M 1 
I  axia£|xiui  xnx  0  ynx|;i 

(lOUlWCP^H  i  ‘  j  [  700  HJfl  550  500  |  ) 
subplot  [2, 1?  ifr 
plot^lLaigy/lafi,  aba  (rw  Ab|  | 
title- ( 1  Ftequen.-ry  Epectrum  [A/D  V-aitagieii H  J 
x  label  ( 1  Frequency  :  :-1H  z  \  '  \ 

ylafc+1,  r  i 

subplot (2j lP2h 

p  lui  1  raqV/  le  t ,  uig  1c  |  FW  AD) ) 

title- f 1  Ftequejity  Spectrum  Phase  (A/D  YcltageJ 1  1 
x label ( 1 Frequency  {MHi| ’  | 
ylabei  t, 1  Fteqnejlcy  Spectrin;  Phase1  (■ 

| - — - - - - - - - - - 

i  Save  data  as  file  that  can  be  read  late  DATMGWC_HE  prrqran-  [real  data  amalpn-iat 

a«vdat  “  v_ftp' 

JAVA  ROP-S-tAt  i:-  14  EAVdjt  t 


7  I  ’ .  i  V  PU  :  Ih  !  :  .n 

! 

f 


VFIg  V.  AD  T.OJip&-FAL  plot 
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